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Es wird die laminare Strömung an einer rotierenden ebenen Scheibe berechnet, die gleichzeitig in Richtung 
der Drehachse angeströmt wird. Die Dicke der durch die Drehbewegung mitgeschleppten Schicht-und damit 
auch das Drehmoment sind stark abhängig von dem Verhältnis der Anströmungsgeschwindigkeit zur U mfangs- 
geschwindigkeit. 


Thhe,laminar flow about a rolaling plane disc is computed in the case that the fluid is flowing against 
the disc in axial direction. The thickness of the layer carried along by the rotation depends strongly on the 
ratio of the velocity of undisturbed flow and ihe peripheral velocity of the disc. 

Le courant laminaire & un: disque plat tournant en rotation, lequel est flanque en meme lemps par un 
courani dans la direction de l’axe, est calcule. L’Epaisseur de la couche emport£e par le mouvement tournant 
et, par consöquence, en mEme iemps aussi le moment de rotation d&pendent consid£rablement de la propor- 
tion de la velocit& du courant non relard& & la velocite d& la p£ripherie. 

BiryncaaeTca NAMUHAPHOe TeyeHue y BPalmmalmımeroch TLIOCKOTO AHCKA, OAHOBPeMEHHO 
HonBepralomerocH HATeKAHHIO B HAIPABIIEHHUHM OCH BPAIIICHHSN. TormmuHa 33XBa4eHHOTO BPAIIa- 
Te/IbHBIM MBH3KEHHEM CIIOH, & BMecTe C TEM U: MOMEHT BPAlleHHA, CHIIBHO saBucHT OT OTHO- 
IHeHHA CKOPOCTH HATEKAHHA K CKOPOCTH OKPYKHOCTH UCKA. „ 


I. Einleitung und Aufgabenstellung 


Die Strömungsvorgänge an einem angeströmten Körper, welcher eine Drehbewegung aus- 
führt, spielen bei Strömungsmaschinen aller Art eine wichtige Rolle. Es ist bekannt, daß die 


. Ablösungsvorgänge an den feststehenden Leitschaufeln und den umlaufenden Laufschaufeln sich 


u. U. stark verschieden ausbilden. Die Untersuchungen von H. Himmelskamp [1]an einem um- 
laufenden Propeller ergaben z.B., daß für die Blattschnitte in der Nähe der Nabe die maximalen 
Auftriebswerte sehr viel höher sind als bei geradliniger Bewegung (bzw. Anströmung) des be- 
treffenden Profiles. Die physikalische Ursache hierfür ist jn einem Abschleudern des Grenz- 
schichtmaterials nach außen infolge der Zentrifugalkräfte zusehen. Untersuchungen von C. Wie- 
selsberger [2] über den Luftwiderstand von Rotationskörpern, die sich um ihre Achse drehen 
und gleichzeitig in Richtung der Drehachse angeströmt werden, ergaben eine beträchtliche Zu- 
nahme des Widerstandes mit dem Verhältnis von Umfangsgeschwindigkeit zu Anströmungs- 
geschwindigkeit, und zwar um so mehr, je schlanker der Körper ist. Ähnliche Ergebnisse er- 
hielten $S. Luthander und A. Rydberg [3] aus ihren Versuchen an rotierenden Kugeln, die 
in Richtung der Drehachse angeströmt wurden. Insbesondere wurde von ihnen eine beträcht- 
liche Verschiebung der kritischen Re ynoldsschen Zahl der Kugel mit dem Verhältnis der Um- 
fangsgeschwindigkeit zur Anströmungsgeschwindigkeit beobachtet, Auch hier ist die Ursache 


“in den Vorgängen in der Reibungsschicht zu sehen, wo infolge der Drehbewegung die Flüssig- 
“ keitsteilchen in Wandnähe mit umlaufen und infolgedessen unter der Wirkung einer starken 


Zentrifugalkraft stehen. Es ist klar, daß hierdurch der Ablösungsvorgang und auch der Um- 
schlag laminar-turbulent stark beeinflußt:werden, und deshalb die Drehbewegung einen erheb- 
lichen Einfluß auf den Widerstand des Körpers haben muß. 

Im einzelnen sind aber diese recht komplizierten Strömungsvorgänge, bei denen es sich 
um dreidimensionale Grenzschichtströmungen handelt, bisher sowohl experimentell als auch 


theoretisch wenig untersucht worden. Lediglich der Sonderfall einer in ruhender Flüssigkeit 
_ umlaufenden Scheibe ist schon frühzeitig von Th. v. Kärmän [4] mit grenzschichttheoretischen 
_ Methoden sowohl für laminare als auch für turbulente Strömung behandelt worden. Während 
Th. v. Kärmän ein Näherungsverfahren zur Lösung der Grenzschichtgleichungen benutzte, ist 


% später dieser Fall von als « ng ( 
- Navier-Stokesschen Gleichungen angegeben worden. Das praktisch wichtigste Ergebnis sind 


W. G. Cochran [5] für laminare Strömung auch als exakte Lösung der 


hierbei die Formeln für das Drehmoment der rotierenden re ve Il. re 
zalischi wichtigste Merkmal dieser Strömung, daß die in der Nähe der Scheibe infolge der 
Be cnmene ird, und durch Zuströmung ın 
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axialer Richtung ersetzt wird. Die rotierende Scheibe übt also auf die Flüssigkeit eine,,,Pump- 
wirkung‘ aus. i 

Für.den allgemeinen Fall eines rotierenden und gleichzeitig angeströmten Körpers sind 
einige allgemeine Ansätze von J.M. Burgers [7] angegeben worden, und neuerdings hat L. F. 
Fogarty [8] für ein rotierendes Propellerblatt einige spezielle Lösungen mitgeteilt. Auch in 
einer Note von G.K. Batchelor [13] werden verschiedene Ansätze für rotationssymmetrische 
Lösungen der Navier-Stokesschen Gleichungen mitgeteilt. Dabei wird insbesondere auf eine 
unveröffentlichte Arbeit von D. M. Hannah [14] hingewiesen, die mit unseren nachstehenden 
Untersuchungen verwandt zu sein scheint. 

Wir stellen uns hier die Aufgabe, eine Verallgemeinerung des Problems der rotierenden 
Scheibe zu behandeln, nämlich die Strömung in der Umgebung einer rotierenden Scheibe, die 
gleichzeitig in axialer Richtung angeblasen wird. Dieses Problem stellt den ersten Schritt dar 
zur Berechnung der oben erwähnten Strömung an einem Rotationskörper, der um seine Achse 
rotiert und gleichzeitig axial angeströmt wird: Ersetzt man nämlich für eine kleine Umgebung 
des vorderen Staupunktes den Rotationskörper durch seine Tangentialebene, so erhält man unser 
Problem der angeströmten rotierenden Scheibe. In der vorliegenden Arbeit beschränken wir 
uns auf den Fall laminarer Strömung, aber wir behalten uns vor, als Erweiterung unserer Auf- 
gabe später die Strömung an der angeblasenen rotierenden Scheibe in turbulenter Strömung 
sowie am axial angeblasenen drehenden Rotationskörper zu berechnen. 

Von unserem Problem der angeströmten rotierenden Scheibe bei laminarer Strömung sind 
zwei Grenzfälle bereits bekannt, die sogar exakte Lösungen der Navier-Stokesschen Glei- 
chungen darstellen: 

1. Für den Fall, daß die Anströmungsgeschwindigkeit verschwindet, erhält man die schon 
erwähnte Strömung an der in ruhender Flüssigkeit rotierenden Scheibe nach Th. v. Kärmän 
und W. G. Cochran. Die Dicke der ‚‚mitgerissenen Schicht‘ (= Grenzschicht) ö, ist hierbei 


proportional Vv/w, wo » die kinematische Zähigkeit und » die Winkelgeschwindigkeit bedeutet. 


2. Der andere Grenzfall ergibt sich für verschwindende Drehgeschwindigkeit. Dies liefert 
die bekannte rotationssymmetrische Staupunktströmung mit einer Grenzschicht für die radiale 
Geschwindigkeitskomponente. In diesem Fall ist nach F. Homann [9] die Grenzschichtdicke 


ö, proportional Yv/a, wo a=dU/dr—= U/r den Gradienten der Außenströmung in radialer 
Richtung bedeutet. In unserem allgemeinen Fall der angeblasenen rotierenden Scheibe haben wir 
nach Bild 1 in der Grenzschicht Geschwindigkeitskomponenten sowohl in der radialen als 
auch in der Umfangsrichtung zu erwarten. \ 3 
Die beherrschende dimensionslose Größe für unser 'Problem ist die Größe 
oo or % 


zT ar 7° 

die das Verhältnis der Umfangsgeschwindigkeit der Scheibe », zur radialen Außengeschwindig- 
keit U darstellt. Da die letztere aber der axialen Anströmungsgeschwindigkeit proportional ist, 
ist somit die Dimensionslose w/a dem Verhältnis der Umfangsgeschwindigkeit zur Anströmungs- 
Er, geschwindigkeit proportional, das nach den oben erwähnten Messungen den Strömungsvorgang 


R: beherrscht. In Abschnitt Ve wird gezeigt werden, daß für eine Kreisscheibe mit dem Radius R 
zwischen diesen beiden Dimensionslosen die Beziehung . 


WIEN 


besteht, wo V,= Ro die Umfangsgeschwindigkeit der Scheibe und W die Anströmungsge- 
schwindigkeit bedeutet. * . 


I. Die Grundgleichungen N 


Für die Berechnung der Strömung legen wir Zylinderkoordinaten nach Bild 1 zugrunde, 
Die Scheibe liege in der Ebene z—0, die Drehachse falle mit der z-Achse zusammen. In der 
Ebene der Scheibe seien die Koordinaten r und @. Die Winkelgeschwindigkeit der Scheibe 
sei w. Die Geschwindigkeitskomponenten seien u in radialer Richtung, v in Umfangsrichtung 
und w in Achsrichtung. Ferner bedeute p den Druck und » die kinematische Zähigkeit. Für 


dieses Koordinatensystem lauten die Navier- Stokesschen Gleichungen unter Berücksichti- 
gung der Rotationssymmetrie = =0): i | 
p 
u u dw A 
Fest: Fo Fre 0 (Kontinuität). Bed), 


2 
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ou v2 Mu 1 op u d 2 u) 
FE 7.0 ae FRE 2 Ir Er a ai Def (eadıal)»r.% 7202) 
MV, w v ev  9/v 22V I 
BT: 5 Hw 9% =, OR, =} (azimutal) . . . (@) 

ow 9w 1 op Beh 0. 08 
v dr Ki TER 0 0% ee zB r or ie se (axial) A (4). 

Die Randbedingungen sind: 

z=0 :u=0; v=u-ro , w—0 | (5) 
= 09:4 = U; V— 0 J . 


Die potentialtheoretische Geschwindigkeits- 
verteilung der in Richtung der negativen 2-Achse 
angeströmten Scheibe . (rotationssymmetrische 
Staupunktströmung) ist 


U=ar; V=0; W=-—2az (6), 


wo a eine Konstante bedeutet. 

Von den beiden Grenzfällen der inruhender 
Flüssigkeit rotierenden Scheibe und der Stau- 
punktströmung an der nicht rotierenden Scheibe 
ist bekannt, daß sie ‚‚Grenzschicht-Charakter“ 
besitzen, d. h., daß sich die Wirkung der Reibung 
auf eine dünne wandnahe Schicht beschränkt, 
deren Dicke proportional /» ist. Das Gleiche dür- 
fen wir deshalb auch für Iunseren allgemeinen 
Fall annehmen. Die Druckverteilung erhält man 
somit aus der für die Außenströmung geltenden 


Bernoullischen Gleichung z 
p+5 (U2-E:W%) = cönst. "....... (7) 
zu ; | 
12 
Par; _ =4ar. 
ml 2: Bild 1, Erläuterungsskizze 
Nach den Annahmen der Grenzschicht-Theorie Geschwindigkeiten: 
prägt sich dieser Druck der wandnahen Rei- _ radial: u, U azimutal: v, 7 axial: ı, W 


bungsschicht auf. 
Wie Th. v. Kärmän [5] und F. Homann [9] gezeigt haben, befriedigt der nachfolgende 
Ansatz die Gl. (1) bis (4), sodaß er auch für unseren allgemeineren Fall versucht werden möge: 


»=rflR);  - v=rgß); v—h(e); P=p(l)—5 NS Le, (8). 
‘ Durch Einsetzen von (8) in (1) bis (4) ergibt sich für die Funktionen f, g, h, pı das folgende 
Gleichungssystem: 
2/40 | (Kontinuität) . . . . (9) 
2 2 . 
2 

p-p+rl=et = ES "(@adial) en. (10) 

d d2g in A 
2fg+h2)= 72 (azimutal) : 22... (11) 
dh idp dh ; 12 
h RE de EV FE Fatal ER er ra (12) 


"mit den Randbedingungen: | 
ee > REN. 
Bewer fu, | 
-  Ausden Gleichungen (9), (10), (11) können die Geschwindigkeitsverteilungen f, g, h bestimmt 
werden; Gl. (12) liefert dann nachträglich die Druckverteilung. Für den Sonderfall a — 0 (in 
> ’ { 3 > f 7x 


D- 
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ruhender Flüssigkeit rotierende Scheibe) ist dieses System von Differentialgleichungen von 
W. G. Cochran [6] gelöst worden. In unserem Fall hätten wir (nach Einführung von dimen- 
sionslosen Größen) das System (9) bis (12) für eine Reihe von Werten der Dimensionslosen a/o - 
zu integrieren. Der Rechenaufwand hierfür würde sehr groß sein. Wie bevorzugen deshalb das 
von Th. v. Kärmän [4] angegebene wesentlich bequemere Näherungsverfahren, ‚welches mit 
dem Impulssatz arbeitet. Hierzu muß zunächst der Impulssatz für die radiale und azimutale 
Richtung aufgestellt werden. i 


II. Grundlagen des Näherungsverfahrens (Impulssatz) 

Wir erhalten die beiden Impulsgleichungen für die radiale und die azimutale Richtung 
durch Integration der entsprechenden Bewegungsgleichungen (2) und (3) über z von der Wand 
z=0 bis zur Grenzschichtdicke z<—= ö. Hierbei werden in den Bewegungsgleichungen (2) 
und (3) zunächst noch die ‚‚Grenzschichtvereinfachungen“ vorgenommen, Dies kommt darauf 
hinaus, daß auf der rechten Seite in den { } die beiden ersten Glieder gegenüber dem letzten 
vernachlässigt werden. 

Für die radiale Richtung erhält man dann aus Gl. (2): 


ö 


8 8 8 8 
ou 1 [ ou 1 ([ & 
behede: e u | Ir Ne. 62 9 2: 
IKzZ. 3K + | 0 dz a + |z,|, 


Durch Elimination der Geschwindigkeitskomponente w mit Hilfe der Kontinuitätsgleichung (1), 
nach Einführung des Druckgradienten aus der Bernoulli-Gl. (7) sowie mit 


7 = 
I pl Sr DR a | 
£ 0 ( z=0 \ ) 

ö % für die radiale Wandschubspannung erhält man nach einiger Zwischenrechnung als Impuls- 
Be gleichung für die radiale Richtung 
5 d oo oo oo 
Be ER FR | uU—u)d + 2 (U —u)d+ [az Re (15). 
= 2 ; 0 iE d dr . 
c% k ö 
Br Dabei konnte die obere Integrationsgrenze ö durch w ersetzt werden, da in allen Inte- 
Bi - gralen für z > ö die Integranden verschwinden. Für den Sonderfall, daß die Scheibe nicht an- 
Be, geströmt wird, U =(), geht Gl]. (15) in die schon von v. Kärmän angegebene Gleichung über, 


‚siehe [4], S. 248, Gl. (34). : 
In analoger Weise erhält man für die Umfangsrichtung aus (3): 


8 BD) ö 

2 ö 
[u + [water | v2 ber, E 
; Or { r 02 82 lo 


Durch eine ähnliche Rechnung wie für die radiale Richtung ergibt sich mit 


'n_,(® | 
: =»( =). eg BSR (16) 


u die Schubspannung in Umfangrichtung für die Impulsgleichung der Umfangsrich- 
ung: ; 


ir (u! = — ne N NER (17). 


un Impulsgleichungen (15) und (17) gelten sowohl für laminare als auch für turbulente 

römung. | 

‚. „Durch Einführen der Ansätze (8) für die Geschwindigkeitsverteilung und unter Berück- 

ders von (14) und (16) erhält man für die Impulsgleichung in radialer bzw. azimutaler 
ichtung: 5 


ar : fe +2 I HM... (18), 


a ER | | 
lt toi Ep a LE ETELERE ANRE BE 3125 k - (19). 
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Dies sind zwei Gleichungen zur Bestimmung der Funktionen / und g, die überdies den 
Randbedingungen (13) zu genügen haben!). 


IV. Durchführung der Näherungsrechnung 
Nach dem Näherungsverfahren der Grenzschichtrechnung, wie es zuerst von Th. v. Kär- 
man [4] und K. Pohlhausen [10] angegeben wurde, werden die Gl. (18) und (19) in der Weise 
befriedigt, daß man für die Funktionen f und g geeignete Ansätze macht, welche die wichtigsten 
Randbedingungen erfüllen. Am besten haben sich hierfür Polynome des Wandabstandes be- 
währt, derart, daß die Grenzschicht in einem endlichen Wandabstand z — ö an die reibungslose 
Außenströmung anschließt?). Wir führen den dimensionslosen Wandabstand 


2 
a oe Er EIER SE) 
s (20) 
ein. Dabei lassen wir vorläufig noch offen, wie die Grenzschichtdicke ö von der kinematischen 


: daU E 

Zähigkeit » und den Konstanten & (= Winkelgeschwindigkeit) und «= rn (= Stärke der 
Außenströmung) abhängt. Wie schon in der Einleitung erwähnt wurde, ist für den Grenzfall 
der ruhenden angeströmten Scheibe ö Y»v/a, und für den anderen Grenzfall der in ruhender 
Flüssigkeit umlaufenden Scheibe ö6 — Yv/o. Für unseren allgemeinen Fall haben wir also zu er- 


. warten, daß ö außer von » sowohl von a als auch von abhängt, derart, daß sich diese bei- 
den Grenzfälle wieder ergeben. - 


Unter Einführung des dimensionslosen Wandabstandes nach (20) lassen sich die beiden 
Impulsgleichungen (18) und (19) in der folgenden Form schreiben: 


lose 2. os 
Br: [sa BE En See (19a). 


Bezüglich der Wahl der Polynomansätze für f und g hat sich nach Th. v. Kärmän [4] 
beim Fall der in ruhender Flüssigkeit rotierenden Scheibe ein Polynom 4. Grades für die radiale 
Geschwindigkeitskomponente f und ein Polynom 3. Grades für die Umfangskomponente g be- 
währt. Wir haben diese Wahl auch für unseren allgemeinen Fall als I. Näherung genommen. 
Darüber hinaus haben wir noch eine II. Näherung gerechnet, bei welcher für f ein Polynom 
5. Grades und für g ein Polynom 4. Grades gewählt wurde. Um die Koeffizienten zu ermitteln, 
hat man als Bestimmungsgleichungen einmal die beiden Impulsgleichungen (18a) und (19a) 
und zum anderen die Randbedingungen an der Wand und am Außenrand der Grenzschicht, 
die wir wie folgt schreiben: 


d ö2 ee) 
en AR F-_. (+0); y=0o; a ERTRAGEN (2125b,c,'@); 
df 9. BEE Ta u Te 
Be: ae =V; a ae: A ß ( 4a4,D,C, ) 


Dabei sollen die beiden letzten Randbedingungen (22e, f) nur für die II. Näherung gelten. 
Im einzelnen ist zu den Randbedingungen folgendes zu bemerken: Die Bedingungen (21a, c) 
und (22a, c) stimmen mit (13) überein. Die Bedingungen (21b, d) erhält man aus (10) bzw. m 
für z=0. Die Bedingungen (22b, d,e,f) dienen zur Erzielung eines sanften Überganges a 
Grenzschicht in die Außenströmung. Beim Polynom g ist für beide Näherungen die ee 
der freien Konstanten gleich der Anzahl der Randbedingungen, so daß nach Erfüllung der Bi : 
bedingungen die Geschwindigkeitsverteilung in Umfangsrichtung bis auf die ea 5 
bereits festgelegt ist. Dies kommt darauf hinaus, daß wir für die Geschwindigkeitspro! er 
Umfangsrichtung Affinität für alle Werte von w/a annehmen. Beim Polynom j dagegen \ = 
Anzahl der freien Konstanten um eins größer als die Anzahl der Randbedingungen. Nac # 
füllung der Randbedingungen bleibt somit in der radialen u ae jun Ä 
' eine Konstante frei [c, in Gl. (23) und (24)], die zusammen mit der Grenzschichtdicke ö a 
den beiden Impulsgleichungen (18a) und (19a) zu ermitteln ist. 


| ä ir natürli u i n d (11) finden können. 
ı) Di ß d. (19) hätten wir natürlich auch durch Integration der Gl. (10) un 
h Wir yählen für cr zadinle Be azimutale Geschwindigkeitsprofil die gleiche Grenzschichtdicke, 
was ansich nicht unbedingt nötig wäre. 
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Unter Beachtung der Randbedingungen (21) und (22) ergeben sich somit die Polynom- 
Ansätze für die Geschwindigkeitsverteilung in folgender Form: 


1. Näherung: f= at(4 31) + cl 3 +28) + 0 P(l—2t4 9) 


Re (23), 
g= 5 @-3:4®) 23) 
II. Näherung: f= at(10 — 15t+ 68) + ct(1 — 61? 81° — 31%) 
| + 0t2(1 —3: +32 —P) + (24). 
g=oa(l —2t+2 — 1%) 
Der Koeffizient c, ist für beide Näherungen durch Gl. (21b) festgelegt zu: 
1 2 ” 62 w2 + a? i = 
a5 (44 REIN, p) a RT Fe Tr nee ae Fe ( 2 
während f 
d ) j 
— FU ea In 2 Da) 
(a r [91 (2 ) 
die oben erwähnte noch freie Konstante darstellt. Ferner setzen wir zur Abkürzung: 
3 E | 
| )jadh=— .Näl 
| ke a Be ee (25b) 
=—20, II. Näherung 


Durch Einsetzen von (23) in die beiden Impulsgleichungen (18a) und (19a) erhält man 


für die I. Näherung?): h 
‘al: ’ 
radial: j 
v 33 19 1 5 1 17 33 4 
— — 2 EI An en _— —— — —— Re 2 re > 
; a AT RT AR ig ® SL i 
se 05.91.19 27 : 
Er 1 £ ee a en 1 Re re A ä E; 
azimutal: = mp Image (27) . 
Für die II. Näherung erhält man mit (24)®): eı 
radial: | 
v 127 52 1 3 9 23 23 ; 
Brut tat ze B), | 
azimutal: Rare 28 54 | % | | 
. 6 u 47 82 un 94 [23 — 47 @ 3 ee EN A RS ES (29). 


Wird die freie Konstante c, aus der Impulsgleichung in Umfangsrichtung, Gl. (27) bzw. (29), 
in die Impulsgleichung für die radiale Richtung Gl]. (26) bzw. (28) eingesetzt, so ergibt sich 
unter Beachtung der Beziehung (25) für c, eine Gleichung für die Grenzschichtdicke 6, die sich 
für beide Näherungen in folgender Form schreiben läßt: 
62 62\8 62\2 62 
2 Be 2 al ce u ze 
er ker ek (103) + (004 Da) are 


j ®) Um die Gl. (26), (27), (28) und (29) zu erhalten, sind entsprechend den Impulssät 18 | 
die folgenden Integrale unter Benutzung der Ansätze der Geschwindigkeitsverteilungen (as) Rs Be 


zuwerten: \ 
Integral | I. Näherung | II. Näherung 
dt Ba 1 er 
Hl Berater zermatt 
2 19 1 181 52 1 | 
— 4? + 08 a he 1 Er ee N f 
fpdı 7 : AR I n ® arzt | 
- 37 A 17 hrs “81 23 
& tt Hast jagt 
1 
2d 33 2 23 | 
{9 i0® i96 
1 9 17. 19 3 
R a (7, ( 47 13 
HELL twitter)” na tmnıtinr)® 


ER ie, u, 
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Hierbei bedeuten A, B, ©, D, E konstante Zahlen, die für die beiden Näherungen in Tabelle1 


angegeben sind. 


Tabelle 1. Koeffizienten der G1. (30) zur Bestimmung der 
Grenzschichtdicke ö 


en 


| I. Näherung | II. Näherung 
A —0,093 89 —0,011 47 
B 0,313 61 0,082 14 
C 1,697 98 0,604 35 
D 8,155 40 3,294 61 
E 1,975 73 1,204 04 


In GI. (30) haben wir eine Gleichung 4. Grades für das Quadrat der Grenzschichtdicke 
ö°/10v. Die physikalische sinnvolle Wurzel ist hieraus in Abhängigkeit von ® und a zu 
ermitteln. 


V. Ergebnisse 


a) Grenzschichtdicke 
Sind und a vorgegeben, so läßt sich Gl. (30) entweder graphisch oder numerisch nach ö 
- auflösen. Wir haben die numerische Lösung vorgezogen. Zwecks Einführung von dimensions- 
losen Größen unterscheiden wir die folgenden beiden Fälle: 


Hal A: 04 
Wir setzen 
rd 62 
p == Pr und 0 0» N EN A ae Di Kein (31). 
Damit ergibt sich aus Gl. (30): 
A1+ Pa + Bol+P)® (CO + Dp)®+Epa=1..... (32). 


Die Auflösung von Gl. (32) für p=0, 1/4, 1/2, 3/4 und 1 ergibt die Werte nach Tabelle 2. 
Sie sind in Bild 2 (obere Kurve) dargestellt. Für den Fall der nicht angeblasenen rotierenden 
Scheibe ist « = 0 und damit p= 0. Der Wert 
der I. Näherung ist mit dem von W. G.Coch- 


o 1.Näherung 


ran [5] korrigierten Kärmänschen Wert 
identisch. 
Fall B: o<Sa 
Wir setzen 
a 20? (33) 
er a 
Damit ergibt sich aus Gl. (30): 
AUT HBÜHFP)P 
+ +DP+Eß=1 (3). 


Die Auflösung wurde für q=0, 1/4, 1/2, 3/4 
und 1 vorgenommen; siehe Tabelle 2 und 
Bild 2 (untere Kurve). Für den Fall der ange- 
blasenen, aber nicht rotierenden Scheibe (ro- 
tationssymmetrische Staupunktströmung) ist —=0 und damit de 0. E 

Das von S. Tomotika [11] angegebene Näherungsverfahren für die Grenzschicht an einem 
Rotationskörper liefert für die rotationssymmetrische Staupunktströmung 


51 29 n *) —2, wobei A=10ß 


Bild 2. Grenzschichtdicke öin Abhängigkeit‘ von a/» bzw. o&/a 


en ce en Ne 


A (05 2520 A 3024 


und hieraus A = 4,716 und somit IE no 2,172. 


ieht aus Bild 2, daß sich die Grenzschichtdicke nur wenig ändert, wenn man die 
Be engsgacheindigkäit, d.h. a konstant hält, und die Winkelgeschwindigkeit Se 
(untere Kurve). Hält ‚man dagegen die Winkelgeschwindigkeit ® konstant, dann nimmt die 
Grenzschichtdicke mit wachsender Anströmungsgeschwindigkeit ab (obere Kurve). De 
hältnismäßig großen Unterschied der Werte von ö nach der I. und II. Näherung kommt keine 


ae a Wr, ei; Ä ik 

ar 

. GER er Tr 
Be” BE N 
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Tabelle 2. 


richtung 7y, sowie axiale Zuströmgesc 
mungsgeschwindigkeit zu Winkelgesc 


Grenzschichtdicke ö, Wandschubspannung 
hwindigkeit cin 


in radialer Richtung z, und Umfangs- 
Abhängigkeit vom Verhältnis Anströ- 
hwindigkeit a/o (vgl. Bild 2, 3, 4) 


7 c 

a @ u SER Yre be 
Ba a & je ö V kin 008 Teo 008 | Yo» 
0%) © 1,307 3,615 0,534 0 ‚553 —0,606 
(0,780) (2,793) (0,544). (—0,537) (—0,55) 
< 1/4 4 0,937 3,061 0,597 —0,65 0,152 
= 1/2 2 0,648 2,546 0,806 —0,7185 0,571 
Eu 3/4 4/3 0,479 2,189 1,120 0,914 0,840 
/ 1 1 0,373 1,931 1,531 —1,036 1,028 

(0,233) (1,526) - (1,500) (—0,981) 

C 

a T Tp 
| ß Pr ö V rea Pi, | Yrea v8 Ya» 
1 1 0,373 1,931 1,531 —1,036 1,028 
m 4/3 3/4 0,383 1,957 1,412 —0,766 1,069 
= 2 1/2 0,390 1,974 1,312 —0,506 1,098 
& 4 1/4 0,395 1,987 1,277 0,252 1,119 
002) 0 0,397 1,992 1,255 0 1,126 

(0,248) (1,575) (1,222) (0) 
or? ar? 
Tu = —— Te4 u 


v 


1) Rotierende Scheibe in ruhender Flüssigkeit. 2) Rotationssymmetrische Staupunktströmung. 
Br Die eingeklammerten Werte gelten für die I. Näherung. Ä 


Bedeutung zu, da die willkürlich festgesetzte Grenzschichtdicke ö als reine Rechengröße anzu- 
sehen ist. Erst aus den Werten der Wandschubspannung (Tab. 4) wird sich ergeben, ob die 
II. Näherung eine Verbesserung gegenüber der I. Näherung bringt. 


b) Wandschubspannung 


Wir wollen jetzt die Wandschubspannung in radialer und azimutaler Richtung und aus 
der letzteren sodann das Drehmoment ermitteln. Für r, und z, ergibt sich: 


er ou d r 

radial: ne) =ur(#) u 0 £ FD 
} ov r 

azimutal: mul), =ur(®) =u3 d, - 37 


Die Konstanten c, und d, sind durch Gl. (27) bzw. (29) und (25b) gegeben. Dabei wird in 
(27) bzw. (29) c, nach Gl. (25) eingesetzt. Wir erhalten dann 


le) ae)” ee re 


Dabei bedeuten P,Q, R, S konstante Zahlen, die für die beiden in Ziffer IV beschriebenen 
Näherungen in Tabelle 3 angegeben sind. . \ 


Tabelle 3. 2 
Koeffizienten der Schubspannung nach GI. (38) und (39) 
I. Näherung - | II. Näherung 
p 0,195 32 0,423 88 
Q —0,25112 —0,363 33 
R 0,294 53 0,218 67 
Ss 0,474 34 0,632 46 


Um die Schubspannung in bequemer Weise dimensionslos 


entweder die örtliche Umfan 


geschwindigkeit der Potentia 
Radius r gebildeten Reynolds-Zahlen seien 


= 


Iströmung U = 


ro 
y y 3 164 


uemer zu schreiben, benutzen wir 
gsgeschwindigkeit der Scheibe = rw, oder die radiale Außen- 


ar. Die mit diesen Geschwindigkeiten und dem 


a ddh): 
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Wir-können dann die dimensionslosen Schubspannungen schreiben in der Form: 
T EENU 
Vre, — oder Vre, —., wobei r=r(r). 
0% oU 


Wir wählen die erste Schreibweise für den in der vorigen Ziffer näher definierten Fall A 
(» > a), und die zweite für den Fall B (» <a). Dies hat den Vorteil, daß man dann in be- 
quemer Weise die beiden Grenzfälle a—> 0 (nicht angeblasene rotierende Scheibe) und — 0 
(nicht rotierende angeblasene Scheibe — rotationssymmetrische Staupunktströmung) erhalten 
kann. 

Im einzelnen gilt mit den Bezeichnungen nach Ziffer Va: 


Fall’ 4: w>a: 


re. = [P+Q:p-a+ Bl + parla Bd 

0 

re = Sail ST ER ee aA]: 
0 


Für den Fall der nicht angeblasenen rotierenden Scheibe (a= 0, p= 0) erhält man mit 
.den Zahlenwerten von Tabelle2 die in Tabelle 4, Fall A, angegebenen Werte. Der 


\ Tabelle 4. Vergleich der Wandschubspannung nach der Näherung und der exakten 
Rechnung für die beiden Grenzfälle 
Fall 4,: a=0, p=0; nicht angeblasene rotierende Scheibe, exakt nach [5]. 
Fall B,: ®=0, qg=0: angeblasene nicht rotierende Scheibe, rotationssymmetrische Staupunkt- 
strömung, exakt nach [9], verbessert nach [12]. 


I. Näherung | II. Näherung | exakt N er fin 
anf Tr - 
Too Gi 0,544 0,534 0,510 = 
Fall A, 
Vreo —0,537 0,553 0,616 je 
_v%, 2 x 
Vre; en 1,222 1,255 1,312 1,283 
Fall B, e 
gr, 
16a 003 0 0 0 0 


Wert der I. Näherung ist mit dem von W. G. Cochran korrigierten Kärmänschen Wert iden- 

tisch. Man sieht, daß unsere II. Näherung dem exakten Wert etwas näher kommt als die 
z I. Näherung. 

Für die II. Näherungsind die Schubspan- 

nungen auch für p = 1/4, 1/2, 3/4 und 1 berech- 

net worden. Sie sind in Tabelle 2 angegeben 

und in Bild3 (untere Kurve) und Bild 4 


a ) vn 
grompeik a 44T | IT 
Bao 
Far eren 
a = 


R - T 
== eo pr 
{7 

« Iren 6 = fl@) ?, 

Wräma, = We 

See BE EN En BE EN 
= > 2 0% Ze 0 
Sr => 
Bild 3. Wandschubspannung 7, in radialer Richtung Bild 4, Wandschubspannung z,, in Umfangsrichtung. 


= 


u (obere Kurve) dargestellt. Man sieht, daß bei festgehaltener Drehgeschwindigkeit beide Schub- 
spannungen mit wachsender Anströmungsgeschwindigkeit erheblich zunehmen. 


Be 


_— 


a 7 04 dr SR 
t Ten WERTEN ER iu 
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Fall 2: o>a: 
Vre,y —-[P+QB + RU +) PB? .... .... (43), 


Bm 


Für die Staupunktströmung (o=0, q=0) erhält man mit ‘den Zahlenwerten von Ta- 
belle 2 die in Tabelle 4, Fall B,, angegebenen Werte. Auch hier kommt die II. Näherung dem 
exakten Wert etwas näher als die I. Näherung. Der Näherungswert von Tomotika ergibt sich 
aus dem Pohlhausen-Verfahren für einen Rotationskörper zu 


Vre, = 2+2) 22 = 1,2835° mit 24462 aa (45). 


Für die II. Näherung sind die. Schubspannungen auch für g= 1/4, 1/2, 3/4 und 1 berechnet 
‚ worden. Sie sind in Tabelle 2 angegeben und in Bild 3 (obere Bl und Bild 4 (untere Kurve) 
dargestellt. 
c) Drehmoment 

Man kann die vorstehenden Ergebnisse, die sich auf eine unendlich ausgedehnte Scheibe 
beziehen, auf eine endliche Kreisscheibe übertragen, wenn deren Radius R groß ist im Ver- 
gleich zu der Dicke der mitgerissenen Schicht ö. Das Drehmoment einer solchen einseitig be- 
netzten Scheibe ist: 


R 
M=-2nr [r PR a IF IESER AE NE ER  e- (46). 
N) 


5 { '  (d r T 
Nach GI. (16) ist in Verbindung mit (8) 7,— ur (2) - T),_RR ‚ und damit 


ee ee et (47). 
Führt man in üblicher Weise den dimensionslosen Momentenbeiwert durch f 
NN (48) 
e 2 Rö x 
ein, so erhält man aus (47) und (48): 
Cu=F [| . _ N RE N (49) 


2 
Dabei bedeutet Re = = die auf den Scheibenradius und die Umfangsgeschwindigkeit 
bezogene Reynolds-Zahl, während 5 


1 a (To)r-R 
„Rl2)=—VRe- en 


die dimensionslose Wandschubspannung darstellt, die durch die vorletzte Spalte von Tabelle 2 


\ 


gegeben ist?). Die Werte von »(2) = (y VRe sind in Tabelle 5 angegeben und in Bild5a_ 


über a/w aufgetragen. Für die in ruhender Flüssigkeit rotierende Scheibe (@ > 0) ergibt die 


exakte Lösung Oy /Re—= 1,935, während nach unserer Näherungslösung dieser Wert gleich 
1,74 ist. 


Tabelle 5. Drehmomentenbeiwert einer rotierenden Kreisscheibe in Abhängigkeit 
vom Verhältnis Anblasgeschwindigkeit zu Umfangsgeschwindigkeit 


ao | 0 er | 3a | 1 | an | 2. 2% 
Cu Re | 1,74 | 2,05 | 2,46 | 2,47 | 3,25 | 3,69 | 4,49 | 6,35 


“) Die Werte der unteren Hälfte der vorletzten Spalte von Tab-2 sind umzurechnen nach der Beziehung: 


Vs mega Ss: 


—— 
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! Wie man aus Bild4 und Tabelle 2 ersieht, nimmt die Wandschubspannung Vre, - BB 
linear mit w/a zu, und zwar gilt 003 
Sr 0) (0) 
Vreu a —— 1,01 m . = < 1) . 
Damit läßt sich der Momentenbeiwert nach (49) auch in der Form schreiben: 
Er lREe 1.01% VE-srJ2 REN CH Age 
0) (0) 


Die Kurve nach dieser Formel haben wir in Bild 5a mit eingetragen. Vollständige Über- 
einstimmung mit den Werten nach Tabelle 5 herrscht für — > = —= 1,27. Aber auch bis zu 
») ' IT 
Werten = > 3 0,64 ist die Übereinstimmung noch recht gut. 


Um dieses Ergebnis für das Drehmoment in bequemer Weise auf die angeströmte Scheibe 
vom Radius R anwenden zu können, ist es schließlich noch erforderlich, die Dimensionslose a/o 
unmittelbar in Zusammenhang zu bringen mit dem Geschwindigkeitsverhältnis W/V,, wo W= 
die Anströmungsgeschwindigkeit und V,= Row die Umfangsgeschwindigkeit der Scheibe be- 
deutet. Zu diesem Zweck muß eine Zuordnung der rotationssymmetrischen reibungslosen Stau- 
. punktströmung gegen die unendlich ausgedehnte ebene Wand zu der reibungslosen Strömung 
gegen eine Kreisscheibe vom Radius R gefunden werden. Eine sinnvolle Zuordnung dieser beiden 
Strömungen erhalten wir aus der Überlegung, daß die beiden Geschwindigkeitsfelder in der Um- 
gebung der Staupunkte übereinstimmen sollen. Dies führt zu der plausiblen Forderung, daß 
in den Staupunkten beider Strömungen die Geschwindigkeitsgradienten senkrecht zur Wand 
- gleich sind. Die Auswertung dieser Forderung ergibt: 


BEE IS re a I 


/// 


3 Ö 
1 NIS 
een 
N Se rum —— roch Tabea 5 
6 en - 317% (Nöherung) 
Ban Ta 
Pag 
pe Laminar (exakt) eg 
F nach Cochran 5 
E 2 
R 
3 
R 
So = 
06 
N 


IS: 
> 


ö 2 
Reynolds-Zahl _ Re= Re _ 


beiwert der rotierenden Scheibe bei verschiedenem Verhältnis von Anströmungsgeschwindigkeit zu Umfangs- 
Ei d ge geschwindigkeit W&/Ro = Wol[Vo 
_ Hiermit ergibt sich dann aus (49a): 2 
= ; ee 2,33. /We Dr: 
. j MH Von 
= VYRe ‚ V, 
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We 
Entsprechend den oben gemachten Angaben gilt diese Beziehung exakt für 1 >» 


aber sie ist auch noch gut brauchbar für a >1. Gleichung (49b) besagt, daß bei fest- 


gehaltener Drehzahl das Drehmoment sich ändert wie die Wurzel aus dem Verhältnis der An- 


 blasgeschwindigkeit zur Umfangsgeschwindigkeit der Scheibe. 


In Bild 5 haben wir den Mo- 


mentenbeiwert über der Reynoldsschen Zahl aufgetragen und W./V, als Parameter gewählt. 
Mit eingezeichnet haben wir die von v. Kärmän [4] gefundene Kurve für die nicht angeblasene 
rotierende Scheibe bei turbulenter Strömung?). Die aufgetragenen Beiwerte gelten für die ein- 
seitig benetzte Scheibe. Wenn wir näherungsweise annehmen, daß hinter der Scheibe infolge 
des Totwassers keine Geschwindigkeitskomponente in Achsrichtung herrscht, dann können wir 
den Beiwert der beidseitig benetzten Scheibe wie folgt berechnen: 


wobei Cm = 


. ®) Für die einseitig benetzte Scheibe gilt: cy = = 


SEE 
NER || 


CM beidseitig = CM einseitig 7 CM; > 


” 


IT LE 
dEENEERBLZE 
= Erz 
2 


e 
IE 


> 


RER LEN SINGEN. 
BERENENDRE 


a 
SE 
re 
= 
Ba 


Ba 
= 
Bi 
5] 
Se 
Eu 
4 


el 


den Momentenbeiwert für We — (0 bedeutet. 


d) Geschwindigkeitsverteilung 

Die Geschwindigkeitsverteilung in 
radialer und azimutaler Richtung wird 
durch Gl. (8) zusammen mit Gl. (23) 
und (24) gegeben. Es seien hier nur 
die Ergebnisse für die zweite Nähe- 
rung diskutiert. Indem Ausdruck für f 


nach Gl. (24)sind die Koeffizienten c, - 


und c, nach Gl (29) bzw. (25) einzu- 
setzen. Führt man für die Polynome 


in t=- nach GI. (24) die folgenden . 


.Ö 
Abkürzungen ein: 


%=#(10 —15t+61) 
—= 1—(1—1)? (143 1+-6 2?) 
G=t(1 —61?+ 81° — 81) 
=t1—N(l +31) 
G=#(1 —3t+32 —?) 
= (1 —t)° 
S=1—2: +21 
u A IE ES, 
so läßt sich die dimensionslose Ge- 
schwindigkeitsverteilung für die in 
Abschnitt Va definierten Fälle A und 
B in folgender Form schreiben: 


(51), 


FallA: ® >; =, 
5 {0} 
RE 5 c | 
IE bergen Fe ae N, @3 
Bild 6, Radiale und azimutale Geschwindigkeitsverteil FallA, G ä EI 
ve 
Kurve: Exakte Lösung für die in Tenaıfe Flüssigkeit” Peila, aa, Gesicht ® Rd g 
Cochran [5] ——=——— = G; RR (53). 
or [73 (0) 
Fall B: o<a; g= 
a 
WET N Ca 
PIE we oe ot G-+ @, ER TE TER ee (54), 
® ® -g 
a u N (55). 
0,073 


\ 
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Die Zahlenwerte der Konstanten c/o, 6,/o, cı/a, ©,/a sind in Tabelle 6 angegeben. Sie 
lassen sich aus Gl. (29) und (25) mit den Werten von ö nach Tab. 2 leicht berechnen. 


Tabelle 6. Die Konstanten c, und c, zur Berechnungder G eschwindigkeitsverteilung 


nach Gl. (51) bis (55) 


| p= ao q=ovla c/o 2910) ER. c,/a 
= ) © 1,929 —6,535 Em 2 
= 1/4 4 1,832 —4,978 ax Er 
= 1/2 2 2,054 —4,050 En = 
Mr 3/4 4/3 2,454 —3,742 = er 
| 1 er 2,961 —3,730 2,961 23730 
= 4/3 3/4 en = 2,764 —2,992 
= 2 1/2 = _ 2,625 — 2,438 
B 4 1/4 2 = 2,535 —2,098 
© f) £ == 2,502 —1,985 


Mit den in Tab. 6 angege- 
benen Werten von c, und c, sind 
die Geschwindigkeitsverteilun- 

. gen für die Parameter p—= 0, 1/4, 
1/2, 3/4, Lund g=0, 1/4, 1/2, 
3/4, 1 berechnet worden und in 


Bild6für den FallA über dem ; 


dimensionslosen Wandabstand 
z/o/v und in Bild 7 für den 
Fall Büber z  a/v aufgetragen®). 
Die beiden Grenzfälle der nicht 
angeblasenenrotierenden Scheibe 
und der rotationssymmetrischen 
Staupunktströmung erhält man 
aus FallA für p=0 bzw. aus 
Pal Peg —U. 
Fälle sind auch die Geschwindig- 
keitsverteilungen der exakten 
Lösungen nach W. G. Cochran 
[5], Bild6, und F.Homann [9], 
Bild 7, miteingetragen. Bei der 
Staupunktströmung ist dieÜber- 
einstimmung mit unserer Nä- 
ausgezeichnet. 
Bei der nicht angeblasenen To- 
tierenden Scheibe weicht am 
‚äußeren Rande der Grenzschicht 
die Näherungslösung um einiges 
von der exakten Lösung ab; 
doch kann auch hier die Über- 
einstimmung beider noch als 
befriedigend angesehen werden. 


herungslösung 


Eine 


anschauliche 


Für diese 


Vor- 


stellung von der Geschwindig- 
 keitsverteilung gibt auch das 
Polardiagramm in Bild 8. Hier 
ist von einem Pol aus der Ge- 
‚schwindigkeitsvektor mit den Komponenten 4 und v für verschiedene Wandabstände z dar- 


gestellt. 


Axiale Zuströmungsgeschwindigkeit:I 
ee neue nach Gl. (6) 


= 


% 


a 
BE 1 | [azimucat | 
Ra SEE 
EINE BERNER 
BEER FA ER 
a ee 

NSS er 
BR SRSEANTESEEN 
Ba ren 
SENSE 


5 YFz 


7 A EEE 
FEEE 


Re 
ee 
DA 
= 
S 
N\ 
35 


Bild 7. Radiale und azimutale Geschwindigkeitsverteilung, FallB, a>w. Gestrichelte 
Kurve: Exakte Lösung für die rotationssymmetrische Staupunktströmung, nach 


Homann [9 


n der reibungslosen Strömung ist die axiale 


a Zur 


6) Die universellen Funktionen 9 bis Qi für die Geschwindigkeitvertehun sind in Bild 9 graphisch 


£ ? men 


R “  Z.angew. Math. Mech. 
110 Schlichtin gu. Truckenbrodt, Die Strömung an einer Scheibe Bd. 32 Nr. 4/5 April/Mai 1952 


——radiale Richtung lee ne : 
05 u/Vo H l 15 
Br 15 
E — 
ee I 
S W < = 70 
R © 05 
.O 05 & 
$ S 
Q 
S 
S 


05: uU 4 
Br z 
er v 
E 
2 05 2-05 


Six 
5 ; 


Bild 8. Polardiagramm der wandparallelen Geschwindigkeitskomponente 
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Durch die Wirkung der Reibung wird die axiale Zuströmungsgeschwindigkeit geändert, 
und zwar‘bei der nicht drehenden Scheibe (Staupunktströmung) abgemindert, weil die radial 


ne 16 Be 
5 : ITd. el er rehenden ’ n 

IN 56 ja 8 

at 

ENETANZER 

A I 

BEN WERE 
x 


‚aus der Kontinuitätsgleichung (9): 


weh) —2 [fd 3.66): 


Scheibe bewirken die Zentrifugalkräfte eine 
Vergrößerung der radial abströmenden Menge, 
und infolgedessen nach der Kontinuität auch 98 
eineVergrößerungder axialen Zuströmgeschwin- 
- digkeit. Für den allgemeinen Fall ergibt sich 

4-4 N 
2DSNZENERE 
EA ATER DZ ISNER ER WS EL 
READER EIER 
VEABEESSNE 
- /} R (2 (2 / ’ 


Sehr weit vor der Scheibe (—oo) gilt: 
vwo)=h(R)=—2azr+c .. (57). 


Dabei stellt das erste Glied den Anteil der‘ 
reibungslosen Strömung dar, und das zweite 
‚die Anderung infolge der Reibung. Wir wollen 


\ ' nicht näher auf die Berechnung von c ein- De 10 
gehen, sondern verweisen lediglich auf die in _ t= z 
Tabelle 2 (letzte Spalte) mitgeteilten Zahlen- & 
s Kal n Bild 9. Die universellen Funktionen @, bis G, für die' 


ergebnisse. > ip Geschwindigkeitsverteilung nach (51) 


x VI. Zusammenfassung | 
_ Es wird die lJaminare Strömung in der Nähe einer unendlich ausgedehnten rotierenden 
Scheibe berechnet, die gleichzeitig in Richtung der Drehachse angeströmt wird. Die beiden 
 Grenzfälle der nicht angeblasenen rotierenden Scheibe und der nicht rotierenden, aber ange- 
blasenen Scheibe (rotationssymmetrische Staupunktströmung) sind bekannt. Die hier mit- 
geteilte Lösung des allgemeinen Falles wird nach einem Näherungsverfahren mit Hilfe des 
Impulssatzes gewonnen. An der Scheibe bildet sich eine Grenzschichtströmung mit radialen 
und azimutalen. Geschwindigkeitskomponenten aus. Die Dicke der durch die Drehbewegung 
 mitgerissenen Schicht (Grenzschicht), die Geschwindigkeitsverteilung in dieser Schicht und das 
Drehmoment der Scheibe sind stark abhängig von dem Verhältnis der Anströmungsgeschwindig- 
keit zur Umfangsgeschwindigkeit. Bei fester Drehzahl nimmt mit wachsender Anströmungs- 


geschwindigkeit das Drehmoment stark zu. 
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Der Sonderschieber für Häufigkeitsrechnung 


Von W.de Beauclairin Stuttgart 


An Hand von kennzeichnenden Beispielen wird gezeigt, daß der neuentwickelte Sonderschieber für 
Häufigkeitsrechnungen ein vielseitig und vorteilhaft brauchbares Hilfsmittel ist, das auch dem mit der Theorie 
der mathematischen Statistik nicht Vertrauten ohne große Mühe verständlich und dienlich ist, um jeden 
Versuchsbericht oder Erfahrungsniederschlag statistisch exakt zu formulieren und zu beurteilen. 


The author describes a newly-constructed special slide rule for stalistical compulations. He demon- 
strates, by several characteristical examples, that the slide rule is a very useful auziliary instrument. It 
enables even persons that are not well acquainted with the theory of mathematical statistics to formulate each 
experimental report in a statistically exact manner and to evaluale them easily. 


Un curseur sp6cial r&cemment construil, par sa maniabilite et son universalit&, facilite V’&laboration de 
beaucoup de problömes du calcul pratique de frequence et du contröle de conclusions stalistiques. Ses Echelles 
logarithmiques sont tonnues et familißres ü l’ingenieur par l’usage du curseur glissant de caleul. Le curseur 
special donne aussi au praclicien, qui ne peut pas s’enfoncer dans la stalistique malh&matique, la possi- 
bilite d’achever ses exp£riences, ses invesligalions el ses conclusions par des r£sultals slatistiques incon- 
testables. S 


BhHoss pasBuran cHemmaıbHaa pacuerHan mHelika oÖneryaer CBoea YAODHOCTBO HU YHH- 
BEPCANBHOCTBE O6PaÖ60TKy MHOTHX 3a1a4Y IIPAKTHYeCKOTO BEIYHCHEHHA YACTOT, A TAKıKe KOH- 
TPONb CTATUCTHYeCKHX 3aK.moyeHul. Bere yIoTpeöseHune ee 10rapuhMHyeckKuX IIKası Odecne- 
yeHo HH:KeHePy AHAJIOTHYHOCTBIO C IIPHBEIYHOH orapuhMmHugeckol AHHeÜKON; TAKHUM O0Pa30M, 
H IPAKTHKy, AeTAIBHO He OCBOHBIIEMYCA C MATEeMATHYeCKON CTATHCTHKOH, 13H BO3MOKHOCTB, 
34BePINHTb CBOH OIIEITEI, HCCHeNOBAHHA H 3AKIIOYeHHA CTATHCTHYeCKH Öe3yIpe4HBIMH PesyJıb- 


TATaMH. 
Einleitung 


Ein neuentwickelter Sonderschieber erleichtert durch seine Handlichkeit und seine Uni- 
versalität die Bearbeitung vieler Aufgaben der praktischen Häufigkeitsrechnung und der Kon- 
trolle statistischer Schlüsse. Seine logarithmischen Leitern sind dem Ingenieur vom Rechen- 
schieber her gewohnt und geläufig, er gibt daher auch dem Praktiker, der sich nicht eingehend 
in die mathematische Statistik einarbeiten kann, die Möglichkeit, seine Erfahrungen, Unter- 
suchungen und Folgerungen mit statistisch einwandfreien Ergebnissen abzuschließen. 


Bezeichnungen: 
Gaußsche Normalverteilung von Fehlerhäufigkeiten: h(z) = h, exp — 0,5 (x — a/o)%. 
Argument, Merkmal, Meßwert, Fehler: x. 
Streuungsvielfache als normiertes Argument: s— 2/0 (a—=0). 
Häufigkeitszahlen: A. 
Scheitelwert der Häufigkeit: A.. 
Scheitelwertabszisse: a (bei Gaußverteilung: = M). 
Prozenthäufigkeit: h/h,—= h%. \ 
Relative Häufigkeit: h/n (für große n: — Wahrscheinlichkeit: ?). = 
Klassennummer der Merkmalklassen: k. 
Summenprozente des Fehlerintegrals: H %. 

: Gesamtzahl der Proben, Messungen, Umfang des Kollektivs: n. 

Arithmetischer Mittelwert der Verteilung: M. 
Mittlere quadratische Streuung: o. x 
Mittlerer Fehler des Mittelwertes: Iu= a 


Herleitung des Sonderschiebers aus dem Gaußnetz 


Eine der wichtigsten Aufgaben der Häufigkeitsrechnung, insbesondere der statistischen 
Fertigungsüberwachung, besteht darin, die Häufigkeitsverteilung einer Reihe von unbeabsichtigt 
und unsystematisch streuenden Merkmalen mit der als ‚Lehre‘ dienenden Normalform der 
Gaußschen Fehlerfunktion | 

h= hu exp — 0,5 (x — a/o)? 

zu vergleichen. Obwohl diese Funktion selbst tabelliert ist und auch die Grundfunktionen y—=er 
und y=e”*in Zahlentafeln vorliegen und auf besseren Rechenschiebern als loglog-Leiter ent- 
‚halten sind, ist die numerische Rechnung für die Praxis doch zu zeitraubend und umständlich. 
Es gibt daher als Hilfsmittel zur bequemeren Auswertung, wenn auch in geringerer Genauigkeit, 
die im Handel erhältlichen sog. „‚Häufigkeits-“ und „Wahrscheinlichkeitspapiere“, die durch eine 
stattliche Folge von Veröffentlichungen der Gruppe um Daeves bekannt wurden und für 
manche Zwecke auch von bestem Nutzen sind, 

Das ‚‚Wahrscheinlichkeitsnetz‘, kennzeichnender „Summenprozentnetz‘ zu benennen, hat 
eine solche Ordinatenteilung, daß das Fehlerintegral oder die „Summenprozentkurve‘ in eine 
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Gerade verstreckt wird. Es ist daher bei der Anwendung dieses Netzes notwendig, die Häufig- 
keitswerte h; zu summieren, die Prozentzahlen h/n zur Gesamtsumme £h=n zu bilden und 
diese ‚Werte in das Netz einzutragen. Durch die sich ergebende Punktfolge ist dann nach Augen- 
maß die bestens ausgleichende Gerade zu ziehen. An dieser können die gesuchten charakteri- 
stischen Größen der die Meßwerte annähernd ersetzenden Gaußschen Glockenkurve abgegriffen 
werden: die Scheitelwertordinate A,, die Scheitelwertabszisse a und die Streuung o. 

Das „‚Häufigkeitsnetz“ ist ein halbes Wahrscheinlichkeitsnetz, in dem auch die Gaußsche 


Fehlerfunktion sich als Dreieck aus Geradenstücken abbildet — allerdings nur sehr annähernd 
nur in den äußeren Teilen. 


Ein anderes Funktionsnetz at % 5% 
ist das von Lorenz!) einge- 2998 
führte „Prozenthäufigkeits- 4, | 
netz“,dasvonihm,‚Gaußnetz“ | a 
genannt wird. Es ist leider noch „ 


nicht im Handel. Darin liegen 


die HäufigkeitswertehY\, —=h,/h, 5 

genau aufeiner Geraden (Bild1), „ N a 

falls sie vorher in Prozentzahlen 5 

des Scheitelwertes A,, d.h. der » n # 
größten Häufigkeit, umgerech- » a we s 

net wurden und falls sie natür- * 

lich der Gaußschen Fehlerfunk- ” £ 

tion genügen. Auch hier werden „41 r } N EEE Kl 

die mehr oder weniger streuen- s | S x 

den Punkte durch eine Aus- » ? R 2 

gleichsgerade vermittelt und an 5 3 60 

ihr die charakteristischen Werte e ’ i 

der Ersatzfunktion abgelesen. ERII TIEEERL hie Y 

Hierzu ist also ebenfalls ein „ 2 

Netzpapier notwendig, das man „ % 


| 

| 

| 
sich selbst anfertigen muß, wei- x | 
ter das Umrechnen in Prozent- 7 
zahlen, das Einzeichnen, Aus- i 


gleichen und schließlich das „ 
Rückrechnen der durch die Ge- „ KeieT. Se SER IF er 
rade verbesserten Prozentwerte > 
in Häufigkeiten. » g 
Allen diesen Verfahren ist ,; Lyr ; 
gemeinsam, daß das Netzpapier NR 
vorhanden sein muß und ein . | - 02 
Rechenschieber außerdem, und 7’- BR 
. daß eben nur diese behandelte « : 005 
Aufgabe damit zu bearbeitenist. Br 
Ein allgemein brauchbares ® 
Rechenhilfsmittelhat demgegen- » = 
über die größten Vorteile, wenn. Bild 1, Kombiniertes Prozenthäufigkeits- und ‚Summenprozent-Netz, 


es auch die sonstigen - oft vor- Lorenzsches Gaußnetz Teilung links, Wahrscheinlichkeitsnetz Teilung rechts 
kommenden Aufgaben der prak- Eingetragenes Beispiel: Höhe von 125 Kiefern. M = 176,2 cm. o = 36,7 cm 
tischen Häufigkeitskontrolle er- 
leichtert und bequem zu bearbeiten erlaubt, wenn es also alle Funktionswerte mit genügender 
- Genauigkeit darbietet, wenn weiter die wichtigsten Abhängigkeiten vor Augen stehen und sich 
sämtliche speziellen Netzpapiere, Zahlentafeln, Rechen- und Zeichenhilfsmittel erübrigen. 
Das Prozenthäufigkeitsnetz bietet die Grundlage zu einer solchen Entwicklung. In diesem 
Netz liegen doch die Prozentwerte jeder Gaußschen Fehlerfunktion auf einer Geraden, und die 
Variation der charakteristischen Zahlen drückt sich lediglich in unterschiedlicher Lage und 
Steigung dieser Geraden aus. Diese Geraden müssen nun durch eine logarithmische Funktions- 
leiter ebenso darstellbar sein wie etwa die einfache Gerade y=az durch die normale logarith- 
mische Rechenschieberleiter. Es ist also nur die Lorenzsche Ordinatenteilung mittels eines 
logarithmischen Maßstabes aufzutragen und auf eine normale logarithmische Argumentteilung 


1) P, Lorenz: Analyse von Verteilungskurven, Technik 2 (1947), 8.8388. 
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als verschiebliche Zunge zu beziehen: damit Be sich alle Lorenz-Geraden durch einfache 
iti i ider Leitern reproduzieren. 

ee ne: des Wahrscheinlichkeitsnetzes logarithmisch um- 
zeichnen und auf die logarithmische Grundleiter beziehen, und weiterhin auch die Leiter der 
Be .. ein „„Sonderschieber für Häufigkeitsrechnungen“ (Bild 2), dessen weitrei- 
chende und vorteilhafte Anwendung im folgenden beschrieben werden soll. 

Der Sonderschieber trägt folgende Teilungen: 

1. Vielfache der Streuung s— x/o im Bereich von s= 0,1 bis 4; i 

9%, Gaußsche Normalverteilung A%Y%, von Häufigkeiten in Prozenten des Scheitelwertes Ro; 
beziffert über A%, = 0,05% ** * 99,4% ; 

3. Fehlerintegral HY, in Prozentwerten, gegenläufig beziffert über H%=4:-49,99% und 
46 -:-0,01%. = 


h% 02,937 s 20 3o 50 a 700 9 50 30 70 $ Fa 


Entwurf: 


45 
Dr-Ing.W.de Beauclair 5 49 
ra 2% IS U 005 
02 1 5 20 2» 2 2 0 0 m 90 95 9,5 007 
H% ie hl 294 


Bild 2. Sonderschieber für Häufigkeitsrechnung 


4. Diese drei Körperleitern sind bezogen auf eine Zunge, die zwei gleichlautende logarith- 
mische Leitern der Merkmalsgruppen- oder Klassennummern X trägt, beziffert über k=0,5:--20. 

5. Oberhalb dieser Teilungen ist eine zweite Zunge mit normaler logarithmischer Teilung 
über zwei Dekaden angeordnet, auf die sich noch folgende Hilfsleitern beziehen: 

6. Eine korrespondierende logarithmische Leiter zum Umrechnen in Prozentwerte und 
andere übliche Rechenschieberaufgaben. 

7. Vielfache t des mittleren Fehlers des Mittelwertes 0, abhängig von der Zahl der Frei- 
heitsgrade m = n — 1, beziffert über t = 3,08 --- 20, d.h. für m = 2----100. j 

8. Eine Teilung 9, zur Beurteilung eines Nullergebnisses: der m-fache Wert der oberen 
(3 o-)Grenze der Grundwahrscheinlichkeit, bezogen auf die Anzahl der Beobachtungen n, be- 
ziffert über 9 = 2: - 5,75, d.h. für n= 2,2: - : 100. 

Ferner trägt der Rechenschieberkörper am oberen und unteren Rande 4 


9. eine als Maßstab dienende Lorenz-Teilung A%, zum Zeichnen der Geraden der Häufig- 
keitsprozente, sowie‘ y a BEE 


10. eine entsprechende Maßstabsteilung 4% zur Verstreckung der Summenprozentkurve. 


Das Ablesen der zusammengehörigen Zahlenwerte wird erleichtert durch einen breiten 
Läufer, der auf Vorder- und Rückseite je den üblichen Ablesestrich trägt. 


Dieser schräg über die Zahlenleitern laufende Läufer bedeutet weiterhin eine ganz wesent- 


liche Neuerung an .Sonderrechenschiebern. Durch diese Anordnung wird je nach der Stellung 
des Läufers über der Körperteilung eine andere Stelle des Läuferstrichs von dieser Leiter ge- 
schnitten, was ermöglicht, in diesen Schnittpunkten senkrecht zum Läuferstrich zusätzliche 
Abhängigkeiten aufzutragen. Diese Funktionskurven (eine oder mehrere) können den üblichen 
Läuferstrich ersetzen oder ihn ergänzen, derart, daß der Einstellung des Läufers auf der Körper- 
leiter ein oder mehrere abhängige Funktionswerte auf gleicher (oder anderer) Körperleiter zuge- 
ordnet werden. Damit läßt sich auch bei Rechenschiebern die sonst nur von Nomogrammen oder 


Big)" erreichbare Vielseitigkeit der Rechenmöglichkeiten einführen (Bild 3, Deckblatt zu 
ild 2). ee 


anzahlen n = 30, 50, 100, 300, 1000 und 5000 an der Körperleiter A%, abhängig von der Ein- 


£ er Sonderschieber trägt also zu beiden Seiten der beiden Läuferätriche auf Vorder- und 
Rückseite noch je sechs Kurven zum Ablesen von Wahrscheinlichkeitsgrenzen für Proben-. 
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stellung %%, mit dem Läuferstrich, wobei jedoch jetzt die Bezifferung h%, als relative Häufig- 
keit h/n oder bei großem n als Wahrscheinlichkeit p zu lesen ist. 

11.-Die kurze Leiter n,/n, arbeitet mit der Marke auf der linken Seite der oberen Zunge 
zusammen und dient zur Umrechnung der mittels der Läuferkurven gefundenen Wahrschein- 
lichkeitsgrenzen für den Fall, daß die Kollektive ungleichen Umfang haben.-Die Marke muß 


5000 
5000 5000 1000. __.300 n=5000 (hl 
700, - 700 1000 300 
300 / En. y 
100 — 30 30 
30 100 
30 50 
5 30 => 
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Inklusionsschluß Transponierungs- und Repräsentationsschluß 
Kurven zum Ablesen von 30-Fehlerbereichsgrenzen auf der Leiter h% 


Bild 3, Deckblatt zu Bild 2. Schrägläufer des Sonderschiebers für Hänfigkeitsrechnung 


dazu auf den Zahlenwert n,/n, des Verhältnisses beider Anzahlen eingestellt werden; dann gibt 
. die Körperteilung den nach n,/n, berichtigten Wert der auf der Zungenteilung beachteten Wahr- 
5 .  scheinlichkeit für gleiche Kollektivmengen (nach Koller)?). 
12. Einige Teilstriche n, an gleicher Stelle dienen zum Umrechnen eines Nullergebnisses bei 
gleichgroßen Kollektiven mittels der oberen Zunge. 
Im folgenden soll von mathematischen Darlegungen abgesehen und besser an Hand von 
Beispielen eine Reihe von Anwendungsmöglichkeiten und deren Vorteile erläutert werden. 


Anwendung des Sonderschiebers für Häufigkeitsrechnungen 


1. Das Aufzeichnen von Normalverteilungen 

5 Das Aufzeichnen von Gaußschen Normalverteilungen von Fehlerhäufigkeiten, die durch 
ikre charakteristischen Zahlenwerte 

 arithmetisches Mittel M, bei Normalverteilung M=a, 

Scheitelwert A, zuxz=a, a) 

Streuung 0 : > 
gegeben sind, wird durch die auf dem Sonderschieber zu beliebigem Argument abgreifbaren 

- Funktionswerte wesentlich erleichtert. Man ordnet zuerst, wenn es nicht schon vorliegt, der 

Merkmalteilung z eine Klassenteilung % bei, die zum Mittelwert « symmetrisch mit k=0 beginnt. 


Wie eng diese Klassen gewählt werden, ist allein die Frage’ der gewünschten Zeichengenauigkeit. 


‚Es empfiehlt sich, ist aber keineswegs Bedingung, einen ganzzahligen Teil der Streuung o als 

+ Klassenbreite zu wählen, was das Aufzeichnen erleichtert. Dann stellt man die untere Zunge 

3 mit dem dem Streuungsmaß o entsprechenden k-Wert auf die Marke bei s=1, d.h. Ay, =60,6%, 

F _ ein und kann weiter zu allen beliebigen Klassenmitten die zugehörigen h‘%, ablesen. Diese trägt 

"man entweder in einem passenden Maßstab gleich auf oder man rechnet sie mittels der oberen 

_ Zunge bei konstant bleibender Einstellung in Häufigkeitswerte x—=h%-h, um. Zu k=0, d.h. 
h%,=100% gehört der in der Schieberteilung ja nicht verwirklichbare Scheitelwert ho. 

Mittels der unteren Körperleiter 4%, lassen sich bei gleichbleibender Zungenstellung gleich 

_ die Zahlenwerte der Summenhäufigkeiten ablesen und danach zum Häufigkeitsbild auch die 

"Summenlinie zeichnen, sowie die anderen etwa interessierenden Zahlenangaben bestimmen: 

etwa die Spannen 7 als Doppeltes der „innen“-Werte, oder die „Gut“- und ‚‚Ausschuß“-Anteile 

als innen- und außen-Werte der H%-Teilung. RL 

ER 2. Das Prüfen gegebener Häufigkeitsverteilung auf ihre Normalform 
Zuerst zeichne man die vorliegenden Häufigkeitswerte h zu ihren Merkmalswerten x auf, 


I ri 


Er er 2 } Ei b 8% 


N > 
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rechnungen. Selbstverständlich darf durch diese Bearbeitung der Häufigkeitszahlen keine Ver- 
fälschung unterlaufen! : 2 
Alsdann verbindet man die Punkte durch eine bestens ausmittelnde Kurve. Ist diese Aus- 
‚ gleichskurve dem Augenschein nach einigermaßen symmetrisch zum Mittelwert undin ihrer Form 
der Glockenkurve ähnlich, dann ist eine weitere Auswertung mit dem Sonderschieber ohne 
weitere Rechnung aussichtsreich. j 

Man teilt nun die Argumentachse nach beiden Seiten von «in Merkmalklassen kein und 
beziffert diese mit ganzen Zahlen k=0 von z=a aus. 

Man lese dann den von der Ausgleichskurve bestimmten Scheitelwert A, ab; er ist, wenn 
nicht gerade die Scheitelwertabszisse a genau auf eine Merkmalklasse % fällt, größer als die größte 
Häufigkeit. 

Zu den Klassenmitten k lese man die von der Kurve bestimmten Häufigkeitszahlen ab; 
nur wenn Klassen- und Merkmalteilung übereinstimmen, sind dies die gegebenen Häufigkeits- - 
werte. Sie sind durch Division durch die Scheitelhäufigkeit h, in Prozenthäufigkeiten A%, um- 
zurechnen. Wohlgemerkt:sind diese Prozenthäufigkeiten A/h, nicht gleichbedeutend mit den 
relativen Häufigkeiten h/Ih = h/n. 

Diese h%,-Zahlen stellt man mit der k-Bezifferung der unteren Zunge an der h%-Leiter ein. 
War die vorliegende Häufigkeit schon einigermaßen genau einer Gaußschen Normalverteilung 
entsprechend, dann treffen alle A%, und die zugehörigen k immer gut aufeinander, und man 
liest ohne weiteres zu k%—=60,6, d.h. s=1 die Klassenbezifferung k, als Streuungsmaß ab. 
Dieses ist zum Schluß noch in Merkmalsmaßstab umzurechnen. Treffen jedoch die zusammen- 
gehörigen h%,- und k-Werte nicht glatt aufeinander und läßt sich auch nicht ein etwa grob 
fehlerhafter Zahlenwert als störend erkennen und auslassen oder sonst durch geringen Ausgleich 
eine für alle k passende Zungenstellung finden, dann muß der Mittelwert und die Streuung o 
entweder rechnerisch auf bekannte Weise ermittelt werden, oder man bedient sich mit Vorteil 
des zeichnerischen Verfahrens des Prozenthäufigkeitsnetzes in der folgend beschriebenen Ab- 


wandlung. 
& 3. Zeichnerisches Verfahren zum Auswerten gegebener Häufigkeiten 
g mit dem Lorenz-Maßstab des Sonderschiebers 
x Erscheint ein Ausgleichen der Prozenthäufigkeitszahlen allein mit der Klassennummer- 
# Teilung der Sonderschieberzunge nicht erfolgreich, so geht man vorteilhaft auf die zeichnerische 
EL Grundlage des Sonderschiebers zurück und zeichnet also die Lorenz-Gerade im Lorenzschen 
E > Prozenthäufigkeitsnetz. 
E Dazu benutzt man den an der oberen Stabkante des Sonderschiebers angebrachten Maß- 


stab h%, undträgt zueinerlinearen Merkmal-oder Klassenteilung die vorliegenden h%, als Ordinaten 
auf (s. Bild 1). War die Häufigkeitsverteilung nur annähernd eine Gaußsche, so streuen diese 
Punkte mehr oder weniger in bezug auf eine einzuzeichnende Ausgleichsgerade; dennoch ist diese 
meist ausreichend sicher zu finden. Durch diese ist dann der Mittelwert a als Argument des 
Durchgangs der Geraden dürch die Netzlinie 100%, abzulesen sowie der Zahlenwert der Streuung o ' 
als Differenz der Argumentwerte zu den Schnittpunkten mit den Netzlinien 100%, und +60,6%. 


Unter-Umständen erscheint ein zweites Aufzeichnen der A%Y,-Punkte und der Ausgleichs- 
geraden notwendig, wenn nämlich der Mittelwert a der Geraden vom Argument des einge- 
zeichneten 100%-Punktes allzusehr abweicht. Dann rechnet man die Prozentwerte in bezug auf 
die-zum Mittelwert a gehörende Häufigkeit nochmals aus und wiederholt das Verfahren, wobei 
sich Zeichenpunkt und Geradenschnittpunkt treffen werden. 2 


Die so gewonnene Lorenz-Gerade, deren Netz ja nicht groß eingezeichnet ist, läßt sich 
auch als Gerade der Summenprozente im Wahrscheinlichkeitsnetz betrachten. Dann ist der an 
der unteren Körperkante angebrachte Maßstab 3%, zum Ablesen der Summenprozente zu 
benutzen, ebenso zum Ablesen der Spannen 7 oder der „Gut“- und „‚Ausschuß“-Anteile zu 
beliebigen Merkmalswerten (s. auch Bild 1). 


4. Die Beurteilung von Häufigkeitszahlen 


In Anlehnung an die schon erwähnten Graphischen Tafeln von Koller wurde der Sonder- 
‚schieber mit den Läuferkurven ausgestattet, die in Zusammenwirken mit der Körperleiter h%, 
(hierbei jedoch als relative Häufigkeit h/nzu lesen) zum Vergleichder in zwei Reihen beobachteten 
Häufigkeiten bei gleichem Umfang der Reihen dienen (entsprechend der Tafel 5 von Kolle r). 
Die Kurven sind beziffert für n = 30, 50, 100, 300, 1000 und 5000. Sie liegen zu beiden Seiten 
des Ablesestriches, mit dem die zu vergleichende relative Häufigkeit eingestellt wird. Dabei sind 
die Kurven der einen Seite des Läufers für den „Iransponierungs‘“- und für den „Repräsen- 
tations“-Schluß bestimmt und so bezeichnet. Die linken Kurven dieser Seite zeigen die oberen 


NEE 
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3 o-Grenzen-des Schwankungsbereiches der kleineren Häufigkeit p,, die rechten die unteren 
Grenzen der größeren Häufigkeit Po: 


Grundfrage ist also: Es sollen zwei Beobachtungsreihen mit gleichem Umfang, aber ver- 
schiedenem Mittelwert daraufhin geprüft werden, ob sie trotzdem als gleichwertig betrachtet 
werden können oder ob ein mehr als zufälliger Güteunterschied begründet werden kann. Man 
sucht zur kleineren relativen Häufigkeit h/n= p,die obere Grenze Pu+-A,des Zufallsbereiches 24, 
oder zur größeren Häufigkeit p die untere Grenze 9,— 4, und bildet die Differenzen A, bzw. A 
Diese sind zwar meistens etwas unterschiedlich, besonders bei kleinen Häufigkeiten und kleinen n, 
das ändert jedoch nichts an der Schlußweise, daß es sich um einen gesicherten Unterschied der 
beiden Beobachtungsreihen nur dann handeln wird,. wenn %— Pu > A bzw. Ay. 


Beispiel 1°): Vergleich zweier Häufigkeiten aus gleichgroßen Proben 
Von je 100 Röhren waren 70 aus einem, 80 aus einem anderen Herstellungsverfahren 
einwandfrei. Ist ein Qualitätsunterschied gesichert ? 
70% als p, mit dem Läuferstrich eingestellt läßt bei n — 100 die Bereichsgrenze p,-+ 4, 
= 88% ablesen, der Zufallsbereich A„= 18% ; da 9, — p, = 80 — 70 = 10%, und also weniger 
als der Zufallsbereich ist, ist an der geringen Stückzahl kein Güteunterschied festzustellen. 


Wegen der logarithmischen Teilung der Sonderschieberleitern ist die Häufigkeit O0 nicht 


‚einzustellen noch abzulesen, ein Nullergebnis einer Stichprobe kann somit auf die besprochene 


Art nicht geprüft werden. Als Ausweg bietet sich die Benutzung der später noch zu besprechenden 
Leiter 9,, die zur Beurteilung eines Nullergebnisses beim Repräsentationsschluß dient. Durch 
Multiplikation ihrer Ergebniswerte mit einem Faktor von im Mittel 1,66 werden erfahrungsgemäß 
annähernd richtige Werte auch für ein Nullergebnis beim Vergleich gleichgroßer Massen gewonnen. 
Davon wird in Abschnitt 5 nochmals zu sprechen sein. Auch ein Beispiel wird dort gebrach 
werden. | 


Liegen nicht gleichgroße Kollektive vor, sondern schließt man von der beobachteten 
Häufigkeit einer Stichprobe mit n=n, auf die Wahrscheinlichkeit bei einer Stichprobe anderen 
Umfanges n,, so sind die wie vorher berechneten Grenzen des Zufallsbereiches im Verhältnis ns/n, 
umzurechnen. Dazu dient die kleine Leiter %s/n, über der linken Seite der oberen Zunge. Man 
stellt die Marke nz/n, auf das vorliegende Mengenverhältnis und liest zu der wie für den Fall 
gleichgroßer Kollektive berechneten Grenzdifferenz A, auf der Zunge eingestellt, den berichtigten 


Wert auf der Körperteilung ab: ‚‚Transponierungsschluß“. 


Beispiel 2: Vergleich zweier Häufigkeiten aus verschieden großen’ Proben. 

Von 50 Patienten sind auf die eine Art 25, d.h. 50%, auf eine andere Art von 20 Patienten 
18, also 90%, geheilt worden. Ist dieses Heilverfahren wirklich besser ? 

Zu 50% und n—50 liest man p, A, = 80% ab, der Zufallsbereich A = 30%. Da das 
Verhältnis n,/n, = 50/20 = 2,5 ist stellt man die Zunge mit ihrer linken Marke n,/n, auf den 
Teilstrich nz/n, = 2,5 und liest zum Zungenwert 30 das berichtigte A = 25% ab. Da die wirk- . 
liche Differenz 9 — p, = 40% wesentlich größer ist als diese 3 o-Bereichsgrenze mit nur 0,3% 
Überschreitungswahrscheinlichkeit, ist ein gut gesicherter Erfolgsunterschied festzustellen. 


Ist im besonderen das zweite Kollektiv sehr’groß, das. Verhältnis n,/n, daher als oo anzu- 
sehen, so entspricht das Vorgehen dem statistischen ‚‚Repräsentationsschluß“, der bei Koller 
mittels der Tafel 4 bearbeitet wird, Beim Sonderschieber ist es nur notwendig, die Marke n,/nı 
der Zunge auf den Teilstrich ns/n = &© einzustellen und im übrigen wie vorher besprochen 
vorzugehen. Für kleine n und extrem liegende Währscheinlichkeiten sind die so gefundenen 
Ergebnisse wegen der geringen Neigung der ‘die Körperleiter schneidenden Kurven ungenau, 
besonders bei dem provisorischen Sonderschieber-Modell, das Bild 2 zugrunde liegt. 


Beispiel 3: Ermittlung einer Grundwahrscheinlichkeit 
Wenn von 20 Kranken 6 geheilt wurden, ist die Heilungswahrscheinlichkeit wie hoch 
einzugrenzen ? | | 
Stichprobenhäufigkeit 30% und n=20 ergeben etwa (über n—30 hinaus extrapoliert) 


> die Bereichsgrenzen 76 und 0,3%, die Schwankungsbereiche sind A= + 46% und A=— 27%; 


"wegen n,/n = © berichtigen sich diese Zahlen auf + 32,5% und — 19%. Der Bereich der 
"Heilungswahrscheinlichkeit reicht also von 11 bis 62,5%. Genaueres läßt sich aus der kleinen 
Versuchszahl nicht entnehmen. | 


Die folgendeh Beispiele wurden mit der Absicht aus Koller entnommen, einen Vergleich der beiden 
Verfahren zu erleichtern. 
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Beispiel 4: Bestimmung des Stichprobenumfanges 

Neuartige Legierungen für Feinwiderstände sollen auf ihrer Bauchbarkeit geprüft werden, 
indem aus Draht von jeder Schmelze Probewiderstände gewickelt werden. Wieviele gleicher Art 
sind jeweils mindestens nötig, um ein einigermaßen gesichertes Ergebnis zu erhalten ? Reichen 
10 Muster oder auch schon 3 aus? 

Es ist recht undankbar, ohne gewisse Erfahrungen und Kenntnis der vermutlichen Güte- 
schwankungen solche Entscheidung treffen zu sollen. Wie hoch jedoch auch der Prozentsatz 
der brauchbaren Probewiderstände werden mag: ganz allgemein läßt sich an Hand der Tafeln 
von Koller im Voraus doch sagen, daß bei n= 10 der Streubereich für alle Häufigkeitsprozente 
rd. +50% beträgt. Zu n—=30 ist auf dem Sonderschieber abzulesen, daß die Streugrenze für die 
verschiedenen möglichen Stichprobenhäufigkeiten immer noch um +20% bis +30% beträgt. 
(Hierbei muß wiederum die mit Läuferstrich und Läuferkurve abzulesende Differenz mit dem 
Teilstrich n,/n,— ©o umgerechnet werden.) Deshalb wurden kleinere Probenzahlen als n—=30 
auf dem Sonderschieberläufer gar nicht erst aufgenommen. Es wäre Verschwendung von Arbeit, 
Zeit und Material, Stichproben geringeren Umfanges überhaupt zu bearbeiten. Selbst wenn 
die Überschreitungswahrscheinlichkeit von nur 3%/% bei dieser Fragestellung als übertrieben 
scharf angesehen und hier zur einfachen Streuungsgrenze übergegangen wird mit der Gewiß- 
heit, wenigstens in 2/, aller Fälle im erlaubten Bereich zu bleiben, wäre der Streubereich bei 
n—30 im Mittel noch + 8%. 

Wenn natürlich schon Erfahrungswerte über die im allgemeinen eingehaltene Gleichmäßig- 
keit der Schmelze oder der Drahtzieherei vorliegen, dann können exaktere Wahrscheinlichkeits- 
grenzen vorhergesagt werden und unter Umständen kleinere Probenzahlen ausreichend sein. 

Die umgekehrte Fragestellung, nämlich eine bekannte Grundwahrscheinlichkeit für sehr 
großen Umfang des Kollektives 1 an einer Stichprobe kleineren Umfanges nachzuprüfen, ‚‚Inklu- 
sionsschluß‘“, ist mittels des gleichen Vorgehens nur annähernd zu lösen. Für Grundwahrschein- 
lichkeiten etwa kleiner’als 15%, oder größer als 85%, sind außer bei recht großem n die um- 
gerechneten Streubereiche hier kleiner als im vorigen Fall des Repräsentationsschlusses. 


s Beispiel 5: Stichprobenkontrolle einer Grundwahrscheinlichkeit 


Eine Lieferung mit zu 20%, vereinbarten Stücken erster Wahl enthält unter 1000 nur 
150 erstklassige, also 15%. Ist dies noch eine erlaubte Abweichung ? 


E Zu 20% und n— 1000 liest man p+4A,=26% und damit A,—=6% ab, oder auch zu 

VE p— Au= 15%, also A,—=5%. Das Mengenverhältnis ist jedoch oo, also reduzieren sich die 

“ Grenzen auf + 4,3% und — 3,55%. Der Bereich der erlaubten Schwankung ist also 16,45%, 

| bis 24,3%, und die gefundene Abweichung um 5% auf 15% ist deutlich größer, also wohl vor- 
sätzlich herbeigeführt. 


Richtiger wäre ein Fehlerbereich von + 3,8%. Es ist ersichtlich, daß dadurch keine 
wesentliche Änderung des Ergebnisses eintritt. 

Sollte wirklich die mit dem Sonderschieber erreichbare, wenn auch nur grobe Annäherung 
an die Fehlerbereiche beim Inklusionsschluß für die praktische Verwendung nicht ausreichen 
so steht nichts im Wege, an Stelle der einen Schar von Läuferkurven noch eine zweite Schar 
anzubringen oder einen ganzen zweiten Läufer vorzusehen. Beim vorliegenden vorläufigen 
Modell des Sonderschiebers wurde die Rückseite des transparenten Läufers mit der Kurvenschar 
für den Inklusionsschluß versehen. Mit ihr ist die in Beispiel 5 erwähnte richtige Streubereich- 
Ablesung vorgenommen. 

Die Genauigkeit dieser Ablesungen ist wegen der kleinen Teilungslängen und der weit- 
abständigen Läuferkurven zwar wesentlich geringer als selbst bei den zum Vorbild genommenen 
Tafeln von Koller, der diese schon gegen den Vorwurf der Ungenauigkeit verteidigen mußte 
Andererseits ist der Sonderschieber durch seine Handlichkeit und die Vielseitigkeit seiner An- 
wendung sehr im Vorteil. Für viele Anwendungsgebiete der Praxis wird das weitaus wichtiger 
sein als eine zusätzliche Dezimalstelle. Und gerade die hier unübertroffen einfache Ablesung 
der Zufallsgrenzen wird den Sonderschieber im Gegensatz zu den-Kurventafeln und insbesondere 
den rechnerischen Verfahren in der Praxis weite Anwendung finden lassen, wenn es mehr auf 
schnelles Ergebnis und beste Entscheidung ankommt als auf einige Prozent an Genauigkeit. 


5. Der Fehlerbereich von Mittelwerten 
Die auf dem Sonderschieber angebrachte Teilung t dient zum Beantwort 
hiebe worten der Frage- 

stellung: Mit welcher Genauigkeit kann man den Mittelwert einer Stichprobe aus einer Bots 
Gesamtmenge als Näherungswert für den wahren, unerreichbaren Mittelwert dieser Gesamt- 


menge ansehen? Innerhalb welcher Grenzen kann die j : 
Elwert abweichen? ieser Stichprobenmittelwert vom wahren 
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Als Grenze des Fehlerbereiches ist wie bei allen diesen Betrachtungen die in der statisti- 
schen Praxis eingeführte 30-Grenze mit 0,3%, Überschreitungswahrscheinlichkeit zugrunde 
gelegt. Ohne weiteres gilt diese Grenze nur für große statistische Mengen mit n = 00; dann ist 

- der Fehlerbereich des Mittelwertes M vom dreifachen Wert des mittleren Fehlers des Mittel- 
wertes oy, begrenzt. Bei den in Wirklichkeit ja nur möglichen kleineren Stichproben verbreitert 
sich dieser Fehlerbereich, weil die Bestimmung von M und oy an kleinerem Zahlenmaterial 
selbst unsicher wird. Die Teilung i gibt nun für jede Zahl der Freiheitsgrade n —1 zwischen 
2 und 100 das Vielfache i des mittleren Fehlers o,, an, also den halben Fehlerbereich. Deshalb 

"ist an der Bezifferung der i das Merkzeichen + gesetzt, um daran zu erinnern, daß diese Fehler- 
grenze nach größeren und kleineren Werten zu gilt: Fehlerbereich=a-t- Oy- 

Die Begrenzung auf Stichprobengrößen n = 3 - - - 101 ist für die Anwendung ausreichend, 
weil kleinere Beobachtungszahlen sowieso zu keiner verwertbaren Aussage berechtigen und weil 


für größere Stichprobenzahlen der Faktor t asymptotisch von 3,08 auf 3,00 abnimmt; es lohnt - 


nicht, diese geringe Änderung in der Sonderschieber-Leiter ausdrücklich zu berücksichtigen. 


Beispiel 6: Fehlerbereich von Mittelwerten 


.. __An 88 Gewinde-Gut-Lehrdornen werden die Abweichungen vom Sollmaß gemessen und 
als Mittelwert M= 1,3875 u, als Streuung o—=1,18 u berechnet. Der mittlere Fehler des Mittel- 


wertes oy ist o/n — 0,1255u. Man liest daraufhin zu n—1=87 den Faktor t—=3,09 ab. Der 


£ Mittelwert der Gesamtmenge dieser Lehrdorne kann dann mit einer Überschreitungswahr- 
) _  ‚scheinlichkeit von nur 0,3%, (falls natürlich alle Fertigungsverhältnisse gleich bleiben) inner- 
halb des Bereiches von + 3,09 - 0,1255 = + 0,338 u liegen, also von 0,9995 bis 1,7755u Ab- 
weichung vom Sollwert als Mittelwert auftreten. Eine bestimmtere Angabe ist nur auf Grund 

einer größeren Stichprobe zu vertreten. 
Die Kenntnis und einfache Ermittlung dieses Fehlerbereiches ist besonders wichtig um 
zu verhindern, daß aus den Ergebnissen von oft zu kleinen Stichproben zu weitgehende und 

. dann unberechtigte Schlüsse auf die Gesamtmenge gezogen werden. 


6. Die Beurteilung eines Nullergebnisses 
Die Leiter p, des Sonderschiebers dient zur Bestimmung der oberen Grenze der Grund- 
_ wahrscheinlichkeit eines Ereignisses in einem großen Kollektiv, wenn eine Stichprobe an einer 
kleinen Anzahl von Beobachtungen die Häufigkeit 0 ergeben hat, dieses Ausbleiben jedoch 
augenscheinlich nur am geringen Umfang der Probe liegt. > j j 
Der Faktor p, ist das n-fache der oberen Grenze der Grundwahrscheinlichkeit, die als 
- untere Grenze eben das Nullergebnis hat. Die Probenanzahl n stellt man mit dem Läuferstrich 
' auf der logarithmischen Körperleiter ein und liest das zugehörige p,—= np ab und dividiert es 


durch n. Der Bereich für n läuft von 2 bis 100; für größere n nähert sich np asymptotisch dem . 


 Zahlenwert 5,9145. 


Bi: Beispiel 7: Treffwahrscheinlichkeit trotz Nullergebnis 
; Bei 5 Schüssen ‘wurde die Scheibe stets verfehlt; wie wahrscheinlich ist es, daß sie 
' auf Dauer doch einmal getroffen wird ? 


. ‚Schießen in höchstens 70%, aller Schüsse doch noch getroffen werden können. . 
Sa Wie schon in Abschnitt 4 vorweggenommen wurde, kann die Leiter p, auch zur annähern- 
E den Berechnung eines Nullergebnisses beim Vergleich gleichgroßer Massen benutzt werden. Die 
_ Zahlenwerte p—= »/n sind dann mit einem von n abhängigen Faktor zu ‚multiplizieren, eine 
- dazu dienende Teilung ist links auf der oberen Zunge mit n, bezeichnet, die zugehörige Marke 
mr Ale ; | 


2 en: x - Beispiel 8: Treffhäufigkeit trotz Nullergebnis bei gleichgroßen Proben 


Eh Br - - Wenn 30 Schüsse die Scheibe verfehlten :wieviele Treffer sind bei den nächsten 30 Schüssen 


erwarten? ey. Nr Be a 
NE Zu n = 30 gehört 9% = 5,37, p = 9/n = 17,9. Multipliziert mit 1,55, indem die Marke 0 
er auf den Teilstrich » — 30 eingestellt ist, wird die Trefferwahrscheinlichkeit 28%, die Treffer- 
zahl höchstens 11. re ee = Re \ 


2 Dieser Aufsatz über ein neues Hilfsmittel zur statistischen Häufigkeitskontrolle darf jedoch 
ee. keinesfalls abgeschlossen werden ohne den ausdrücklichen Hinweis, daß alle Verfeinerungen des 
 yissenschaftlichen Apparates letztlich nur dann berechtigt sind, wenn durch die daran gebildeten 


= 


& : Be Zu n—= 5 gehört np = 3,46; p = 3,46/5 — 0,692. Die Scheibe wird demnach bei emsigem 
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exakten Erfahrungen und Kenntnisse eine Stärkung und Sicherung des eigenen unbefangenen 
Urteils bewirkt wird. ne 

Man sollte auch nicht vergessen, daß die Gaußsche Normalform sich als Lehre für Häufig- 
keitsverteilungen nur deshalb empfiehlt und einführte, weil sie mathematisch die meisten ‚Ver- 
einfachungen bietet, nicht etwa weil sie den tatsächlichen Verhältnissen bestens entspräche. 
So läßt die Normalform mit ihrem asymptotischen Verlauf beider Äste auch große Abweichungen 
noch zu, wenn auch mit äußerst geringer Wahrscheinlichkeit. Dagegen ist es doch absolut un- 
möglich, nicht nur sehr unwahrscheinlich, daß z.B. beim Abstechen von dünnen Stiften be- 
stimmter Maßtoleranz aus Rundstahl selbst bei noch so großer Produktionsmenge einmal ein 
starker Bolzen aus dem Automaten fällt. Sollte es gelingen, wirklichkeitsnähere und plausiblere 
Fehlerfunktionen ebenso weit aufzubereiten und für die Praxis leicht anwendbar zu machen wie 
die Gaußsche, dann würde zum Teil begründete Abneigung gegen eine wirklichkeitsfremde 
Statistik ihren Sinn verlieren. 

Eingegangen am 25. Juli 1951. 


Graphische Lösung von gewöhnlichen Differentialgleichungen. 
mit nomographischen Hilfsmitteln”) 
Von W. Richter in Frankenthal 


Es wird ein allgemeines Verfahren angegeben, das gewisse gewöhnliche Differentialgleichungen ersier 
und höherer Ordnung mit nomographischen Hilfsmitteln zu lösen gestattet und welches die bekannten Ver- 
fahren für deren graphische Integration (Alferov, C’zuber, Heinrich, Mehmke, Sobotka, Willers) 
als Sonderfälle umfaßt. An Beispielen wird die praktische Anwendung gezeigt. Es wird auf Sonderfälle 
de. allgemeinen Verfahrens sowie auf die Möglichkeiten, nomographische Hüfsmittel noch weitgehender zu 
verwenden, hingewiesen. 


The author gives a general method for the solution of certain ordinary differential equations of the first 
en F2 or of a higher order, by aid of nomograms. T'he method includes the known graphical methods of integration 
R (Alferov, C'’zuber, Heinrich, Mehmke, Sobotka, Willers) for the equations of the special type as special cases. 

| The practical application is demonstrated by several examples. Special cases of the general method are 
explained, and the possibility for a far wider use of nomographic methods is pointed out. 


L’auteur donne une methode generale pour la solution de certains &quations differentielles ordinaires 
du premier ordre on d’ordre superieur, & l’aide de nomograms. La methode comprend les me£thodes bien 
commes de l’integration graphique (Alferov, Ozuber, Heinrich, Mehmke, Sobotka, Willers) de ces &qualtions 
comme cas special. L’application practique est dEmonströe par plusieurs exemples. Plusieurs cas sp6ciaux 
de la methode generale sont expliqu£s ei il est reuvoy& a la possibilit d’un usage plus &iendu des möthodes 
nomographiques. . 

Naerca oöımml MeTon pemeHun HeKOTOPHIX IPOCTEIX AHbhepeHLHassHHX YPaBHeHHÜ ITePBOTO 
H BTOPOIO NOPAAKOB HOMOTPAßHYeCKHMH CPEACTBAMH; HTOT METON CONEePKHT OÖINEHSBECTHEIE 
CHOCOOH HUX TPaduyeckoro HHTErPHpoBaHns (CHocoöst A,rheposa, UYydepa, Treipuxa, Mamke, 
Coöorrku, Bunnepca) Kak YacTHsie cııyyau. Ilpaxrtuyeckoe IIpuMeHeHHe IIOKA3aHO Ha IPHMepax. 
YkaskIBaMTCH YACTHEIE CHYYAH OÖINETO METONA, A TAKIKE BOSMOMKHOCTB ÖOIEE IIMHPOKOTO HCHONb- 
80BAHHA HOMOTPAhHYecKHUX CPEenCTB. 


Im folgenden sollen £, stets in’Richtung der «-Achse und n; in Richtung der y-Achse auf- 
zutragende Koordinaten bedeuten, welche als Funktionen von x und y aufzufassen sind. Für 
ihre Differenzen sollen die Symbole 


(&)i,r = 0%; — (=E bzw. n) 

und für zweireihige Determinanten sollen die Symbole 
NR, _|% ßi 
| (&, Pix = [%, ],. or & Bi 
verwendet werden. Wenn ein i. A. konstruktiv festzulegendes Element P (Gerade oder Punkt) 
in einem speziellen Fall durch die Konstruktion nicht bestimmt ist, soll von einer P-Unbe- 


stimmtheit gesprochen werden. Fallen jedoch zwei solcher Elemente P und Q zusammen, 
so soll das eine PQ-Koinzidenz heißen. 


1. Das Verfahren von Helmut Heinrich 


Das Verfahren von Helmut Heinrich!) zur graphischen Integration von gewöhnlichen 
Differentialgleichungen erster Ordnung kann folgendermaßen beschrieben werden. Der Dif- 
ferentialgleichung y’ = f(&, y) ist ein Richtungsfeld tg a=f(x, y) zugeordnet. Das Linienelement 


*) Über einen Teil der hier gegebenen Entwicklungen habe ich auf der T' f 
u. Mech. im März 1951 in Teen. berichtet. E Be Fe 


1) Mitt. techn. Inst. Staatl. Tung-Chi-Universität (Woosung-China) II, 1: (1935). 


en. 


“4 bezifferte Funktionsleitern) 
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(2; yi, tg %) des Punktes P; (&,%;) wird durch seinen uneigentlichen Punkt Pr, also durch 
seinen Schnittpunkt mit der uneigentlichen Geraden 2, festgelegt. Dies geschieht mit Hilfe 
der Ablesegeraden g; eines Fluchtliniennomogramms, das außer aus 2, noch aus zwei Feldern 


I ad (z, Y) 
N N (x, y) 


I fe &,(®, Yy) 

\ 72102 (%, Y) 
besteht. Das Verfahren gestat- 
tet die graphische Integration 


von Differentialgleichungen der 
Form 


und 


‚_ Ma 
ee) 


(Bild 1). Beide Felder I und II 
können in Sonderfällen in Kur- 
ven (mit den Werten von # bzw. 


oder in Punkte (Pole) ausarten. 
Dieses Verfahren kann verall- 
gemeinert werden. 


2. Eigentliche Gerade als 
Zapfenlinie 

Wird anstatt der im Ver- 
fahren von Helmut Heinrich 
ausgezeichneten uneigentlichen 
Geraden z, eine eigentliche Ge- 
rade z, die Zapfenlinie genannt 
werden soll, gewählt, so erhält 
man das in Bild 2 dargestellte 
Verfahren. 

Ist-die Zapfenlinie z paral- 
lel zur Ordinatenachse, so er- 
gibt sich eine graphische Lösung 
für die Differentialgleichungen 
der Form 


(Ons (Mu2 | 
a—dı Y-Ma 
NT Er F 2 
YT@0)- On " 

(a= const). 
Ist jedoch die Zapfenlinie z nicht parallel zur Ordinatenachse, so folgt dafür eine graphisch 
Lösung für Differentialgleichungen der Form 


(2 (Muz 


Bild 2. Eigentliche Gerade als Zapfenlinie 


2 — — m (En)u2 
ee (m,b= const). . . . (8) 
| % Mr |-Enın 


Die Anwendung des Verfahrens ist naheliegend, wenn eine eigentliche Gerade, welche dann als 
Zapfenlinie. z zu wählen ist, zu den Lösungen einer Differentialgleichung y’=f(x,y) gehört. 
Für die praktische Verwendung sind Sonderfälle besonders geeignet. Beide Kurvenfelder Iund II 
können in Kurven, speziell in Gerade (von denen eine die uneigentliche Gerade sein kann), 
oder in Punkte (von denen einer ein uneigentlicher Punkt-sein kann) ausarten. 


2,1. Beispiele a) Die Bernoullische Differentialgleichung 
9 ee, | 


0 
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in welcher A, B und (© Konstante sind, ist mit 
a=A, nel, Mes &=1LT;, 


ein Sonderfall von (2). Für ein spezielles Beispiel werde A=B=—4 und (=—2 gewählt, 


: Sry Fp 


RE 
+2" 2+4 


is 
25 X) +25 
Bild 4. Differentialgleichung y’-2 a1 +2) +&@—y) ("—y)=0 


A(— 1, — 1) folgen, entsprechen im Richtungsfeld Knotenpunkte. Die PP*-Koinzidenz B—1,1 
ist ‚Schnittpunkt der beiden zu den Lösungen der Differenzgleichung gehörenden a 2 
Bild 4 wurde außer der durch den Punkt P(3,2) gehenden Lösungskurve eine Reihe weiterer 


Lösungskurven gezeichnet. 
3. Kurve als Zapfenlinie 


Anstatt einer Geraden kann als Zapfenlinie z auch eine Kurve, deren Eleschune n=eo(£) 
sei, gewählt werden (Bild 5). Man erhält dann eine graphische Lösung für Differentialgleichungen 


der Form 


En a u — (&nus 
(#— 8) Hua - 


2) Über Singularitäten s. Abschn. 6. _ 7 
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N2—= 0 


also die Differentialgleichung 
Be. y 
Ye 
vorgelegt. Die Zapfenlinie, 
die ebenso wie die. Gerade 
€&=(0=—2 zu den Lösungen 
der Differentialgleichung ge- 
hört, ist parallel zur y-Achse 
im Abstand E&=4=—4. 
Das Kurvenfeld I artet in 
den Pol $ (—2, —4) und das 
Kurvenfeld II in die z-Achse 
mit ihrer Bezifferung aus. 
A(—4,0) ist eine PP*- 
Koinzidenz ?), ihm entspricht 
ein Knotenpunkt des Rich- 
tungsfeldes. In Bild 3 wurde 
außer anderen Lösungskur- 
ven jene durch den Punkt 
P (6, 5) gehende gezeichnet. 


b) Die Differential- 
gleichung 


y-222(142) 
+ @-9) (0 
hat die Form (3) für 
m—=—1l1, b=0, 
= — a, mn=Y% 
=Mm=E. 
Außer der Zapfenlinie „= —£ 
gehört noch die Gerade 
£=—1 zu den Lösungen der 
Differentialgleichung. Das 
Kurvenfeld I besteht aus. 
zwei zu den Koordinaten- 
achsen parallelen Parallel- 
strahlenbüscheln, und das 
Kurvenfeld II artet in die 
mit den Werten von & be- 
zifferte Gerade n=E£ aus. 
Den beiden P*-Unbestimmt- 
heiten, welche aus den P’P’’- 
Koinzidenzen O (0,0) und 
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worin & jene Funktion von & und 
% bedeutet, welche durch 


au) ey) 1 
&lYy) may) 1|=0 (a) 
5 ed) 1 
bestimmt ist. Für o(&)=m&-b 
erhält man den in Abschnitt 2 be- 
sprochenen Falleiner zur Ordinaten- 
achse schiefen Geraden zals Zapfen- 
linie. Auch hier sind die im Ab- 
schnitt 2 angeführten Sonderfälle 
von Bedeutung. 


3,1. Beispiele 
Die Differentialgleichung 


Bild 5. Kurve als Zapfenlinie 


entspricht für „=gd)=-2, 4 —=a, n=0, &—=0, = y sowohl (4) als auch (4a). Als 
Zapfenlinie tritt die Parabel 7 = & auf, während die Felder I bzw. II in die mit den Werten 
von & bzw. % bezifferte 
‚x- bzw. y-Achse ausarten. 
Für den Koordinatenur- 
sprung O(0,0) als P’P"- 
Koinzidenz folgt eine P*- 
Unbestimmtheit, der ein 
Knotenpunkt des Rich- 
tungsfeldes entspricht. In 
Bild6 wurden die durch 
P (10,10) gehenden beiden 
Lösungskurven gezeichnet. 
Die für die Differential- 
gleichung- kennzeichnen- 
den Bereiche des Rich- 
tungsfeldes sind aus dem 
. die Geraden g festlegenden 
Fluchtliniennomogramm 
anschaulich daran zu er- 
kennen, ob g die Zapfen- 
linie z in zweireellen Punk- 3 5 
ten schneidet, sie berührt Bild 6, Differentialgleichung ya y’ u =0 
oder sie nicht reell trifft. a © 


Mn 10 


4. Niehtausgeartete Feldtransformation 


Dem Verfahren von Helmut Heinrich und den daraus abgeleiteten, in den Abschnitten 2 
und 3 besprochenen Verallgemeinerungen ist gemeinsam, daß das zweidimensionale Feld der 
Punkte P(z,y) auf Punkte P* abgebildet wird, welche eine eindimensionale Gesamtheit (un- 
eigentliche Gerade, eigentliche Gerade, Kurve) erfüllen. Es liegt also eine ausgeartete Feldtrans- 
formation vor. Diese einschränkende Annahme soll nun aufgegeben werden. Es werde au 
gefordert, daß die Punkte Pf als Schnittpunkte zweier Geraden 9; — ( PePiund’, = (Pi =) 
festgelegt sind. h;, die zweite von ihnen, die an Stelle der Zapfenlinie 2 tritt, sei durch die Felder 
III (P/’) und IV (P;‘) auf dieselbe Art bestimmt wie bisher die Gerade = (Pi Pi) durch 
die beiden Felder I(P;) und II(P;’) (Bild 7). 


Nach diesem Verfahren können Differentialgleichungen der Form 


y_ 310, Ola HI). (& m134 
- all, ME + I, El 


- h . Z. angew. Math. Mech. 
124 Richter, Graphische Lösung von gewöhnlichen Differentialgleichungen 34.32 Nr. 4/5 April/Mai 1952 


i öst w i = —0, &=a, m=1 (a=const), so geht (5) 
graphisch gelöst werden. Ist speziell &3;—=a, 3 BR 
in (2) über. Ist aber speziell &=0, n—b, &—=1, n=b+m (b,m =const), so geht (5) in 
3) über. ; TIDE 
“ Wieder werden bei Anwendungen die in Abschnitt 2 angegebenen Sonderfälle für die 
Felder I, II, III und IV besondere Bedeutung haben. 


4.1. Beispiel. Die Differentialgleichung: 


‚_@=-WW+re—29Y) 

’TW-a)@—2z+y) = 
hat die Form (5), wenn &=z, =, ,=yW, n=y &=Yy n=%, K=m=l/k (k>0) 
gewählt wird. Die Felder I bzw. Il sind die mit den Werten von x bzw. % bezifferten Parabeln, 
deren gemeinsamer Scheitel der Koordinatenursprung 0. (0,0) und deren Achse die y-Achse 
bzw. die z-Achse ist. Das Kurvenfeld III besteht aus zwei zu den Koordinatenachsen parallelen 
Parallelstrahlenbüscheln. Das Feld IV artet in einen uneigentlichen Pol, den unendlichfernen 
Punkt der zu den Lösungen der vorgegebenen Differentialgleichung gehörenden Geraden 7=& 
aus. Durch ihn gehen alle Geraden h;. Dem Koordinatenursprung 0 (0,0) und dem Punkt A(4, 4) 
als P’P’’-Koinzidenzen entsprechen g-Unbestimmtheiten. In Bild 8 wurden einige Lösungs- 
kurven der Differentialgleichung darunter jene durch P(—5, 36) gehende gezeichnet. 


5. Erweiterungen 


Durch Ketten von Fluchtliniennomogrammen mit Zapfenlinien können auf verschiedene 
Weise alle angeführten Verfahren erweitert werden, wie das in Bild 9 für den in Abschnitt 2 
besprochenen Fall einer eigentlichen Geraden als Zapfenlinie z angedeutet ist. Für die prak- 
tischen Anwendungen werden allerdings solche umfangreiche Zeichenarbeit erfordernde Verfahren 
nur in ganz besonderen Fällen wirklich verwendet werden. 


5.1. Beispiel. Die Differentialgleichung: 
Pe yy—Ca) i 
(—4A)(y—Ca)+ B(y— Da) 

in welcher A, B, C und D Konstante sind, ist mit m=b=0, Kein=322; 
.De—y UE—Y 
| Aa SE SEE 
ein Sonderfall von (3). Für ein spezielles Beispiel were A=0, B=—2, C=1D=-—1 
gewählt. also die Differentialgleichung 


wa) 


(k+ 0) 


’TayW-)—2(e+y) 
gelöst. Die x-Achse, die ebenso wie die Gerade n— — € zu den Lösungen der Differential- 


gleichung gehört, wird Zapfenlinie. Das Kurvenfeld Iartet in den Pol 8 (0,2) aus und das Kurven- 
feld II geht in die uneigentliche Gerade als Doppelleiter über. Auf ihr können die Punkte P‘/ 
durch ein Nomogramm festgelegt werden, das außer der uneigentlichen Geraden die beiden 
Leiten „=&=x (beziffert mit den Werten von x) und n„=—£=y (beziffert mit den 
Werten von y) enthält. Der Koordinatenursprung stellt eine P‘/-Unbestimmtheit, also eine 
P*-Unbestimmtheit dar, der Punkt A (— 2,0) ist eine P P*-Koinzidenz. Beide Punkte sind 
Knotenpunkte des Richtungsfeldes. In Bild 10 sind einige Lösungskurven der Differential- 
gleichung, darunter jene durch P (—1, —4) eingezeichnet. 


6. Singularitäten 


Das durch die Differentialgleichung y’—=f(&,y) bestimmte Linienfeld hat Singularitäten 
an jenen Stellen, in denen y’—=tg& unbestimmt ist. Bei den besprochenen Verfahren gibt es 
dafür zwei Möglichkeiten. Wenn P;=P? ist, soll von einer P P*-Koinzidenz gesprochen werden. 
Wenn jedoch Pf selbst unbestimmt ist, so soll das eine P*-Unbestimmtheit heißen. Bei den 
in den Abschnitten 2 und 3 besprochenen Verfahren können nur die Punkte der Zapfenlinie 2 
solche P P*-Koinzidenzen darstellen. Ist eine eigentliche Gerade Zapfenlinie (Abschnitt 2), so 
folgt aus 9,=z (also aus jeder gz-Koinzidenz), aus jeder P’-Unbestimmtheit (die. beiden mit «; 
bzw. y; bezifferten Kurven des Feldes I sind identisch), aus jeder P’’-Unbestimmtheit (die 
beiden mit &; bzw. %; bezifferten Kurven des Feldes II sind identisch) und aus P}=Pi’ (also 
aus jeder P’ P'-Koinzidenz) eine P*-Unbestimmtheit. Dasselbe gilt für den Fall einer Kurve 
„als Zapfenlinie (Abschnitt 3), wobei allerdings eine 92-Koinzidenz nur dann auftreten kann, wenn 
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eine Gerade als Teil der Zapfenlinie existiert. Für den Fall der nichtausgearteten Feldtrans- 
formation (Abschnitt 4) wird anstatt der gz-Koinzidenz möglicherweise eine gh-Koinzidenz 
auftreten und es sind auch P’’-Unbestimmtheiten, PIV-Unbestimmtheiten und P'" P!Y-Koinzi- 
denzen denkbar. 

Jeder Singularität einer Differentialgleichung muß bei Anwendung der besprochenen Ver- 
fahren eine der angeführten Koinzidenzen oder Unbestimmtheiten entsprechen. Die Umkehrung 
gilt jedoch nicht. 


S 7 
ta 8 We, 
% „8 
SS 
SI 
||s 
8 ja 
+|| 
SE 
38 
I! 


y= 


Bild 8, Differentialgleichung 


Bild 7. Nichtausgeartete Feldtransformation 


y(y — x) 
z(y —- 2) — 2(& + y) 


Bild 10, Differentialgleichung y’ = 


Bild 9. Kette von Fluchtliniennomogrammen | 
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7. Differentialgleiehungen zweiter Ordnung 


BI), EHE ” 
ale), M)lsı +1), &n)lsi 


definiert? ®3 4 als Funktion von a, b und jenen Veränderlichen, von denen $; und n; (i 1,208: 4) 
abhängen. Die Gerade, welche den Punkt P(a,b) mit dem Schnittpunkt P* der beiden Geraden 


9=[Pı(&» m); P2(&2 N)] 
/ und 
h=[Pz(&,n); Pı(&u 7)] 
verbindet, schließt mit der &- 
Achse den Winkel & ein, der 
durch tga—=3®3;; bestimmt 
ist, 

Die im Abschnitt 4 bespro- 
chene Methode der nichtausge- 
arteten Feldtransformation ist 
auf Differentialgleichungen von 
der Form (5) anwendbar, also 

auf Differentialgleichungen 
’ Y 12 
Er D3;. 
Fürdie Erweiterung auf Differen- 
tialgleichungen |zweiter Ord- 
nung?)gibtes'zwei Möglichkeiten. 

7.1. Eine Differentialglei- 
chung zweiter Ordnung habe die 
Form 


b a la FE 
PT 


Bild 11, Differentialgleichung zweiter Ordnung von der Form y’’ = “03 Pr. dy' BT, fr) 12 7 
Yan dx sa; REF ( ) 
worin i 
ee op); 
= m (& P) 
i=1,2,'3, 4 


Funktionen von & und ß 
bedeuten, und & und ß für 
jedes % jeweils zwei Verän- 
derliche aus dem Tripel «, 
y, y' sind. Zu den An- 
fangsbedingungen %,, %o, %0 
erhält man (Bild 11) die 
Punkte Pi, P3_o, Ps, 
Pi_,. Die Geraden 
90 = [Pı-0, " Ps_o] 
und 
ho = [P3-0,. Pi-0] 
bestimmen den Schnitt- 
, punkt Po*. Seine Verbin- 
Bild 12, Differentialgleichung zweiter Ordnung von der Form y’’= y’ ;E o1? dungsgerade pp mit P% 
ER “ } s %, %Y,) legt für einen will- 
kürlich gewählten kleinen Schritt e= x, — z,den Punkt P} (z,, yı) = ne dessen Ordinate yj 
der dem Punkt P, (&, Yo) benachbarte Punkt P, (x), yı) bestimmt werden kann. 


7.2. Die Differentialgleichung: 


a N ee N 


. ,°) Siehe auch: Helmut Heinrich: Deutsche Mathem. 3 _ wi ie 
Ing.-Arch. 8 (1937), 10 (1939), 11 (1940), 13 (1943). N Wh ee 
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in der fürdie Funktionen &; und n; dasselbe wiein 7.1 gilt, ergibt wegen = I_ "die in Bild 12 


r 


dargestellte graphische Methode. Für einen willkürlich gewählten engen Schritt e= y, — % 
folgt zu dem auf der Geraden p’ liegenden Punkt Pi (y,, yı) der dem Punkte P, (20, %) benach- 
barte Punkt P, (,, Yı). 

7.3. Beispiel. Die Differentialgleichung: 


7 _Y (3 TE Ye y) 


J (2% — y) 
entspricht der Form 7.1 für 
1 
R H=—% a Na 3% = 0% 
ir! k—0 N7=0 nmn=—Y. 


Das Feld I besteht aus zwei zu den Koordinatenachsen parallelen Parallelstrahlenbüscheln, das 
Feld II geht in den uneigentlichen Punkt der Geraden n=£ als Pol über, so daß die Geraden g; 


ya? 


VER —y—y)) 
(2% — Y) 


Bild 13. Differentialgleichung y’= 


—£ paralleles Parallelstrahlenbüschel bilden. Die Felder III und IV sind die mit den 
ne bzw. y' bezifferten Koordinatenachsen. Für die Anfangsbedingungen 


1; y=)d, yp=—2 
wurde in Bild 13 die Lösungskurve gezeichnet. 


8. Differentialgleichungen dritter Ordnung 
Für Differentialgleichungen dritter Ordnung sollen die Veränderlichen &, ßin den Funktionen 
| =, mem, dee 
für jedes % jeweils zwei Veränderliche aus dem Quadrupel x, Y, y’, y’’ bedeuten. 
Be, Für die Differentialgleichung: | 
i v"=Y6 | 
kann das in Bild 14 därgestellte Verfahren verwendet werden. Den Anfangsbedingungen %0, Yo» 


Pe Porbe ich wie in 7.1 findet man auf p’ zue=n—% 
4, y6 entsprechen die Punkte P,, Po, Po. Ähnlich wie in 7. 
| Punkt Pi zu’ welchem Pi und P, bestimmt werden können. 
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führt die Differentialgleichung 


day" 
dy 


Wegen 


8.2 


VL TE N 12 
Y —Y 5 D34 


auf das in Bild 15 gezeichnete Verfahren. 


1. 


grrr 


0A) u, + (R— A) RS, Bunyoporätepgusropa "ZT PIIE zı®,n ‚%=..A wıog 1op uoa SunupiO 1OpLıp Sunyops[äferyuszsirtet "ST PIIT 


75, .R 
sI ” A „R= ,,,f% wIiog Ip uoA ZunupıQ I99ylıp Bunyorojäferyusasjfidt "9T PIIA 


eh. / f 5 g x 
KLRanT ‚„‚f wIog 19p uoA ZunupIQ Iaylıp Sunyprajäferzusssg lad "FLPIT 


Dioaene 
ee, 


x 


x 
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d 2) 
8:3. y = m ergibt für die Differentialgleichung: 


2 


rn rn y! 12 
ee D34 
das in Bild 16 dargestellte Verfahren. 


8.4. Beispiel. Die Differentialgleichung: 


yayly —Y)+yY'y?(a+y')—=0 
hat mit 
=Yyı (0), =), » u=y 
NM=0, Nn=Y', Na ee 0 u) 
die Form 8.2. Alle vier Kurvenfelder des Nomogramms arten in die mit je einer der Veränder- 
lichen &, y, y' und y’' bezifferten Leitern auf den Koordinatenachsen aus. 


Für die Anfangsbedingungen & = 10, „=5,; y=—1, y = —.0,5 wurde in Bild 17 
die Lösungskurve gezeichnet. 


9. Erweiterungen 


Das Verfahren kann entsprechend auf gewisse Differentialgleichungen von höherer als 
dritter Ordnung erweitert werden. Auch spezielle Systeme von Differentialgleichungen können 
auf diese Weise behandelt werden. 


9.1. Beispiel. 


‚_?@—-Yy)—y2 
2 (2 — 8) — x2 

Bee: 2 (Y Es 2) E 
2.(y— a) + ay 

Die erste Gleichung 


dieses Systems er- 
hält man aus 


mit 
H=n- 
e 2 
& — gl 
,=Mm- 
Na 
N %s 
die zweite aus 
r z 12 
2» = a#s314 
mit 
& u 
ZU Pe 0 
u=0 
N {0} = . #Z dr 
2 —y) — y2 h 2:W—2) 
E =n—Y Blid 18, Differentialgleichungspaar y’= nn rzl= zw —o+ ey 
SEERgS 2 


In Bild 18 wurde für die Anfangsbedingungen % =—10, „= 4, 2, = 10 das Paar der Lösungs- 
kurven gezeichnet. _ 


P2 


Eingegangen am 20. Juli 1951. : 
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Die Ermüdung der Metalle*) 
Ergebnisse der von der EMPA von 1925—1950 durchgeführten Untersuchungen 
von M. Ros in Zürich 


Körper ohne und mit Kerbe. Spannungszustand homogen und heterogen, ein- und mehrachsig. Ein- 
heitliche Fassung der Bruchtheorie für statische und Ermüdungs- Beanspruchung. 


The author considers the fatigue of metals. He examined solids under various conditions: with and 
without notches; homogenous and inhomogenous state of stress, uni-awial as well as multi-azial. He gives 
a uniform version of the theory of breaking for static and dynamic stress. 

D’auteur discute la faligue des me&laux. Il a examine des solides sans et avec enlaille, d’&tat de tension 
homogene et het£rogene, d’un- ou plusieurs awes; il donne une formulation uniforme de la th£orie de cassure 
pour pression stalique et dynamique. 

Tera c Haceykoi HU Öe3 HaceykH. ÜnHoponHoe H P&3HOPOAHOe, O]HOOCHOe H MHOTOOCHO® 
COCTOAHUA HanpsokeHus. Ernnunasa BepcHha TeopkHu paspkIBa Aula HATpyskuH CTATHyecKof H 

° Harpysku Ha YTOMJIeHHe. 


EEE 
>. 


Voraussetzungen. Werkstoff praktisch fehlerfrei, quasi-homogen und quasi-isotrop, 
= voruntersucht und sorgfältig ausgewählt. Verhalten zähe. Verformungsmechanismus vorwiegend 
auf Gleitungen beruhend. Versuchstemperatur normal. _ 


I. Statische Anstrengung 
y Verlauf des Zug-Dehnungs- und Druck-Stauchungs-Diagramm gleich (Bild 3). 
Festigkeitsversuche mit quasiisotropen Metallen bei ein- und zweiachsigen homogenen 
Spannungszuständen haben im Gebiet der Zeitfestigkeit nach nur wenigen Lastspielen — Grenz- 
fall: statische Beanspruchung — gezeigt, daß im Falle zähen Verhaltens mit reinem Gleitungs- 
bruch die Guest-Mohrsche Hypothese erfüllt ist: 


1, = 9 _ konstant (Bild 1 und 2). 
g Demgegenüber gilt für den Fall spröden Verhaltens mit reinem Trennungsbruch, der 
Be auf der Zugseite die Mohrsche Hüllkurve begrenzt, die Spannungshypothese, wonach die wirk- 
=s liche, größte Zug-Hauptspannung o, im Bruchzustand, in sehr weiten Grenzen, unabhängig von 
Be: den übrigen Spannungskomponenten ist: 
B | Oyaz = 03 — konstant (Bild 1 und 2). 


‘ Für die Fließgrenze, als Grenzzustand der statischen Anstrengung, gilt die Theorie 
der resultierenden Schubspannung in der Oktaederebene des Hauptspannungs-Würfel- 
elementes 


v2 
 — ar Vo +05 -+ 08 — 0,05 — 0,05 — 030; 


bzw. die Vergleichspannung 
Ga —Voj + 05 + 05 — 0,05 — 0,05 — 0301 


welcher die resultierende Vergleichsdehnung 
Er a ed Ze Fe a NER 
Di; =Ö6, == 3 ; ö + 0 7 03 — 6,6, AG 0,0, m; 0,0, 


zugeordnet ist (Bild3und 4). Diese von der EMPA verallgemeinerte Theorie von 
Mohr ist mit der Gestaltänderungs-Hypothese identisch (Bild 5). 


Bei gekerbten Körpern (Bild 6-8) stellt sich zwischen den statisch höher angestreng- 
ten Stellen und den weniger beanspruchten Gebieten ein derartiger Spannungsausgleich ein, 
‚daß sich der Einfluß der Spannungsspitzen verliert; der gekerbte Zugkörper (Bild 7) weist zu- 
mindest die gleiche, infolge behinderter Querdehnung überwiegend eine etwas höhere spezifische 
Zugfestigkeit auf wie der nicht gekerbte, volle, einachsig beanspruchte Zugstab (Bild 8). 


aM *) M.Ros u. A. Eichinger: Die Bruchgefahr fester Körper bei ruhender — statischer — Beanspru- 
chung. EMPA-Bericht Nr. 172. Zürich, September 1949. — Die Bruchgefahr fester Körper bei AR 
Be ernehung 7 as ; Sera Sach Nr. 173. Zürich, September 1950. 
. Roß: Fatigue des Mötaux. Rev. Mötallurg. Paris, Mai/Juin 1947. — La Fatigue d ud 
Rev. Mö6tallurg. Paris, Novembre 1948. — La Fatigue des Metaux. Riviete La ee 


No. 1, gennaio, 1950. — Static Failure and Fatigue of Steels. i i iati 
re rer gue of Steels. Sheet and Strip Metal Users, Technical Association. 
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II. Erschöpfung durch 
Ermüdung 


1. Versuchskörper 
ohne Kerbe 
Spannungsverteilung 

homogen. 

Die Grenzwerte der 
versuchstechnisch ermit- 
telten einachsigen Er- 
müdungsfestigkeiten 
auf Zug 


der Wechselfestigkeit 07 
Ursprungsfestigkeit 0; 
Schwellfestigkeit 6,, 
(Bild 9) für verschiedene 
normale und hochwertige 
Kohlenstoff- als auch le- 


‚gierte Stähle mit Zug- 
festigkeiten 


ßz=03=40—115kg/mm?, 


sind dem Bild 10 zu ent- 
nehmen. 

Die Spannungs-Am- 
plitude o,—o, nimmt von 
der wahren statischen 
Druckfestigkeit gegen die 
wahre statische Zugfestig- 
keit (A/ B=1) ab (Bild 11). 

Handelt es sich dem- 
gegenüber um oft wieder- 
holte mehrachsige Be- 
anspruchung, so ist zur 
Frage der Ermüdungs- 
Bruchgefahr bei mehr- 
achsigen homogenen Span- 
nungszuständen festzustel- 
len, daß die Coulomb- 
Mohrsche Hypothese — 
wenigstens für die Ur- 
sprungsfestigkeit o, : — 


0 
nahezu genau erfüllt ist 
(Bild 12, 13 u. 1). Dieses 
Ergebnis wurde an Fluß- 
stahl, Stahlguß, Schweiß- 
gut, Aluminium und seinen 
Legierungen gewonnen 
(Bild 13). Lediglich im 
Falle zweiachsigen Zuges 
mit 0,03 ist sicherheits- 
haiber mit etwa 90% 
der einachsigen Zug-Ur- 
sprungsfestigkeit zu rech- 
nen. 
Die wahre Druck- 
Ursprungsfestigkeit ist bei 
glatten Stäben aus Fluß- 
stählen mindestens das 1,6- 


fache der Zug-Ursprungs- 


festigkeit (Bild 13 u. 14), 
a 5 E 9* 


Druck 


Olängs 


Bild 1. Grenzzustände der statischen und Ermüdungsbeanspruchungen für 


verschiedene Verhältniswerte' von = und 


Oquer 


5 . Spannungszustände 
Tü 


NIS 
zweiachsig, die dritte Hauptspannung ist gleich 0. 


£ i FA 
Torsion miF Innendruck 


0XWu® 
7% 


eh 
[AM 


L® 


% Gl (rascher laslabfalı, 
6 Zug h angs X 0 eher niedriger) 


era sr 
dl ER 
% 
% 
pas 


N 


\ 


ggf a 
\ zn | 
2 Hr URELITERIEO) 


© Zug \quer 
Tr (geschichleh),” 


7 


7% 


bl (schr&g, nicht längs u.guer) 


% 
% Z 
Torsion voll (unter 45°) 
92 
Ws GI-Gleilungsoruch 
Drucklöngs 7 Ir-Irennungsbruch 
BASE: 


Bild 2. Flußstahl „SMI“. Zweischsige Spannungszustände an der Fließgrenze, beim 


Erreichen 


der Höchstlast und für die Bruchlast. Gemessene Werte. 
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' ' ä j der Verhältniszahl d/b (Bild 14). 
h en bei gelochten Stäben nur das 1,1—1,2fache, je nach 
Dieses Verhalten der gelochten Stäbe ist auf das während der wiederholten Druckbeanspruchung 


Stahl „N” 


Schubspannung in der Oktaeder-Ebene : 
Vergleichsspannung : 
Vergleichsdehnung : 


32 d5 -4 6, + d, +63 -6,6,-6,63-0364 


Normalspannun 


der Oktaeder-Ebene. pa 


0,0 ‚03 = Hauplspannungen 
d,, 6, , ö, = plastische Haupt 
dehnungen 


Verlängerung 
40 30 20 105 


Bild3. Normaler Flußstahl ‚‚St N“. Zug-Dehnungs- und Druck- Stauchungs-Diagramme weisen gleichen Ver- 
lauf auf. Schubspannung und Normalspannung der Oktaederebene. Vergleichsspannung. Vergleichsdehnung. 


_ 


sich ausbildende ungünstige System von Eigenspannungen zurückzuführen. Aus diesem Grunde 
darf die zulässige Druckbeanspruchung nur wenig höher als bei der Zugbeanspruchung ge- 
wählt werden. Die reine Schub- 
Ursprungsfestigkeit beträgt bei 
Flußstahl, Stahlguß, Reinalumi- 
nium und Avional etwa 70% von 
der Zug-Ursprungsfestigkeit (Bild 
13), bei Schweißgut sogar etwas 
-über 80 %. 


Bei technischen Werkstof- 
fen sind noch in der Regel die 
Festigkeits- und Verformungs- 
eigenschaften von der Orientie- 
rung der Hauptspannungsachsen 
0,0505 zu den Hauptachsen der 


2 = Vonsuiche EMPA 
zug 
a Fe 


Huber- v Mises 
- Hencky 


-@—[A Druck f — Zu -{Al-+% 


£% 
Y IN 
Zn 


4 


Guest-Mohr 


Guest-Mohr ; Tax = 2 z 9 _ konstant. 


Huber - u. Mises - Hencky: Oye5=/0, +0, -0;- 
(bei 0,=0) 

Bild 4. Normaler Flußstahl ‚‚St N“. Zweiachsige Beanspruchungen. Grenz- 

zustände an der Fließgrenze. Ergebnisse der EMPA-Versuche. Theoretische 


Werte entsprechend der Gestaltänderungs-Hypothese von Huber- 
v. Mises-Hencky und der Theorie von Guest- Mohr, 


0, = konstant, 


Anisotropie des Materials mehr 
oder weniger stark abhängig. _ 


Bei Stählen ist die Zug- 
Ermüdungsfestigkeit in der Quer- 
richtung — selbst bei relativ ein- 
wandfreien Schmiedestücken — 
um etwa 20% geringer als in der 
Längsrichtung und Bleche weisen 
bekannterweise winkelrecht zur 
Oberfläche oft nur ganz geringe 
Dehnbarkeit auf, wodurch ein vor- 
zeitiges Versagen im Falle einer 
Zugbeanspruchung in dieser Rich- 
tung verursacht werden kann. Dies 
ist auch der Hauptgrund dafür, daß 
sich der dreiachsige Zug so gefähr- 
lich auswirken kann. 


m 


HRS 


Z.angew. Math. Mech, 
Bd. 32 Nr. 4/5 April/Mai1952 


Ros, Die Ermüdung der Metalle 


2. Versuchskörper mit Kerbe 


Spannungsverteilung hete- 
rogen. 


Im Falle gekerbter Stäbe gel- 
ten für ein bestimmtes Element 
das/dy/d2 der heterogenen Span- 
nungsverteilung grundsätzlich die 
gleichen Betrachtungen wie für ein 
Würfelelement eines homogen be- 
anspruchten Körpers ohne Kerbe. 
Die Betrachtungsweise führt aber 
zu verwickelteren und weniger über- 
sichtlichen Verhältnissen als die 
einheitliche und gleichzeitige Er- 
forschung des gesamten Spannungs- 
bzw. Anstrengungsverlaufes im am 
meisten gefährdeten Querschnitt 
des auf Ermüdung beanspruchten 
. Körpers. 

Zweckdienlicher und klarer 
ist daher der nachfolgende Weg: 


Bei heterogener Bean- 
spruchung in gelochten und ge- 
kerbten Stäben wie auch an Quer- 
schnittsübergängen (Bild 15—18) 
darf an der Stelle der Spannungs- 
spitze, gemäß den Versuchsergeb- 
nissen im Pulsator, mit einer hö- 
heren rechnerischen zulässi- 
gen Beanspruchung auf Zug 
gerechnet werden als in- einem 
glatten, homogen beanspruchten 
Versuchsstab (Bild 19 u. 20). Hier- 
bei handelt es sich um einen rech- 
nerischen ‚‚elastischen Maßstab“ 
und nicht um’eine wahre Span- 
nungsgröße. Bei einer .'Beanspru- 
chung durch Zug nimmt im Er- 
 müdungszustand die rechnerische 
Spannungsspitze O,.. mit anstei- 
gender Formzahl &, zu, und zwar 
‘für die füblichen Flußstähle und 
Aluminiumlegierungen von o,..— 1 
sure lzbisiehwao,.=17 für 
%—=6 (Bild 19 u. 20). Es deutet 
dies daraufhin, daß sich ein Span- 
nungsabbau in der Zone der 
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Statische Anstrengung und Bruchgefshr. 


4. Hypothese der 
konstanten Gesfalt- 
änderungsenergie 


Huber -Mises - Hencky 
Mohr - EMPA 


0,0, ,05= Haupkspannungen 


0 =V0,*.0,405-0,0,-0.0,-0;0, 


Von der EMPA verallgemeinerte MohrscheTheore 


a Gleitungsbruch 


Theorie Guest-Mohr = Tas; = konsk. 


Bild 5. Gestaltänderungs-Hypothese identisch mit der von der EMPA 

verallgemeinerten Theorie von Mohr, gültig auch für die plastischen Ver- 

formungen bis kurz vor Bruch. Für den Gleitungsbruch selbst trifft die 
Theorie von Guest-Mohr zu. 


4 3 2 


I 
= 


0x Re 


Bild 6. Einachsig beanspruchter gelochter Stab. Spannungsverteilung. 


P . d d 
Bild 7. Normaler Flußstahl „St N“. Gelochte Stäbe, — = 0,016 — 0,06 — 0,257. 
“ Statische Bruchflächen. 


60 


Festig keiten 


Zugfestigkeit B, 


Verhaltniszahl 


Ren] 
>) 


2,0 


1,0 


Loch $ 0,8mm 3,0mm 36mm 
9= & = 0,046 2 - 0,06 0,257 


Bild 8. Normaler Flußstahl ‚‚StN“. Gelochte Stäbe — siehe Abb.7 — Zug- 
festigkeit Bz. Zug- Ursprungsfestigkeit eO& = 9, gp Spannungsgefälle (Gradient). 


Ermüdungsfestigkeit 


Sinusformige Wechselbeanspruchung 


Fliessgrenze 


a 
Ursprungs - 
5 Festigkeit 


A festigkeit 


0, = untere Grenze 
z 0, = obere Grenze 


‚Bild 9. Ermüdungsfestigkeitsdiagramm. Beanspr h 
Schematische Darstellung. Veen einachsig, 


- HE: 
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Spannungsspitze vollzieht, 
wodurch sich die Beanspruchung 
aus dem Gebiet der Zug-Ursprungs- 
bzw. Wechselbeanspruchung gegen 
das Gebiet der Druck-Ursprungs- 
beanspruchung des Ermüdungsdia- 
grammes verlagert (Bild 21—23)., 
so daß trotz einer Dehnungsampli- 
tude, die praktisch der rechneri- 
schen Spannungsspitze entspricht, 
der gekerbte Stab diese bis 1,7fache 
bzw. 6,4fache rechnerische Span- 
nungsamplitude zu ertragen ver- 
mag (Bild 19). 

Im Falle der Biegebeanspru- 
chung ist bei Rechteckstäben 


[07 
%—=1,5 MAL 
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ne N gg Tue HN 2 EEE a ARE Er > BEE EEE 2 EEE En 
kalmm® 7,2 Gelochte Stabe 
70 Uusp. > Spannungsspitze 


M/W 
und bei Rundstäben 
Bu Ta 


zu setzen (Bild24). Die Formzahl a, 
wird nämlich so definiert, daß sie 
die vom Zustand des Erreichens 
der Grenzspannung an der Stelle 
der Spannungsspitze Oyar—= Ak * On 
bis zur völligen Ausnützung der 
im engsten Querschnitt vorhan- 
denen Reserve angibt, und zwar 
in der Form: 


Belastungsgröße in 


ee Belastungsgröße bei rein 


vollplastischem Zustand 
elastischem Verhalten 


worin die Belastungsgröße P,M gieo.» 
Mrors., 9 usw. sein kann. 

Dies galt im Falle einachsiger 
Beanspruchung an der Stelle der 
Spannungsspitze. Treten demgegen- 
über an den Stellen der höchsten 
Beanspruchung mehrachsige Span- 
nungszustände auf, dann muß die 
Grenzlinie der Ermüdungszu- 
stände für das jeweils vorliegende 
Verhältnis A/B (Bild 19) herange- 
zogen werden, wozu noch die Form- 
zahl &, unter Berücksichtigung der 
Anstrengung bei mehrachsigen 


Spannungszuständen, 
o 
Oo en eY’ mar, res 
eIn res 


. zu ermitteln ist. Auf diese Weise 


kann der Konstrukteur mit Hilfe 
der in den Bildern 19, 20 und 24 
gegebenen &% —0,,,-Beziehung den 
Sicherheitsgrad auch bei sehr ver- 
wickeltem Spannungsverlauf im 
Körper mit einer für die Praxis aus- 
reichenden Genauigkeit beurteilen. 
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Legierte und nichtlegierte Stähle. Cr -Ni-Mo- Stahl 


%02 = 100 (1) 


Ermüdungsfestigkeiter 
im kq /mm? 
Stäbe voll, nicht gelacht. 


Ni- Stahl 


%=58(1) 


Cr-Mo-Stahl 
0.-50.@) 


ötratische Festigkeitswerte 
Mittelwerte 


Legierte ötahle: 


ötahlart Bz 05 Aro 

Cr-N-Mo 15.100 0 
Ni 75 58 14 
h Cr-Mo n 50 #6 

Nicht bezw. schwach legiert: 

Stahlark: Br 05 Ay 
Bz = Zugfestigkeit in kg/mm?® ‚sT7o” a 32 
05 = Streckgrenze in kg/mm? ‚9T32° Var dere 20 
No“ Bruchdehrung in %o „ST44" Am = 30 72 
‚STN’ 08 35 

Toleranz 10% 10% 20% 


Bild 10, Unlegierte und legierte Konstruktionsstähle. Diagramme derjeinachsigen Ermüdungsfestigkeiten 
jängs, unlegierte Stähle f, — 40, 47, 55 und 75 kg/mm?®. Legierte Stähle , = 70, 75 und 115 kg/mm?, 


Reissfestigkeit bzw. 
Zug -Kriechfesti gkeit 


„R=80 kglmm? 


(AlB=1) % 
Zug Schwellbereich 
D n e N 
Wechselbereich Ermudungs 
festigkeil 
BEER Ey on: 
Ä rSprungstestigKel 2 Festigkeit 
Druck -Sehwellbereich ? bi: =] Ä Mm Bi rk Host, 
Wechselfestigkeit I] I 02= Bene Streckgrenze 
gÜ+= 185,4 
x Hin Ba 2} 
x ESEL 7 


Kain or en 


Pu 
+ 
ee Druck- 
A i ae 
Eee BETT 480 =.05 
2 (AlB=0) 
ee 
1 
ı 
m der Ermüdungsfestigkeit. Normalstahl 
| Ei wei Sana slett. allseitig bearbeitet. Dauerfestigkeit in 
> 
Abhängigkeit von der Verhältniszahl 
ı h A ou min 
\ _- Ermüdungsfestigkeit | | Deren 
r — f 4 
1440-2 (AB=) 
DApR. Be a 


Wahre a6 bzw. 


Längsspannung 0% 


Bild 13. Stahlguß. Grenzzustände 
von Hohlstäben. Beanspruchung 
Innendruck ein-, zwei- und dreiach 

eigenschaften des 8 


oO) 
KD) + 26,5/ 26,5 
25/25 


0) 
ee 
+ 28/14 A 


der Ursprungsfestigkeit 
durch Längskräfte und 
sig. Statische Festigkeits- 
tahlgusses: 


Fo er *2 
Et 5% DEE 
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j 
: | 
= 
| ; Br“ 2 J 
a = BR 
Bild 12. Hohlstäbe nach zweiachsiger Ursprungs-Beanspruchung. Oben: ‚‚Perunal P“, Ir = 3 Hang‘ 
Unten: Reinaluminium Ige = — Hang‘ 
28,5 kglmm® 
Bi: 
+ 


Verhältniszahl: 
5 

latter Stab = 4,6 
BE 


05 ® 
gelochter Stab =4,2 
405 287 kglmm? e' 


Bild 14. Ermügungsdiagramm von glatten und gelochten 
Flachstäben (d/b = 1/5). Flußstahl ‚‚St 60/70“ von 10 mm 
I Dicke. 


Streckgrenze o Zug- 


Ein- 


obere unt festigkeit | schnürung 
kg/mm?|kg/mm2| P, kg/mm® 


Hohlzylinder . 30,7 29,8 42,7 
Vollzylinder . . 29,1 27,0 


"u 
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$ = 10mm 


d=|72 


b: 80 


1 2 4 
Verhältnis der Abmessungen 


au Se b: 160mm 


4 
Verhaltnis der Abmessungen 


grosse Probestäbe 


160 |mm 


r= 40, 1,6 und Omm 
41 64 » oo 


r=1 ‚04und Omm 
A=t1;64 m 0 


kleine Probestäbe 


Bild 15. Der Erschöpfung durch Ermüdung unterzogene Versuchsstäbe mit flachen, halbkreisförmigen und 
tiefen spitzen Außenkerben. Gelochte Stäbe. Verhältnisse der Abmessungen. EMPA-Versuche 1945—1950. 


Bild 16. „Avional M“ b = 80mm, d = 72mm. 
Flachstab mit Bohrung. Ermüdungsanriß nach 
Druck-Ursprungsbeanspruchung. ; 


Bild 17 (rechts)! „‚Reinaluminium“ Flachstäbe mit Boh-- 
rung b = 50mm, d = 0,8 — 2,0 —6,4mm. Zug-UT- 
- sprungsbeanspruchung, Bruchflächen. 
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DEE SEHE EEE nn 


Ok e% = Om Br 


al 
ETF RP, SET Wen >>} 
+8 


| 
4 in 
r oO 
h Be 
100 = 
f 0, ! = E 
bzw. 100%) #15 = 
oO 
ee Salss 
r Era vr 
> = 4 
u: agalm 
2 0 .E 
ZEES e0+ =, > 
en 
ER: 


ım 


(glatter Stab) Formzahl Ock 


Bild 19. Flußstahl ‚‚St N“ und ‚‚St 60/70“. Zug-Ursprungsbeanspruchung gelochter und ‚gekerbter Stäbe. Blechdicke s = 10 mm. 
Ermüdungsfestigkeit. Zar Verstehe mit großen, allseitig bearbeiteten Stäben in niederfrequenten Pulsatoren Amslerscher Bauart. 


Bild 18 (links). Normaler Flußstahl, ‚StN“. Gelochter Stab, Breite 
b = 99,8 mm. Dicke s = 9,5 mm. Ermüdungsbruchfläche nach 
200 000 Lastwechseln zwischen 100 und 2600 kg. 


aus gemessener Gesamtdehnun | 
mit E= 24000 kg/mm® 


1-30 kgImm? 


Bild 23. Gelochter Stab der Bildes21. Eigenspannungen eo, und 
Spannungsgrenzen o, und 0, beider Ermüdungsbeanspruchung 
zwischen 1+24t — siehe 


23 000 
B — = ———= % 
ild 22 + 1250 18,5 kg/mm . 
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i SEN  5H60 


Zug -Ursprungsfestigkeit .04 — - o 
Zug-Druck-Wechselfestigkeil 0% — © _ © 


1,60 4165 4 158 
160 sfpsz]8 164 ee 
0153 447 Ok’ Ag 158 Ba] 19 9 + 
440 1,43 0147 + ' ©1447 
kr 138] Onag tl m KH Spannungsspitze Re du 
' j I : 
: ut N | k 
| 
(bzw. 100%) 44 | ' | 
100 | | 
| 
N 08 | 
03 0,88 | | l 
ER 087 . 
= SORT 087 0,73 N ® | 
= Z 83 = 0,65 J 
& & ‚| 07511 0,75 2 } 0,66 
| Ka 3 He Ye 
£ | | E | PRnSpannung 059 Tu 066 
22 |. 18|els 2 8 | :|e 21 
8 In =D ıPD|:S|Em > m — Ex BES zB ! 
3) 2) 8 ı 215 |5 er ı 88 | 15 2 | 
S 8 S 15|8|15 & 53 el ıle> 5 | 
EZ ul 1.818 — | —_ 5— ® 
cc vr 22|2|88 u 25 > BRCI= j 
nee ee : bee ser 2 
Er ms Ss Me s| ıl®S S l 
ie Rallis u °| 5||= © WISS & IISS & | 
© 1 = ı=&|= = = —< Kg = 
= Ss j = I ER I EB j 
4 a 1 "- 156 161162 186 5% 2,90 3 
ı {) 
‘glatter Stab) Formzahl OXx 


Bild 20. Flußstähle ‚‚StN“ und ‚‚3t 60/70. Ermüdungsfestigkeiten „o gelochter und gekerbter Stäbe und die dabei erreichten 
rechnerischen Spannungsspitzen &p' g0- Niederfrequente Pulsatoren Amslerscher Bauart. 


b =150 mm 
(5 = 9,Bmm) 


Messlänge : 


e0; 16205 


-28 kglmm* ; 


Bild 22, Gelochter Stab des Bildes 21. Ermüdungsdiagramm und die 
Spannungsgrenzen an den Meßpunkten (4) ; (2) (6) und () bei der 
Ermüdungsbeanspruchung zwischen P=1 + 24t. 


ap elastische Formanderung 
15 26 2,47 3 2,38, in oo D 


bleibende. Verformung 


4-26 496 492 nl 
0,85 IN #00 » 
050 
1=24t 400 0,1 
004 0,06 
| 1-46 04 008 
1:18: 008 e 


1 x 0,08 
20mm, mit Hilfe elektrischer 


ungs r (Strain gages). Erster Ermüdungsversuch: P=1 + 21t, Versuchsstab nach 10° Lastspielen intakt — Last erhöht. 
es : Se ans ee dungerereten: P=1-26t, Anriß an der Stelle (4) nach 100 000 Lastspielen. 


x 
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SEN 5850 Cr-Ni-Mo-Sh 5%Ni-SH. 


Rundstäbe — «® o d u ; "194 
Flachstäbe — = 6 Pr es 1874 114 Flachstab Rundstab 


07 72a [0,79] 
0,61 ae! 
u“ 047 


I 
| 1 
! l 
: l 
n j 
| ı 
| | 
! j 
! j 


G 
25 ee. 7 


188 2,05 2,13 


x 2,85 
/ 4 TE - 


latter Stab [| o x = 
i #Zug-Drunk glatter Stab Formzahl org = mas „3. en _ _ 8%, für lachsäbe, bu.47-EÜ_ = 470C%-Für 
| *Biegung eOn 2 MjWnetto 2 M/Wnetto Rundstäbe. 


N 


Bild 24. ‚‚StN‘‘, ‚„St50°, „5% Ni-St‘‘ und „Or-Ni-Mo-St‘‘. Dauerbiegeversuche mitTaußen gekerbten Rund- und Flachstäben. 
Anzahl der Lastwechsel beim Bruch 10 -10®, 
M 
W, 


Rechnerische Biegewechselfestigkeit ee + = ne, zur (Flachstäbe) bzw. + = 
netto 


ax "ken 


=1,7 PR (Rundstäbe) 


netto 
und die dabei erreichten rechnerischen Spannungsspitzen PL 


r 32,5 kg Imm® 


Nennspannun \ A 
0m? 23 kam 


+ - 116 %o 


gemessene Rückfederung 
nach dem Herausschneiden 
der Messpunkte. 


37,5 mm 


-24,0 kglmm? 


Bild 25. Normaler Flußstahl ‚‚StNes. Gelochter Flachstab: b=150mm, d= 75mm. Eigenspannungen o,; und Spannungsgrenzen 
©, und o,, nach 600 000 Lastspielen zwischen P=1 16,5t. Anordnung der elektrischen Dehnungimeasen, 
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Bei gleichzeitiger Wirkung einer Längskraft P (Zug), eines Biegemomentes M, und eines 
Drehmomentes M, (z.B. kreisrunde Welle, D=2r mit umlaufender Kerbe) in einem bestimmten 
Lastintervall A/B, kann die resultierende Formzahl ax ,.. gemäß der Gestaltänderungshypothese 


(Os 0,— Vo?-+-372) angenähert wie folgt definiert werden 


Er EEE 
Vorne] 


3. Eigenspannungen 

Eine Folge des Spannungsabbaues sind Eigenspannungen im entlasteten Versuchs- 
stab (Bild 25), die jedoch in diesem Fall mit Verfestigung, Zerrüttung und Alterung verbunden 
sind. Ist der Körper schon vor dem Gebrauch mit ungünstig verteilten Eigenspannungen behaftet 
gewesen, ohne daß damit eine ‚Schädigung des Materials verursacht wurde, so werden sie während 
der wiederholten Beanspruchung derart verändert, daß. für ein bestimmtes A/B-Verhältnis 
im Ermüdungszustand praktisch stets dieselben Eigenspannungen vorhanden 
. sind. Der normale Flußstahl wie auch die einwandfreie Schweißung verfügen über ausreichende 
dynamische Zähigkeit — beispielsweise durch die Dehnung an der Stelle der Spannungs- 
spitze ausgedrückt —, um die ur- 
sprünglich vorhandenen Eigenspan- 55 ‘ 
nungen umzubauen, a Schon Kahn 
wenige Promille bleibender Deh- N sr. o|70, 820 32,2 34% 
nung genügen. 


AK res = 


Druck- Ursprungsfestigkeit 
Ar 


[e) 


[EN 
[ >} 


4. Einfluß des Spannungs- 
gefälles - 

Ein Einfluß des Span- 
nungsgefälles konnte bei tech- 
nisch normalen Abmessungen der 
Querschnittsübergänge, in einem 
weiten Größenbereich (1:8) nicht 
festgestellt werden (Bild 26—28). 
Opnaz 
di 
ist nur dann von Einfluß auf die 
Größe der Ermüdungstestigkeit, -. 
wenn die Dicke der von der Be- 
arbeitung herrührenden, kaltver- 
festigten Randschicht, im Verhält- 
nis zur Kerbgröße bzw. den Ab- 
messungen des maßgebenden Quer- 
schnittes ins Gewicht fällt (Bild 27). 


5. Einfluß der Körpergröße \ 
Ein Größeneinfluß (Bild 5 
:28) fängt sich erst bei äußerst 
scharfen Kerben an bemerkbar zu 
machen, und zwar zeigt mittelharter 
- Flußstahl im Falle großer Versuchs- h) 
körper mit tief und so scharf wie 12222 
möglich ausgeführten Außenkerben 


ei ug- rungsfestigkeit, die Bild 26. „,St N“ und., ‚St 60/70‘. Einfluß der Gradienten der Spannungsspitze 
en Ze at 8 8 x auf die Ermüdungsfestigkeiten ‚oo, „OT und ‚ot. Flachstäbe. Homo- 


nur etwa die Hälfte re SE loge Reihe d/b = 1/5. Verhältnis der Abmessungen 1 + 8. Versuche in: nieder- 
reicht, welche an kleinen homologen . frequentem Pulsator Amslerscher Bauart. 


Stäben gewonnen wurde (Bild 29). 


6. Einfluß der Verschärfung der Kerbe 
"_ Bei sehr kleinen Kerbradien wurde eine Grenze gefunden, von der an kein Einfluß der 


. 


weiteren Kerbverschärfung festzustellen ist, so bei Flußstahl etwa o—=0,5 mm (Bild 29) und bei 
_ Aluminiumlegierungen ein Viertel davon, etwa e—=0,1 bis 0,15 mm (Bild 30), vorausgesetzt, daß 


nn ——-. 


(Tafel 1) 
| 23,5 +4% 
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Das Spannungsgefälle 
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Verhältnis der Abmessungen 
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keine Anrisse in der Kerbe schon vor der Bearbeitung vorhanden sind. Es scheint demnach, 
daß dem Ermüdungsanriß die Zerrüttung einer Zone um die Stelle der Spannungsspitze herum 
vorausgeht, so daß sich auf diese Weise eine wirksame, wenn auch nicht direkt wahrnehmbare 
„Ausrundung“ von selbst ausbildet, die dann von einem Einfluß sein kann, wenn die vorhandene 
Kerbe schärfer als die ‚„Ausrundung“ der zerrütteten Zone ist. 


III. Zusammenfassung 


Die Tatsache der Ermüdung von Metallen, die sich nur als Folge wiederholter Bean- 
spruchungen bestimmten Größenintervalls einstellt, kann nicht in Abrede gestellt werden. Sie 
wird durch die Erfahrung und das Experiment bestätigt. Der Begriff ist der organischen Welt, 
unserem eigenen Leben entnommen. 


Das Problem der Ermüdung gehört zu .den verwickeltsten Problemen der Werkstoff 
mechanik. Zu beachten sind: 


die elastischen, plastischen und überelastischen Verformungen, 

die Hysteresis, 

die inneren Spannungen, Eigenspannungen erster, zweiter und dritter Art, 

die Fähigkert des Spannungs- 
2 Ac. N” „1947° b= 50mm Spitzen-Abbaues ausnahms- 
'» 1949” b= 12mm (--) b=40mm weise deren Stützung, 

die Dämpfung, 
die Erholung, 


die Verfestigung und 
Entfestigung, 


die Alterung und 
die Zerrüttung. 


» 


Zu all dem kommt noch der 
in der Natur der Entstehung, Er- 
zeugung und Formgebung der 
Werkstoffe liegende Umstand, daß” 


PT Homogenität 
r Stab: 1909,29,7 | Isotropie 
Ouy ar | Ideal-Begriffe sind, welche die 
Wirklichkeit nicht kennt. 


Der Ermüdungsbruch ist 
weder ein Gleitungs- noch ein 
Trennungs- und auch nicht ein 
Verschiebungsbruch. 


4-0 Mo1 ho2 057504 05 06 %07 W009 4 
L=0 08 087 0351475 262070 35 4381,05 140 (50) 50 100 25120) 6,0) Omm 
Bild 27. Normaler Flußstahl ‚‚StN“. Grenzwerte der rechnerischen Ermüdungs- 


‘ Die Erschöpfung durch 


 festigkeiten @,y = 505 Und Gy man = efman =%r "40x. Einfluß des Ga- Ermüdung stellt sich, bei gesun- 


dienten bzw. des Verhältnisses d/b. Ergebnisse der EMPA-Versuche 1947 und 1949. den, praktisch fehlerfreien Werk- 

stoffen als Folge allmählich, mit 
der Anzahl und Höhe der Spannungswechsel fortschreitender Störungen zufolge über- 
elastischer und plastischer Verzerrungen und Verbiegungen des Raumgitters, welche eine 
Schwächung der Kohäsion und schließlich Zerrüttung und örtliche Trennungen — Phase 1: 
erster Anriß, Phase 2: Ausbreiten des Anrisses ohne wesentliche Verformung — zur Folge 
haben, ein. Der Widerstand des dadurch geschwächten Querschnittes wird dann durch 
Gewaltbruch überwunden. Die durch Ermüdung erschöpfte Bruchfläche weist überwiegend 
intra-kristallinen und nur ausnahmsweise inter-kristallinen Verlauf auf. 


‚ „Die Erschöpfung durch Ermüdung ist grundsätzlich eine Wiederholung sta- 
tischer Beanspruchungen, welche sich aber auf einen in seinem kräftemechanischen 
Verhalten, durch die wiederholte Inanspruchnahme, jeweils veränderten Werkstoff auswir- 
ken. Hier liegt die Schwierigkeit der theoretischen Fassung des Ermüdungs- 
problems, indem sich der jeweilige wirkliche Spannungs- bzw. Anstrengungszustand unserer 


genauen Kenntnis entzieht. Die überelastische Verformung.scheint anderen Gesetzen als die 
plastische zu folgen. | 


* 


Pe 
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Unter diesen Umständen erscheint es ge- Modell Et _ GO g max _ g5 kalmm 
rechtfertigt für die Erschöpfung durch Ermüdung 4 Be > da ne 
in der Spannungsamplitude der stationä- % 9=$- 0,257 
ren Ermüdungsphase einen ‚elastischen u 
Maßstab“ zu wählen. den 
Bestimmend für die statische Anstren- . &=2,42 


gung und Ermüdung sind: 


1. bei Körpern ohne Kerbe mit homo- Ausführung 
gener Spannungsverteilung: 10 x 


a)im Falle einmaliger Belastung bis 
zur Fließgrenze die resultierende Schub- 
spannung 7,, der Oktaederebene des Haupt- 
spannungs-Würfelelementes bzw. die Vergleichs- 
spannung 0,;—0,, somit die von der EMPA er- 
weiterte Theorie von O. Mohr (Bild1 u. 5); %=07 

b) im Falle einmaliger Laststeigerung 
bis zum Bruch die maximale, wirkliche 
Schubspannung 7,.. für den Gleitungsbruch 
und die maximale wirkliche Hauptzugspan- 
nung Oyas für den Trennungsbruch, in 
Übereinstimmung mit der Theorie von ©.Mohr 
(Bild 1u.5) und schließlich 

c)im Falle wiederholter Beanspru- 
chung auf Ermüdung bis zum Erreichen der 
Ermüdungsfestigkeit (Dauerfestigkeit) selbst 


Om: 1,0 


u 4 
[i 


die Hypothese von Coulomb-Mohr zu- Bild 28. Gelochter Stab. d/b = 0,257. Homologe 
. .. . ’ . Abmessungen 1:10. Einfluß des Gradienten 

mindest für die Wechsel- und Ursprungsfestig- da 

keit (Bild 1 u. 13). 270 


Dazwischen muß man für den jeweiligen vorliegenden Fall die maßgebende Anstrengung 
schätzen, um so die Sicherheit eines Werkstückes beurteilen zu können. 


Breitseiten der Probestäbe 


im Schnitt 
bei SF,N s=1«44Imm 
n 5560 s=K«3>1 mm 


x) his 10:40° gleich viel 


| 
| 
| 
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D 
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Ermüdungsfestigkeit = Nennspannung e Echt in kg /mm® 


0 
A -In5 m 3 4 6,4 9 27 © 
- . e 0 
djb =o21 04 02. M 0,04 0,02 0,01 Ela 
möglich 


3 { m. -Ermüdungsversuche g+ mit außen ge- 
Bild 29. Stähle ,‚StN“ und ,‚,‚8t60/70°, Blechstärke 20 mm Zug-Erm 0% 
* kKerbten Flachstäben =1+m, Amsler-Pulsatoren hoch- und niederfrequent. Anzahl der Lastspiele 10°, 
So Einfluß sehr scharfer, tiefer Kerben auf kleine und große Versuchskörper. 
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Bild 30. Reinaluminium und Leichtmetall-Legierungen ‚‚Avional M‘‘ und ‚.Avional PT‘‘. Blechstärke 15 mm. 
Zug-Ermüdungsversuche mit außen gekerbten Flachstäben a, =1-+ »: e05 und ‚o. Hochfrequenter 


Amslerscher Pulsator. Anzahl der Lastspiele 10°. Einfluß der Kerbenverschärfung. 


Rund - Stab: gekerbt ungekerbt 5 
auf Zug beansprucht : : 


Ermüdungsfestigkeiten 
I 0. O) 


RE Effektive Zug-Ermüdungs- 
1 Festigkeit des Materials 
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Nennspannung 


Rechnerische auf den 
netto Querschnitt bezogene 
Ermüdungsfestigkeit 
des gekerbten Stabes. 
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Bild 82. Rund- und Flachstäbe gekerbt. Kerbe tief vo 
hyperbolischer Form. &,-Größe nach Neuber für Fr N 


a Y - 2 
N. 3 verschiedene Werte von = - 
Inar E . 4 f je b 
Bild, 31. Gekerbte Körper. Ermüdungsbeanspruchung. Grundsätze der k . \ > 


X : EMPA-Richtlinien. 
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_ 


2. bei Körpern mit Kerbe, mit heterogener Spannungsverteilung: im Falle 
der Ermüdungsbeanspruchung die Größe der maximalen Spannungsspitze 0,4. bZW. 
der maximalen Anstrengung o,,—=0, (Bild 31) und die Verteilung der Anstrengung 
über dem Querschnitt der Ermüdungs-Bruchfläche im Sinne der &,-Zahl als 


__ _Belastungsgröße in vollplastischem Zustand 
Belastungsgröße bei rein elastischem Verhalten 


Eigenspannungen, welche sich als Folge des Spannungsabbaues im entlasteten Ver- 
suchskörper einstellen, sind im Ermüdungszustand für ein bestimmtes A/B-Verhältnis praktisch 
unabhängig von den Eigenschaften, die vor dem Versuch vorhanden waren. 


+ IC’ man do ,es 
| en bzw. en] 
Einfluß, wenn die Dicke der von der Bearbeitung herrührenden, kaltverfestigten Randschicht, 
im Verhältnis zur Kerbgröße bzw. zu den Abmessungen des maßgebenden Querschnittes ins 
Gewicht fällt oder wenn die Randschicht durch eine Sonderbehandlung der Oberfläche ihren 
materialtechnischen Habitus grundlegend verändert. \ 
Große Körper mit tiefen und sehr scharfen Kerben weisen eine ganz wesent- 
lich tiefere Ursprungsiestigkeit o, auf als kleine, homolog ausgebildete Versuchskörper, weil 


0 
. sich die infolge Zerrüttung im Kerbgrund ausbildende ‚Ausrundung‘“ als scharfe Kerbe 
am großen Körper ungünstiger auswirkt als an einem kleinen Körper. 

Eine weitere Verschärfung sehr kleiner Kerbradien: bei Flußstahl von o=0,5mm 
und bei Aluminiumlegierungen von o=0,1 bis 0,15 mm übt beinormalen Werkstoffverhältnissen, 
keinen weiteren Einfluß auf die Größe der Ermüdungsfestigkeit (Ursprungs- und Wechsel- 
festigkeit) aus. _ r 

Die Hypothese von Coulomb-Mohr behält auch im Falle gekerbter Körper 
mit heterogener Spannungsverteilung ihre Gültigkeit. 

Dreiachsige Grenzzustände für statischen und Ermüdungsbruch werden 
‚auf zweiachsige zurückgeführt. In der weitaus überwiegenden Anzahl von Fällen ist der zwei-'_ 

- achsige Spannungszustand maßgebend, wodurch sich eine sehr beachtenswerte Vereinfachung 
ergibt. - | ; | 
; Ist die Kerbziffer a, bekannt (Bild 32) dann müssen für eine bestimmte- Werkstoff- 
ttung, einerlei in welchen maßgebenden Intervallen sich die Grenzkräfte — Spannungen bzw. 
Anstrengungen — bewegen, die entsprechend der Elastizitätstheorie berechneten o,4.- DZW.. 
G,.,-Werte=0;*,0, bezogen auf den glatten, nicht gekerbten Stab, aus entsprechenden sinn- 
gemäßen Laboratoriums-Versuchen, bekannt sein (Bild 19) wodurch sich dann für den jewei- 
ligen Fall der zugeordnete, maßgebende ‚„-Wert ermitteln läßt (Bild}19), welcher 
durch den Konstrukteur der Bemessung des Werkstückes selbst zugrunde zu legen ist. 
Das Phänomen der Erschöpfung an der Fließgrenze, der Erschöpfung 


[0,77 


— (Bild 24). 


Das Spannungsgefälle — der Differentialquotient ist nur dann von 


durch statischen Bruch als auch der Erschöpfung durch Ermüdung wird durch die Fr 


Theorie von Coulomb-Mohr von einem einheitlichen Gedanken erfaßt. 
Eingegangen am 23. April 1951. 


Die in Achsenrichtung rotationssymmetrisch belastete dicke 


Kreisplatte auf nachgiebiger und auf starrer Unterlage 


- Von Istvan Szabö in Berlin 

Im Anschluß an frühere Untersuchungen wird das Problem der in Achsenrichtung rotationssymmetrisch 
belasteten dicken Kreisplatie auf die Fälle nachgiebiger bzw. starrer Lagerung ausgedehnt. Für die an der 
EN Berührungsseite auftretenden Schubkfäfte wird das CO oulombsche Reibungsgesetz angenommen. Für 
die eingespannte Platte gelingt es, exakte Lösungen anzugeben, während für die am M. antel kräftefreien ‚Platie 
(freie Lagerung) Näherungslösungen im Sinne des de Saint- Venanischen Prinzips gegeben werden. 


ar: In continuation of previous investigations, the problem of the thick circular plate under rotationally- 
ar  symmetrical sirain acting in awial direction is extended to the case that the plate is bedded elastically or 
rigidiy. The Coulomb friction law is assumed for the shearing forces acting at the surfaces in contact. In 
the case of clamped borders, it is possible to give ewact solutions, while, for free borders, approximative 
> solutions are given, by use of Saint-Venant’s principle. 
En continuation d’investigations anterieures le probl&me de la plaque circulaire Epaisse, chargde dans la 
direction de l’axe symmötriquement & la rotalion, est iendu sur les deplacement ölastique et rigide. Quant 
aus forces poussant du cöl& de contact, on suppose qu’elles agissent selon la loi de friction de Coloumb. Pour 
la plague pinote il est possible d’ indiquer des solutions ewactes, tandis que pour les bords de la plaque 
depourvus de forces (placement libre) des solulions approwimatives sont donnees selon le principe de Saint- 
25 Venanl.— =. EEE HL nr | 
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B npono:mzeHne IpensIayImuX HCcTeNOBAHHH IpoÖ.leM3& TOAICTOÜ KpyToBOH INIACTHHEI, BPAIIa- 
TeIBHO-CHMMETPHYHO HATPY>KeHHOH B OCEBOM HAIPABJIeHHH, PACHPOCTPaHNeTcH HA cayaaHh 
NONNAMINEHCA, HIH IKE IKECTKOH ONOPEI. st CKAJIEIBAIINUX YCHAHH, BA3HUKAWIIHX H& CTOPOH® 
COIMPHKOCHOBEHHA, IIPHHAT 3aKOH TpeHun KyoHa, ToyHsie pemeHua MOrTyT ÖbITb HaÄNeHBI 
3a3KaTOH INTACTHHBI; HA CIIyyall >Ke INIACTHHBI C CBOÖONHOK OT CHI 0OKOBOH IIOBEPXHOCTEI, 


AAIOTCA NPHÖ.HNMKEHHLIE perueHus B CMblcıe upunmuna ÖeH-BeHana. 


1. Einleitung 

In zwei früheren Arbeiten habe ich das Problem der rotationssymmetrisch belasteten 
dicken Kreisplatte auf elastischer Unterlage ohne Berücksichtigung des Eigengewichtes behan- 
delt?). Diesen Untersuchungen lagen die Hookeschen Gesetze zugrunde, während A. Meffert 
dasselbe Problem für eine nachgiebige und vollkommen glatte Unterlage gelöst hatte?). Unter 
„nachgiebig‘“ wird hier im Anschluß an K. Marguerre?) eine solche Unterlage verstanden, 
bei der an Stelle des Hookeschen Gesetzes die Annahme getroffen wird, daß die Verschie- 
bung dem örtlichen Normaldruck proportional sei. 

Das Ziel der folgenden Ausführungen ist es, das Problem der nachgiebig bzw. starr gelager- 
Be: ten Platte unter Berücksichtigung einer rauhen Unterlage einem ähnlichen Abschluß wie bei 
er der elastischen. Unterlage!) zuzuführen. Die Berücksichtigung der Oberflächenrauhigkeit soll 
in der Weise geschehen, daß die Schubspannung an der Plattenunterseite im Sinne des Coulomb- 
schen Gesetzes dem übertragenen Normaldruck. proportional gesetzt wird, wobei der Propor- 
tionalitätsfaktor die Haftreibungszahl ı, ist. 

Einige Resultate einer, dieser bald folgenden Arbeit vorwegnehmend, wird zum Schluß 
noch die exakte Berücksichtigung des Eigengewichtes angedeutet. 


Be; 2. Problemstellung 

Eine Kreisplatte von dem Halbmesser a, der Höhe h, dem Schubmodul @ und der Poisson- 
schen Zahl » ruhe auf einer nachgiebigen Unterlage von der Bettungszahl X und werde durch 
den achsensymmetrischen Druck f(r) belastet (Bild 1). Die Randbedingungen lauten wie folgt: 


Ta 0er er Br Rd ER REE (1), 
ONE. 2 ee ee ee (2), 
0 N) = Re 
o,(h, N) = KR Chr) ae Ta ee Gy 


Bild 1. Die achsensymmetrisch belastete dicke Kreis latte auf 5 
a (bzw. starrer) Unterlage pP nachgiebiger 


(Für eine starre Unterlage ist diese Randbedingung durch die Forderung 
zu ersetzen.) NS NS ee ee . (4a) 


UNE) ee ae ee (5), 


= (0) ADYN AR N er 
) Ing.-Arch. XIX (1951), 8.128 und | = 
2) A. Meftert: Diss. T. h ee Be Re 
3) Ing.-Arch, IV (1933), 8.341. 
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Die Forderungen (4) und (5) bedeuten, daß zwischen Platte und Unterlage stetiger Kontakt 
besteht, was in Wirklichkeit solange der Fall ist, wie o,(h, r)< 0 (Druckspannung) ist; wie man 
vorzugehen hat, wenn Zugspannungen auftreten, wird uns später noch beschäftigen. Die Rand- 
bedingung (6) bezieht sich auf die in Bild 1 dargestellte freie gelagerte Platte; für eine am Rande 
(r=.a) fest eingespannte Platte ist sie durch 


NEN RE RE ee (6a) 
zu ersetzen. 


Die Spannungen o,, o,, o; und r (Bild 1) lassen sich durch die Vertikal- bzw. Radialver- 
schiebung { bzw. o bekanntlich folgendermaßen ausdrücken: 


| 
=20(& Re 3)’ e 
G eo = > a[& x) Eu ne ma ie 9 
n=20| 4.) Kr er - 
wobei 
el 8 
a BR (8) 
bedeutet. 
Aus (7), (8) und den Se 
00, , 0, — 0% & 9 
Fr r +; rt, 3 9) 
ergeben sich für die Verschiebungen die RN 
0) 0 
een 1 a (10), 
wobei 4=- tz ER) Br ar Are Laplaceschen Operator bedeutet. 
r r 


3. Lösung für die am Rande eingespannte Platte 
Nach meiner ersten Mitteilung) sind die der Forderung (1) genügenden Lösungen 


DE [Bee + Quer? — or (pe + a) Jar) +Ae+B - - A), 


S fı (3—4) 21]. —ı - an Be 
0e= 2 rtaaan)a, Zr ann ar 
ir 2e 452 N J y\ r 
2t— an Pt re ) ı( k ) 
h 

Hierbei bedeuten Jy(A,r) und Jı(Arr) die Besselschen Funktionen erster Art, A, den nac 
Tu a)—=0 es an P,, Q&; Pr, (x, A und B noch freie Konstanten. (Bezeichnen 
wir er die unendlich vielen und positiven Nullstellen von J,(t) mit t,, so sind die A, durch 


= = gegeben.) 
Aus (2) fließt mit (7), (11) und (12) 


| (Pet At 2 m OL I ge ed (13), 
während die Randbedingung (3) auf ähnlichem Wege zu 
20 3(-%-3- Bit rn): In + Az en = — fl) 


führt. Multiplizieren wir diese Gleichung mit rdr bzw. rJo(Anr) B u 1,2,...) und integrieren 
von r—=0 bis r=a, so ergibt sich wegen 
% % : 0 für ken 


frane=©: [aarana= Io) fir ken 
ar 0 
1— 2» 


er lt > | EN S yarE (14) 


3 a) A.Nädai: lastische: Platten, 8. 310. 


3) Ing.-Arch. XIX aus), 8.128. Er 
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% 


bzw. 


— [1 In Qyr) rdr 
-ZmRtmW+-AtmWa=-— an  & 0. (B5). 


Setzt man die aus (13) und (15) berechneten P, und Q, in (11) und (12), so erhält man: 


d,2 A.2 er En 2] 
= Sn El Mur? k Ho (349 2322) e'**] 


k=1 


+2 Sin As! Il) +42 B 2.2.06), 


= 22 2 ÄLz Az — 4.2 
Sn 1,0 [8 Ay—2%,2)e "% Eee a3 [8 49+24,2)ek’+e "%*] 
— 7,6] near) ERERL N (17). 
Mit diesen Verschiebungen liefert die Randbedingung (4) 5 
er 77 yh 6; 
20,- Ze 5 3) he AR _ SinA;h) — 5 3) (ühet Sin A;h) 
nr en “ a) =—K = en [84,422 He] | (18). 


+ 4 Ay ger a E= (3—4 vw—2 A; h)e'*"] + „Sin Ahllr)- — K(Ak+ B) 


Nach Multiplikation mit rdr bzw. rIo (An r)dr und Integration zwischen r—=0 und r=a er- 
halten wir 


EN 20») 
Palit anne a 2 Re (19) 
bzw. 
pre" "— Einah gl Ay+2%,h)e ee) 
+4- Anhe'? + Sin Art z;, Br 21) ee ) (20). 
= 21 29)0u(Eojah + an 2 Ein An) 


Aus der wre (5) ergibt sich mit (16) und (17): 


er SinA,h)+ he Sin 4,h)—2 o, Sin u J(Arr) 


1 = > BE 
BE 52 + I na re HR Sin Ach). — er „het (21). 
— Sin Ah) + 0 Cof A, ı Tl) | 


/ - 


Multiplizieren wir diese Gleichung mit rJy(A„ r)dr bzw. rdr und integrieren zwischen r=0 und 
r=4, so bekommen wir: 


> [Pe (Arhe "+ Sin A,h)+ qu(Ay he’*" + SinA;h) Sal Sin EACH) Jı(d2r) dr 


— [Du (Anhe Sind h)—qu(duh en" — Sind, DEE ORTE) te) (2), 


bzw. 
Zeh "+ Sin un) + qu(Auheir”+ ©in A, h) | i ze we H 


—20(1— 29) &inyh] rFylgr) dr = mal —»)A (a). 
ö | SE eT 
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Setzen wir.nun aus (20) g; in (22) bzw. (22a) ein, so bekommen wir off i : 
’ enbar ein dliel 
Suchwesyptem für die Koeffizienten p, in folgender Gestalt: a 


3 [P(a)ps + O(A)oı Fo [R(A,)pn + 8(,)o,] Fi 2,0) 


(23). 
en Pr +9 (Ar) or] ER = ma? (1 —») A 
Hierbei bedeuten: 
= wir K- IT a) 
PiA)=— 2 (Sin?A,;h — A; h?) + (24,;,h + Sin2 ” n) RER (24a), 
Q(A) = 2(1—2 ») Arne SinA,h = — no er" Sin Ah] WERE RT (24b), 
F(A)= he" — Siny4h— u [8—4v—2uhlet" Le" ...... (24c), 


S(A,)=2(1—2») S ,h+e "Sind, — = SinA„hl(—v) en" Sina, 2) (25a), 


R(4,)= 


K 3 2.72 a 
amt ein; EN, fr gyaun)dldr)dr (250), 
N 0 


b; = IrIı (A; r) dr A Ne We: (25€): 
0 


Das Integral für a;, läßt sich zwar mit Hilfe der bekannten Formel für Produkte von 
Besselfunktionen®) auf eine unendliche Reihe zurückführen, die aber schlecht konvergiert, 
so daß man das Integral bequemer mit der Simpsonschen Regel bestimmt. Auch b, wird 
man auf diese Weise bestimmen. 

Da die durch (17) gegebene Verschiebung o=o(2, r) wegen J,(A,a)—=0 der Bedingung (6a) 
genügt, ist damit das Problem der am Rande eingespannten Platte gelöst. Der u. U. notwen- 
digen Forderung der vertikalen Unverschiebbarkeit der Stelle z2=h, r=a, kann man in der 
Weise gerecht werden, daß man in (16) &=£(2,r) durch &(2,r) —£(h, a) ersetzt. 

Über die näherungsweise Auflösung des unendlichen Gleichungssystems (23) und der 
weiteren ähnlichen — in dieser Mitteilung noch auftretenden — Gleichungen läßt sich dasselbe 
wie in meiner zweiten Mitteilung”) aussagen: Die dort angeführten Gedanken lassen sich fast 
- Wort für Wort auf die hier behandelten Probleme übertragen, so daß ihre Wiederholung über- 
flüssig erscheint. 


4. Lösung für die freigelagerte Platte 

Für die Erfüllung der für diesen Fall geforderten Randbedingung (6) bleibt nach Auflösung 
des Gleichungssystems (23) keine Möglichkeit mehr, da sämtliche Konstanten durch die vor- 
angehenden Randbedingungen (1) bis (5) festgelegt sind. 

Der noch nicht erfüllten Randbedingung (6) kann man im Sinne des de Saint-Venant- 


schen Prinzips folgendermaßen genügen: 
Für r=a erhalten wir aus (16) und (17) eine mittlere Ba 


+ eenenadg Er = (he %#* 4 Sin 4;h) 


(26) 
A vr) 4: a 
a 3025) he »® I Ein A,h) — co; Sin A;h .- a, 
und ein el ouen 
nn [Inc TE rel Ne Anz, 3d— 25 u an re ER 1 24h Cord 
—2 Sinn) te : „met — 26h Ban hen] 


JA 
+ 2 Ein };h) — 9. [Ach —ı A;h+2(1 —&oj 2) ee, 


6) Siehe z.B. Watson: N of Bessel Functions, S. 148. 
?, Ing.-Arch. XIX Bo 


. ut R & EN Pe 9 x 
. 4 x : £ u EEE ar 
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Nun können wir der Forderung der Kräftefreiheit des Plattenmantels bezüglich der Resultieren- 
den so gerecht werden, daß wir den Plattenmantel durch —o, auf Druck bzw. Zug, und durch 
das Kräftepaar — mr auf reine Biegung belasten. Die zu diesen Belastungen gehörigen Ver- 


schiebungen sind bekannt®): Zu —0,: 


— v0, (4 ); A 
eier, Ian Et 
zu — Mr: 
R Be en Fr N AA 7 
3mn|2»(2—) Han male —z)r 
a A-+)@R EEE Ut)am ö 
x Die entsprechenden Spannungen ergeben sich nach (7) zu: > 
5 ,=0; %=0I; =h——0l, . 3 
& 12m Be 29). 
2 ? 0: n=03 og, =: = — | ae ® ) 
% > Aus den letzten Gleichungen ersehen wir, daß von unseren aus den Randbedingungen fließenden 


Gleichungen lediglich die aus (18) hervorgehenden Gleichungen (19) und (20) entsprechend der 
zusätzlichen Belastung zu modifizieren sind. 


An Stelle von (18) tritt jetzt die Gleichung 


477, 5 Pe 3 —A; h Ik or: 
2 (A, he kon a aa k — Ein};h) 
E I ech. & 
+ aa an Sa rt ne 4 
Ik —ı.h HR re (30). 
a 4%,(1—2») [e k ae k I+7, Sin Ah Jo (Ayr) 
o,h 3 mp [vh?+2(1— ») r2] 
RA DIR tet eRer ee 
eg 2(1+»)@ 2(1-F»)@R® 
Multiplizieren wir diese Gleichung mit rJ,(A,r)dr und integrieren zwischen Null und a, so be- 
kommen wir wegen r » 
f 2.a2 | 
Sr Ir) dr = — Io (ma) 
N) 42 
en 9 En St 3? Unzer Ah : 
20525, he — Sin Ah) — pt, Anbei? — Gin Anh) 
m — et N Ze —A,h Anh A 
+0, Go an | A Art 2 net 4 ah] Rn 


I TA De Ein ER 
44,(1—2») „h) +7, n (+ »)GR372 T,(A,a) 
Denken wir hier aus (27) mr eingesetzt, so stellt (31) für n=1,2,3... ein unendliches, li 
; Ran. > j ae ae lineares 
Gleichungssystem für die Koeffizienten p;,und q, dar: Sie kö Be 
und (31) berechnet werden. 2 ee ie können aus den Systemen (22), (222) 


N? Multiplizieren wir (30) mit rdr, so liefert jetzt die Integration: 


Gar 0 of RS or ap ea 
ee Fe RER v®+(1 2 
Te ee. Rode Pe .S=—EBS, 


- 


“ 


womit nach Auflösung der Gleichungssysteme (22), (22a d j in © 
mr die Konstante B gefunden ist: ur ne = re Bis 


B-_ a 208 1, __ eh _3mbR+l—ne 3 
Ar ee rin 


Damit ist auch dieses Problem gelöst. 


8) Ing.-Arch. XIX (1951), S. 136. 
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5. Die Kreisplatte auf starrer Unterlage 
Die Randbedingung (4a) liefert aus (16): 


oo 


En ck M en 
Dan le Ar+ 2). EN Er en 


Er 33). 
(9 —4y— 2A, h)e'#"] + Sin ih Jo(A,r) + Ah+ B=0 = 
Mit rdr bzw. rJ,(A,r)dr multipliziert und integriert, erhalten wir 
RA Pe TEE (34) 
bzw. 
er ah A, Ir ee - 
Se Ey LIA,he ee ne 
IRA N FEN I gnd—2s).° | (35) 


+ (3—49— 24, h)e*kr] = Siny,h=0 | 
1 


(Selbstverständlich hätte uns eine andere Überlegung bezüglich des Bestehens von (33) zu dem 
gleichen Resultat geführt.) 

‚ Denken wir aus (35) g, in (22) u. (22a) eingeführt, so erhalten wir ein lineares unendliches 
een für die 9, womit —da (6a) erfüllt ist —für die eingespannte Platte das Problem 
gelöst ist. = 

Für die frei gelagerte Platte ist auch in diesem Falle die Lösung nur im Sinne des de Saint- 
Venantschen Prinzips herzuleiten: Wir belasten die Platte am Rande wieder durch —o, und 
— MR (wobei 0, und mz.durch (26) und (27)gegeben sind) und haben aus der Randbedingung (4a) 
mit den aus (28) entnommenen Verschiebungen {, und £, an Stelle von (33) die Gleichung 


2 mean) [84-422 h)e "+ 4] 
k=1 


= 4 ren [e 4" Zn (3—49—2 2, h)e’«®] rer Sin in) Jo(A;r) (36). 


vo,h 3mrph+2(—v)r]) 
Da SE SE Eee 


Mit dieser Gleichung verfahren wir genau so wie mit (30): Mit rJ,(A,r) bzw. rdr multipliziert 
und integriert, erhalten wir 


ee | 


0 


47,(1— 2») x (37) 
RR, DmÄl—r) 
na © 
bzw. 
Be en vo,h _ 3mPR+(A— re] 
BA sg ren. (38). 


Nach Einsetzen von mr in (37) bekommen wir zusammen mit (22) u. (22a). das zur Ermitt- 
lung von p; und q, notwendige Gleichungssystem. Die Gleichung (38) liefert dann den Wert 


von B. 


6. Eine Bemerkung für den Fall, daß zwischen Platte und Unterlage ein 
„Abheben“ stattfindet 
Die hergeleiteten Lösungen haben zur Voraussetzung, daß zwischen Platte und Unterlage 
stetiger Kontakt besteht, was in der Wirklichkeit offenbar nur solange der Fall ist, wie zwischen 
ihnen Druckkräfte herrschen. Dies wird bei den in der Praxis vorkommenden Belastungsfällen 
meistens der Fall sein, um so mehr, da ein ‚Abheben‘ bei schweren Platten schon durch das 
Eigengewicht verhindert wird. (Eine Berücksichtigung des Eigengewichtes der Platte kann in 


“erster Näherung dadurch geschehen, daß man das Eigengewicht als gleichmäßige Last zu f(r) 


addiert; eine exakte Behandlung erfolgt in der nächsten Ziffer.) 

Ergibt nun die Rechnung, daß zwischen Platte und Unterlage Zugkräfte auftreten, so kann 
man auf iterativem Wege folgendermaßen zu einer der Wirklichkeit entsprechenden Lösung 
kommen: Neben den anderen Randbedingungen werden innerhalb des Druckbereiches_ für 
z—h die vorangehend errechneten Spannungen (mit geeigneten Ergänzungen zur Wahrung des 


- 
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Gleichgewichtes) vorgeschrieben und dann aus den Verschiebungen Z an der Stelle 2=h gemäß 
0,—=—Kt die neuen Druckspannungen o, errechnet, die für z—=h wieder der Platte vorgeschrie- 
ben werden, usw. Auf diese Weise kann man nach einigen Schritten zu einer guten Näherungs- 
lösung kommen. Man wird freilich von dieser mühsamen Methode kaum Gebrauch machen, _ 
da man in der Praxis solche Voraussetzungen schafft, daß ein Abheben gerade nicht stattfindet. 


7. Die Berücksichtigung des Eigengewichtes 
Die exakte Berücksichtigung der Massenkräfte, im Gegensatz zu der in der vorigen Ziffer 
erwähnten näherungsweisen, ist bei dicken, insbesondere Fundamentplatten von Bedeutung. 
(So ist z.B. bei einem Schornsteinfundament das Eigengewicht der Fundamentplatte von der 
Größenordnung der aufzunehmenden Last.) Wie dies geschehen kann, soll hier kurz angedeutet 
werden. 
Bezeichnet y das spezifische Gewicht des homogenen Plattenmaterials, so sind die Diffe- 
rentialgleichungen (9) durch 
00,0 ul BU 00, 
Or r Re 


zu ersetzen, so daß aus (7), (8) und (39) die Differentialgleichungen 
ee Yan). | 
1-2)4+24+ 0-0 | { 
| Sn a I ER (40) 


5 
+ +-+y=0 27 EUR ER NE (39) 


0 de 
(1—2»)(Ao = oe 0 
AN | 
hervorgehen. 2 
E In einer dieser bald folgenden Mitteilung, welche die achsensymmetrisch belastete dicke 
Kreisplatte auf elastischer Unterlage und bei anderen Stützungen unter Berücksichtigung des 
‚Eigengewichtes ausführlich behandelt, wird gezeigt, daß die den Differentialgleichungen (40) 


= genügenden und die Randbedingungen (1) bis (3) befriedigenden Lösungen durch 
E = 3 mia ar +2) + + [Br 2a ee] 
= u 
4 (41) 
R BRETT SIR ELLIC FEN ET ; 
kr +4%(1 u N 10. +») (+ vr), 
e =2\ m [IB — 49 — 22) e Er + er’ + [3 —4v+ 222)’ + e *] 
Jı(Aer) (42) 


41 20)C0j4; en > 
| + aeldre PyRgRENFra 20641») 
dargestellt werden. Mit diesen Verschiebungsfunktionen werden nun zunächst die Randbedin- 
gungen (4) und (5) erfüllt. 

Aus (4) und (5) ergeben sich die an Stelle von (18) und (21) tretenden Beziehungen 


1—» (1—2»)yh & + 
A el Be. + Imre ’r* _ Sin A,h) — gy (Ah ek — EinA;h) 


Be 420-2) 0x Cof;h| a z ES ß 42h) + ee] 
Wr a k=1 k 


Jo (Ar r) 
4, (1 — 29») 


(43) 

DE | — 49 — 24h) ER" + | +4(1— 29), Sinz,n! 

BEAHETER) 
46(1+9) f 


und N 


— K(Ah-+B) 


gImrhe" +6inA;h)-+ 95 (Ar eir' + EinA,h)—2c,(1 —2») Ein; h] e =. 
r » Ph 


De 1ap)yB] Sa | | | 
| + Sinne wenn 303 (44). 


— 95 (Ay he’ — SinA,h)-+2(1— 2») c,&ofA;h] er 


7: sit Mech. 
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Mit diesen-Gleichungen verfährt man genau so wie mit (18) und (21), d.h. man multipliziert 
sie mit rdr bzw. rJ,(A,r)dr und integriert zwischen Null und a, wodurch man B und die Koeffi- 


zienten ohne Erfüllung von (6) ermitteln kann. Der letzten Forderung wird man auch hier im 
Sinne des de Saint-Venantschen Prinzips gerecht. Die Ausführung im Einzelnen bedeutet 
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- a um so weniger neue Schwierigkeiten, 


unde verzichten wir hier auf nähere Einzelheiten. 


da (26) und (27) ungeändert bleiben; aus diesem 


Für die eingespannte Platte ist in (41) das letzte Polynom zweiten Grades, d.h. 


r 


= 


—.durch\ 


3 oa) (2? + v2?) 


el —2y),y22 
4G(1—») 


laufen wie gerade angedeutet wurde. 
Eingegangen am 16. Juli 1951. 


\ 


Runge-Kutta-Verfahren unter 
höherer Ableitungen. 


In neuerer Zeit hat man bei der numerischen Inte- 
gration gewöhnlicher Differentialgleichungen die Ge- 
nauigkeit dadurch gesteigert, daß man außer der mit 
der Differentialgleichung gegebenen Ableitung, z.B. 
y’=f(z,y) im Falle der Gleichung 1. Ordnung, noch 
die durch Weiterdifferenzieren gebildete nächst höhere 
‚Ableitung, z.B. y’ =f’(z,y) = g(z,y), oder noch 
weitere Ableitungen zur numerischen Rechnung heran- 
' zieht. Bei den Interpolationsverfahren führt dies u.a. 
. auf sog. Hermitesche Interpolationsformeln, die sich 
durch hohe Genauigkeit auszeichnen!). 
"us liegt nahe, den Gedanken der Verwendung 
höherer Ableitungen auch auf andere Integrations- 
verfahren auszudehnen, insbesondere auf das Runge-. 
NK Verfahren, das sich für daspraktische Rechnen 


hi. en Sa 


Verwendung 


- 


> 


en hat. Ein besonders einfacher Weg der 
ung höherer Ableitungen besteht hier 
n die weiterdifferenzierte Gleichung als 
lgleichung höherer Ordnung auffaßt und 


rmeln für Gleichungen »-ter Ordnung — 
eachtung einiger hier anfallender Besonder- 
behandelt?). Beim Übergang auf Glei- 
dnung findet hier, wie man weiß, 

j keitssteigerung statt, so- 
ung von 3. oder höherer 
rgang von de ır 2. Ord- 
enauigkei' rändert, 


t3). Dafür 
eidende 
iesem Falle 


ischen Inte 
Math. Mech. 30 
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otz einiger Nachteile immer wieder als höchst brauch- 


sprechenden verallgemeinerten Runge- 


gration von. % 


x undaus Ahneh die gesuchten Endwerte yı, 41» .:, Y 


Zu ersetzen, während in (42) das letzte mit zr behaftete Glied zu streichen ist; diese Verschie- 
bungsfunktionen erfüllen dann auch die Randbedingung (6a). Die sonstigen Rechnungen ver- 


$ 


Gegeben sei allgemein eine Differentialgleichung 
n-ter Ordnung von der Form 


yn) = f(a,yy,..., ym-1)) FE AL 


von der wir annehmen, daß f genügend oft nach den 
unabhängigen Veränderlichen differenzierbar sei. 


Durch p-maliges Weiterdifferenzieren entsteht daraus _ 


die Differentialgleichung (n + p)-ter Ordnung 


wPrr=rm) ay,y,...ya-y) 

= 9, yYy,..,yR—U)).. © (2), 

eine Gleichung, die in jedem Falle unabhängig ist von 
der zweithöchsten Ableitung y(n+2—1). Dafür aber 
vereinfacht sich, wie man weiß, die Runge-Kutta- 
Rechnung von vier auf drei Zeilen an den drei Argu- 
 mentstellen Er Be | 


ı7=%, 2U=%oH+h2, Sit h, ; 
Pr denen die entsprechenden Werte 
ET HL DD NE Da 
sowie die hiervon abhängigen Funktionswerte 
A ne" +P s u 
k=gi @Fp)! = 


r 


2 : ’ 2% 5 nr+P : £ 
NE (4) 
zugehören. Außer diesen drei Funktionswerten k; 


benötigt man noch am Schrittanfang 2=x,, also an 
der Stelle I, die p Funktionswerte 


ym)=f, yn+y=f,... ym+-) =fa) 
rer | ..(b), 


die hier nicht, wie beim normalen Runge-Kutta- 
Verfahren (n + p)-ter Ordnung als Anfangswerte vor- 
liegen bzw. sich nach einer Runge-Kutta-Vorschrift 
Er n Daten des vorangegangenen Schrittes er- 
sondern als Funktionswerte zu den gegebenen 
ngswerten &9, Yos.--, Yo(®—1) zu berechnen sind. 
ellen II, III werden sie nicht benötigt. 
lrei k-Werte ky, kır, kırr ermittelt, so er- 
n aus ihnen nach: der allgemeinen Runge- 
Vorschrift die Mittelwerte 


REICH “ W, . 0 „n—1) £ 


(n—1) 
T 


z x . Mech. 
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Rechenschema für Differentialgleichungen 1. Ord- 


4 
an der Stelle x, = 2, + h.als Anfangswerte des neuen 1 . 
nung, einmal weiterdifferenziert: 


Schrittes. 
Im Falleder einmal weiterdifferenzierten Gleichung 
1. Ordnung (p=1, rn =1) lautet der einzige hier k; 


A ei in ine 


zur Errechnung von y, dienende Mittelwert ; = y — nl=f = 
1 £ 
k= (kr +2kı); —_— 1 114 
h I To Yo 1 1 vo | kr 
der Funktionswert kyrı wird also gar nicht gebraucht. L 
Damit aber reduziert sich hier die gesamte Rechnung I | zo+Al2| o+5”%t7 kı kı E 
N 


auf nur zwei Zeilen /7 und II und erfordert dabei die 


Berechnung nur dreier Funktionswerte f,, ky, kır. Der | 2% th | Yı=YyHttıotk | | 3 
Fortfall erhöhter Genauigkeitsordnung ER Fr in A h| 
diesem Falle bei gleichem Arbeitsaufwand durch ent- : f en san Ene „ 
sprechend verkleinerte Schrittgröße wettmachen. MitVohwgen: 3 (#72 kız) F 
Für die beiden praktisch wichtigen Fälle der Diffe- j n 4 x 
rentialgleichung 1. und 2. Ordnung sei das Ergebnis Bechenschema für Differentialgleichungen 2. Ord- ö | 
nochmals zusammengestellt. nung, einmal weiterdifferenziert: ’ { 
ee ee Rt ee 
a h? n 
i | en | Yy v=y’h u=/Z k=-f7 
? 
I To Yo v0 Yo kı 
1 1 1 3 
II %+ 5) ytzautgargh Yottotzkı kır 
5F = II | th | Yo+!ıo+ %ao+ kır | v0 +2%+3 kır kırı 
F; \ 
= 
23 | z0+h | Yı=YyHtPottont+k | met 2itk | 
er: } 1 
I: Gleich 1. Ord 3 MEAE = 
= eichung rdnung Miltchrertäl P 309 kr+12kıı —kııı) 
33 y’ = f(®, y)- k’=krı+2kır. 


v 
Ei 


1. Mitnahme von f’: Verfahren für Differential- 


Rudolf Zurmühl. 
gleichung 2. Ordnung. 


Darmstadt. 
Zwei Zeilen für die Schrittrechnung, drei Funk- 
tionswerte. 

Keine Genauigkeitssteigerung, Genauigkeit von 
der Ordnung ht. 


2.Mitnahme von f(®” mit p>2: Verfahren 
(p + 1). Ordnung. 


Resonanzlösungen inhomogener Mathieuscher 
Systeme. 
1. Aufgabenstellung 

Bei Untersuchungen der Stabilität erzwungener 
Schwingungen [1] wird man auf ein inhomogenes 
Differentialgleichungssystem mit periodischen Koeffi- 


Drei Zeilen für die Schrittrechnung, p+3 Funk- zienten geführt: 

tionswerte. PR . N F H 

Genauigkeitssteigerung auf die Ordnung hr +3, Bin, 73 Dil tn PS BD 23 RESTE me 
Gleichung 2. Ordnung: ee QM). 


Hierin haben die periodischen Koeffizienten und die 
Störglieder die gleiche Periode 2x, da die Zeitvariable 
durch s=wt dimensionslos eingeführt worden ist. In 
(1) bezeichnen ®„ die Eigenfrequenzen, Fyyn, Hm sind 


y’=f(e, Y; y') . 
Mitnahme von f(P mit p>1: Verfahren (p-+2). 


Ordnung. Gleichungskoeffizienten, während ® und s<1 als 
Drei Zeilen für die Schrittrechnung, Parameter für Störfrequenz und Amplitude aufgefaßt 
p-+3 Funktionswerte, _ werden. Die Koeffizienten Fnn, Hm sind im allge- 


meinen Funktionen des Parameters . Man kann 
daher das System (1) in gewissem Sinne als eine Ver- 
allgemeinerung der inhomogenen Mathieuschen 
ERS betrachten und die Lösungs- 
methoden aus der Theorie dieser Gleichung auf E® 


Genauigkeit für y von Ordnung AP +# 
Genauigkeit für y’ von Ordnung AP +°. 


Es folgen die Rechenschemata für die einmal 
weiterdifferenzierte Differentialgleichung 1. und 2. 


Ordnung als Beispiele, die sich an Hand der allge- 
meinen Runge-Kutta-Schemata (vgl. Anm, 2) 
leicht beliebig fortsetzen lassen. Das zweite Schema 
ist beispielsweise auch gültig für die zweimal differen- 
zierte Gleichung 1. Ordnung, wobei lediglich die 
Wertefürv, nichtnach der Runge-Kutta-Voschrift, 
sondern als Funktionswerte berechnet werden, und 


zwar nur für die Zeile I, ind wobei der zweite Mittel. 


wert k’ entfällt. Es bedeuten 
h2 
v=yh, 2 


u=—Y 


System (1) übertragen. Da überdies die allgemeinen 
Eigenschaften der Integrale des homogenen Systems 
hinreichend bekannt sind [2], hat man nur noch die 
Partikularlösungen des vollständigen Systems zu 
berechnen. 


2. Eigenschaften der Lösungen 


Die allgemeinen Eigenschaften des partikulären 
Integrales der inhomogenen Mathieuschen Einzel- 
gleichung sind von Gertrud Kotowski [3] mit der 
Methode der Variation der Konstanten ausführlich _ 
untersucht worden. Auch hier ist dieser Lösungsweg 
zweckmäßig. ER 
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Die Gleichungen (1) haben für 


Hm=0 (m=1,2,..., N) 
die Lösungen 


2N 
Um(s)= % &nvmn(8) . (2), 
n=1 
wenn die vmn(s) 2N Teillösungen sind, die ein 
Fundamentalsystem bilden. Betrachtet man nun die 
&n (n=1,2,...,2N) als Funktionen von s und be- 
stimmt sie aus 


2N : 
en An(s) mn (8) =(, 


2N 
22 (8) umn(s) =e Hm cos 3, (mel, 2,0,.,-N), 


so kann man auch die Lösungen des vollständigen 
Systems angeben (vgl. z.B. [4], 8.61). Mit der 
Cramerschen Regel berechnen sich zunächst die 
Anls) aus: &,-A=4A,„, worin die Wronskische 
Determinante A = Det |vmn, ümn| und die Deter- 
minante 4, vorkommen. 4, entsteht aus A nach 
Ersatz der n-ten Spalte durch die Funktionen der 
rechten Seiten von (1). Die Nennerdeterminante A 
kann für keinen Wert von s verschwinden und ist hier 
überdies konstant, da (1) keine Ableitungen erster 
Ordnung enthält. Auch über die Zählerdeterminante 
A„ kann man einige Aussagen machen. Nach E. 
Mettler [2] kennt man den Aufbau der Lösungen des 
zu (1) gehörigen homogenen Systems im einzelnen: 


N 


m= 3 (ane!!n® pmn(s) 
n=1 


(3) 
T&y ne en! mn(—8)), (m=1;2,0,.25N) | 


Hierin sind die &, die 2N Integrationskonstanten, 
+ 0n und — on im allgemeinen 2N komplexe charak- 
teristische Exponenten und die 91n(s) periodische 
Funktionen mit der Periode 2x. Bildet man nun A,, 
so kann man aus (2 N—1) Spalten einen Exponential- 
faktor ausklammern und diese Faktoren heben sich 
bis auf einen weg. Deshalb läßt sich 4, auf die Form 


(4) 


bringen, worin entweder das positive oder negative 
Vorzeichen im Exponentialfaktor gilt. Entwickelt 
man nun die in (4) vorkommende Determinante nach 
den Elementen eH„coss, so bleibt eine (2 N—1) 
reihige Determinante mit gleichperiodischen Funk- 
tionen Ymn und 97mm übrig. Mithin ist (4) aus Sum- 
manden von der Form 


e =ien® ge Hm cos 8: By (8) (4) 


zusammengesetzt. Der Aufbau der periodischen 
Funktionen ©, in diesen Summanden interessiert 
dabei im einzelnen nicht. Es wird jetzt vorausgesetzt, 
daß die mn und ömn durch Fourierreihen dar- 
zustellen sind. Dann lassen sich die Summanden (4) 
in der Form 


53 e a d. an „et -149)) 
(4”) 


An=e "®n® Det |e HmCoS 8, Ymn Omn| 


& 

=E 5 an BR 
schreiben, worin die Werte der Fourier koeffizienten 
c, aus den Funktionen ®,„, im weiteren nicht bekannt 
zu sein brauchen. 

Da A= const. ist, hat man zur Berechnung der &, 
nur noch Quadraturen über die Determinanten A, 
auszuführen: 


nz [ Ant) ae 


J ee 
Hierin sind durch s, die Integrationskonstanten ein- 
. geführt. 


re - 


Dadurch erhält man wieder die gleichen 


Exponentialfunktionen wie in (4), wenn nicht die 
Resonanzbedingungen 


Er lim nE+NRe on F0, | 
(6) 


verletzt sind. Wenn die Ungleichungen (6) für jedes r 
und n bestehen, dann ist das partikuläre Integral des 
vollständigen Systems (1) reinperiodisch. Denn nach 
Einsetzen der eben berechneten Funktionen in (3), 
bekommen die a, gerade den zu (4°) jeweils gehören- 


den reziproken Exponentialfaktor e’?®n°, so daß 
diese Faktoren ganz aus der Partikularlösung ver- 
schwinden. Die so entstehende Lösung wird später 
auf einfachere Weise berechnet und diskutiert. — 


Der Resonanzfall mit nichtperiodischen Lösungen 
tritt zufolge Ungleichung (6) dann ein, wenn zugleich 
für eine oder mehrere spezielle Indexziffern 


Sm _n=0, | 


Neon=P, (P Z 0 ganz, Bl) m) 


ist. Da die charakteristischen Exponenten +0, und 
—_n (n=]1, 2,...., N) nur bis auf einen beliebigen 
ganzzahligen Realteil bestimmt sind (vgl. Gl. (3)), be- 
deutet die, Resonanzbedingung (7), daß in der Lösung 
(3) gleiche charakteristische Exponenten auftreten. 
In einem solchen Fall zusammenfallender Exponenten 
treten nach der Floquetschen Theorie im allge- 
meinen nichtperiodische Partikularlösungen mit Poly- 


nomen in s bis zur Potenz s*! in die Lösungen (3) 
ein, wenn u die Vielfachheit des charakteristischen’Ex- 
ponenten bezeichnet. In den Sonderfällen allerdings, 
in denen die Determinante zur Bestimmung der 0, 
einen kleineren Rang als 2N—1 hat, wird man nur auf 
Polynome niederer Ordnung geführt ([4], S. 238). 
Hier soll nur der praktisch wichtigste Fall des Be- 
stehens der Resonanzgleichung (7) für einen speziellen 
Zeigerwert n weiter verfolgt werden. Ein so spezieller 
Wert zweier charakteristischer Exponenten ist natur- 
gemäß nur für ganz bestimmte zusammengehörige 
Parameterwerte ® und e möglich und die Bedeutung 
der so auftretenden Kurven w(e) als Grenzkurven in 
einer Struttschen Karte wird nochmals erwähnt 
werden. Wenn also +o, und — og, in diesem Sinne 
zusammenfallen, ist u—=2 und ein säkulares Glied in 
der ersten Potenz zu erwarten. Tatsächlich liegt ein 
Ausnahmefall nicht vor, denn genau wie bei der 
Mathieu-Gleichung gibt'’es zu Parameterwerten w(e), 
die gemäß (7) auf eine periodische Lösung führen, keine 
weiteren linear unabhängigen periodischen Lösungen. 
Da nämlich für eine allgemeine, weder geradenoch un- 
gerade partikuläre Lösung @mn(s) (m=1,2,...,'N; 
n=const) auch mn (—s)9 Lösung ist, kann man 
daraus durch Summe oder Differenz eine gerade oder 
ungerade Lösung bilden und sich daher auf die Be- 
trachtung solcher Lösungstypen beschränken. Da 
weiter ein Lösungsvektor eindeutig die Anfangsbe- _ 
Omn(0) =0 erfüllen können muß, ist eine Fundamen- 
tallösung entweder gerade oder ungerade und zwei 
linear unabhängige gerade oder ungerade Lösungen 
gibt es nicht. Weiter kann man für kleine e durch 
eine einfache Entwicklung nach dem Störungspara- 
meter e leicht erkennen, daß zu einer geraden oder 
ungeraden Lösung jeweils verschiedene Parameter- 
kurven w(e) gehören. Demnach können zwei perio- 
dische Lösungen im allgemeinen nicht gleichzeitig 
auftreten und für die daher in (3) fehlende Partikular- 
lösung muß eine Säkularlösung mit der ersten Potenz 


"von s eintreten. 


Zur Berechnung der partikulären Lösung der voll- 
ständigen Gleichungen hat man über die Determinan- 
ten A, zu integrieren. Da nach Obigem diese Deter- 
minanten aber in einer Spalte lineare Potenzen von s 
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enthalten, treten in die Lösungen Summanden mit s 
und stein. Diese sind daher für unbegrenzt wachsende 
Zeit (s> oo) nicht mehr beschränkt und sollen daher 
als Resonanzlösungen bezeichnet werden. 

Die bisher gefundenen Ergebnisse lassen sich zu 
folgender Aussage zusammenfassen. Die partiku- 
lären Lösungen der vollständigen Gleichungen (1) sind 
im allgemeinen periodisch und damit beschränkt. 
Unbegrenzt wachsende Resonanzlösungen sind nur 
dann möglich, wenn die Parameter w, & spezielle 
Werte annehmen. Diese Resonanzlösungen werden 
mit der Zeitvariablen wie s” (n >2, ganz) unendlich. 
Es ist noch nachzuweisen, daß der Resonanzfall nur 
für solche Parameterwerte wirklich eintritt, für welche 
das homogene System eine gerade periodische 
Lösung zuläßt. \* 

Die im allgemeinen vorhandenen periodischen Par- 
tikularlösungen sollen jetzt einfacher mit einem 
Fourieransatz untersucht werden. Der Ansatz 


Um = I’ (Gm ,r 60818 + bm,rsin rs), (m=1,2,...,N) 
r=0 € 
.& 


wird in (1) eingeführt und liefert zwei voneinander 
unabhängige Gleichungssysteme für die Fourier- 
koeffizienten an,, und bn,r. Mit den Abkürzungen 


Om Finn Hm 
m= > Amn- po: WE 
ergeben sich für die a, ,: 
N) ) 
2 N 
gm Am,o + I) Amnan,.ı=0, 
= n=1 
| 
N 3 (9) 
(gqm— 1) am,ı te I Amn(an,a+2an,0)=E Ems 
n=1 
EI Are) 
R N 
(qm —r?) Am,rt& I) Amn (Un, r-ı + @n,r+ı)=0 
n=1 
und für die bn,, mit bm,=0: 
r=1 
B N 
(gam—1)bm,, te 3 Amnbn,e=0, 
n=1 
(10) 


N en 


N 
(am r?) bm,r + & I Amn(bn ‚r-ıtbn,r +1) =0. 


n=1 


Die Gleichungssysteme enthalten die Parameter w 
und e und abhängig von diesen lassen sich -die 
Fourierkoeffizienten etwa mit Hilfe unendlicher 
Determinanten berechnen. Das System (9) ist in- 
homogen im ea zu den homogenen Gleichun- 
en (10). Die ra ne verschiedenen Lös- 

rkeitsbedingungen sollen jetzt diskutiert werden. 


3. Lösbarkeitsbedingungen 


a) Homogene Gleichungen. Betrachtet man zu- 
‚ nächst das Differentialgleichungssystem (1) ohne die 
Störfunktionen und macht wieder den Lösungsansatz 
(8) für ganzperiodische Funktionen, so wird man auf 
die Koeffizientengleichungen (9) und (10) geführt, in 
denen die ®,=0 (m=1,2,..., N) zu setzen sind. 
Für N=1,d.h. für die Mathieu-Gleichung allein, 
gehen (8) und (9) in die bekannten Gleichungen für 
die Fourierkoeffizienten der Mathieuschen Funk- 
tionen Ce, und se, über. Die voneinander unab- 
hängigen, homogenen Gleichungssysteme (8) und (9) 
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sind nur dann nichttrivial lösbar, wennn die zuge- 
hörigen Koeffizientendeterminanten verschwinden: 


Ara) (w, e)=0 (%), 
Ard) (wo, e)=0 (10°). 


Es gibt also nur für ausgezeichnete im allgemeinen 
voneinander verschiedene Parameterwerte w, & ent- 
weder eine gerade oder eine ungerade periodische 
Lösung mit der Periode 2x. Damit eine derartige 
periodische Lösung bestehen kann, muß ein go, in der 
Lösung (3) genau den Resonanzbedingungen (7) ge- 
nügen. Da weiter Sm 0, =0 die allgemeine Bedingung 
für die Grenze zwischen stabilen und instabilen Lö- 
sungen ist, definieren die zugehörigen Werte w(e) die 
Grenzkurveeines Instabilitätsbereiches in der Strutt- 
schen Karte. Und zwar handelt es sich um die Grenzen 
der Instabilitätsbereiche erster Art, die für 2=0 von 
den Frequenzen &(e=0)=2w,,/p (p>0 ganz) aus- 
gehen [2]. Hier speziell sind es jedoch nur Bereiche 
mit gerader Ordnungszahl p. Denn die Determi- 
nantenbedingungen (9) und (10°) führen für &>0 auf 
die Beziehungen 


gm—r=0(m=1,2,..,N; r=1,2,..., ©), 
d.h. om=rw(e=0), 
wie sie nur in den Fußpunkten der Bereiche erster 


‚Art und gerader Ordnungszahl erfüllt ind. Von diesen 


Fußpunkten gehen zwei getrennte Grenzkurven aus, 
für deren Parameterwerte ®,s je eine gerade oder 
ungerade periodische Lösung mit der Periode 2 
bestehen. 

b) Vollständige Gleichungen. Sind die Parameter & 
und e so gewählt, daß (9) erfüllt ist, so kann man die 
Fourierkoeffizienten der Partikularlösung der voll- 
ständigen Differentialgleichungen nicht berechnen. 
Eine periodische Lösung besteht nicht, denn dieser 
Fall entspricht dem in Abschn. 2 untersuchten Reso- 
nanzfall. Wenn die Parameter ® und z hingegen 
anders und zwar so gewählt sind, daß (10) erfüllt ist, 
dann ist für die gleichen Parameterwerte A/a) (w, &)£0 
und die Fourierkoeffizienten der partikulären 
Lösung der vollständigen Gleichungen können be- 
rechnet werden. Durch diese Ergebnisse sind die 
Resonanzmöglichkeiten der partikulären Lösung ein- 
geschränkt. 

4. Zusammenfassung. Die partikulären Lö- 
sungen des vollständigen Differentialgleichungs- 
systems (1) sind im allgemeinen periodisch mit der 
Periode der Störglieder. Resonanzlösungen mit unbe- 
grenzt wachsender Amplitude (Wachstum sr, n> 2) 
gibt es nur für diejenigen Parameterwerte w, &, welche 
gerade, gleichperiodische Lösungen der homogenen 
Gleichungen zulassen. Dies ist nur auf jeweils einer 
der beiden Grenzkurven der Instabilitätsbereiche 
erster Art und gerader Ordnung möglich. Da die 
Grenzkurven von den Frequenzen 2w„/p (P=2, 4, 
6,...; m=1,2,...,N) ausgehen, kann man diese 
Resonanzfälle als Subresonanzen bezeichnen. Die in 
der Nähe der Schwebungsfrequenzen om +@n/p 
liegenden Bereiche zweiter haben in diesem Zu- 
sammenhang übrigens keine Bedeutung. 
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Über positiv definite Matrizen 

In letzter Zeitist in verschiedenen Arbeiten, die mit 
Fragen der Stabilität in Zusammenhang stehen, der 
Ausdruck „positiv definit‘‘ in bezug auf nichtsymme- 
trische, reelle gqüadratische Matrizen gebraucht wor- 
dent). Dabei wurde als kennzeichnend für diesen Be- 
griff bald die eine 
schaften 
M': sämtliche Hauptminoren sind positiv, 
R: sämtliche Eigenwerte haben positive Realteile, 
S: sämtliche Eigenwerte des symmetrischen Teils sind 


positiv 
ee gelegt und offenbar der Umstand außer 
acht gelassen, daß die drei durch M bzw. R bzw. S 
definierten Klassen von Matrizen zwar bei Beschrän- 
kung auf symmetrische Matrizen identisch, im allge- 
meinen aber verschieden sind. Die folgenden Zeilen 
sollen einer Klarstellung dienen. Es wird.darin gezeigt: 

I. Aus S folgt sowohl R als auch M. 

II. Im übrigen ist keine der drei Eigenschaften 
R, 8 eine Folge von einer der beiden andern. 


Auch M und R zusammen sind für die Eigen- 
schaft S nicht hinreichend. 


Will man die Bezeichnung ‚positiv definit‘“ bei 
nichtsymmetrischen reellen Matrizen überhaupt ver- 


‚ bald die andere der drei Eigen- 


M, 
III. 


bei symmetrischen Matrizen geläufigen Verhältnisse 
zweckmäßig, darunter die Matrizen mit der Eigen- 
schaft 8, also solche zu verstehen, deren symmetri- 
scher Teil positiv definit ist. 

Wir beweisen zunächst die Behaup 


ptung I 

Il. Aus Sfolgt R. Diesist eine einfache Folgerung 
aus dem Satz von I. Bendixson?): . Der reelle Teil 
jedes Eigenwertes einer reellen Matrix liegt in dem 
abgeschlossenen Intervall, das durch den kleinsten 


begrenzt wird. ; 
 I2. Aus S folgt M: Um dies durch vollständige 
Induktion zu beweisen, setzen wir voraus, daß die Be- 
hauptung für alle Matrizen, deren Ordnung unter n 
‘liegt, ricehtigist, und zeigen zunächst, daß sie dann 
auch für alle Matrizen n-ter Ordnung zutrifft. Es sei 
‚also die Matrix X von der n-ten Ordnung und der 
symmetrische Teil X, von X positiv definit. Ist aber 
eine symmetrische Matrix positiv definit, so gilt das- 
selbe für die Matrix jedes ihrer Hauptminoren. Folg- 


positiv definit. Nun sei ® die ®, entsprechende Unter- 
matrix von W. ®;ist dann der symmetrische Teil von 
8 und daher hat ® die Eigenschaft 8. Liegt also die 
Ordnung von ® unterhalb von n, so hat ® n. V. auch 
die Eigenschaft M, insbesondere ist also det ® >0. 


wegen Il die Eigenschaft R, nämlich nur Eigenwerte 
mit positivem Realteil. Somit ist auch das Produkt 
‚ller Eigenwerte von X positiv d.h. det U >0. 

Für'n—=1 ist die Aussage 12 trivial. (Es ist auch 
sicht, ihre Richti 


olglich gilt sie allgemein. 


ee 


folgenden Beispielen erkannt werden. 
l. Sist weder Folge von R noch von M. 


er ei 1. 5\ Hauptminoren: 1, 1, 6_ 
ee ) 14:8 


ices definies positives. Ann, 
(29. Parodi M., A propos 


129. Parod Ap 
Acad. Sci. Paris 232 (1951) 


; : xelles 

a ter “ st 

204206. 
Acta math. 
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wenden, so erscheint es demnach im Hinblick auf die. 


und den größten Eigenwert ihres symmetrischen Teils 


‚lich ist die Matrix jedes Hauptminors ®, von Us 


Hat aber ® die Ordnung n, so ist B=AX und Ahat 


gkeit für n = 2 direkt zu bestätigen.) 


_ Die Richtigkeit der Behauptungen II und III kann _ 
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II 2. M ist nicht Folge von R. 


-) 


—l 2 
113. Rist nicht Folge von M. (n>2) 
Hauptminoren: 1,1,1; 10,10, 1; 100 


= | 12:0 ) 
3(—1+ 1/99), 4. 


8.13 
—1 3 1 

Hier mußte n>2 genommen werden, da im Fall 

n=2 aus M die Eigenschaft ER folgt. 


Hauptminoren: —1, 2, 1 
Eigenwerte: 3 (+1 +i /3) 


Eigenwerte: 


III. 8 folgt nicht aus R und M. 
Dies ergibt sich aus dem zu II1 gegebenen Beispiel. 


Betrachten wir daraufhin noch kurz die von Herrn 
M. Parodi angegebenen Verfahren zur Bildung 
„positiv definiter‘‘ Matrizen! Es ist bekannt, daß 
jedes der beiden Systeme von Ungleichungen 


DS N 
ED IR de ec) 
k=1 
N 
Ai 22, |axi | RN ET (2) 
k=1 


(Fhi=l2,...,n) Y 


für die Eigenschaft R hinreichend ist.. Beide Systeme 
zusammen reichen für $ (und damit wegen I auch für 
M) hin, denn sie haben entsprechende Ungleichungen 
für den symmetrischen Teil 4, der Matrix:X = (a;x) 
zur Folge und daher ist X, positiv definit. 4 
In der zweiten der genannten Arbeiten stellt 
Parodi aus einer Matrix X, welche (1) und (2) ge- 
nügt, neue Matrizen her, indem er alle Diagonal- 
\ elemente um ein und dieselbe positive Größe ö ver-. 
; mindert, welche kleiner ist als der kleinste Eigenwert 
x, von X; Da die reellen Teile der Eigenwerte von X 
auch nicht kleiner sind als x,, so bleiben sie bei diesem 
Prozeß positiv. Die neuen Matrizen haben also die 
Eigenschaft R. Sie haben sogar die Eigenschaft &, 
‚denn bei dem beschriebenen Prozeß erfahren auch die 
Eigenwerte von W,; lediglich eine Verminderung um ö, 
_ bleiben also ebenfalls positiv. 


Die in derselben Note gemachte Aussage, Matrizen 
mit den Eigenschaften min «;; >maz |ajz | und det (a;x) 


+ I r 
>0 seien für n <3 ‚‚positiv definit‘“, muß, wenn 
man diesen Begriff durch 8 festlegt, auf n 2 einge- 


„ schränkt werden, wie das Beispiel 


192 62 N de 2132 
A-1612 lee 
EN ers ia), det X, = —111 zeigt. 


Für n =3 haben solche Matrizen noch ‘die beiden 
Eigenschaften M und R, für n>3 jedoch im allge- 
meinen keine derselben mehr. 

In der anderen Arbeit addiert Parodi zu der Matrix 
M,(a) mit a; =%>1 und a; = 1(i>£k) eine Matrix 
(cr), in der \ 
%—1 & -- n—1 

n +23 —3 


ist, und zeigt mittels eines Satzes von Ostrowski, 


elle — du (a) 


daß die so erzeugten Matrizen die Eigenschaft M 


haben, Dieselben Matrizen haben aber sogar die 
_ Eigenschaft 8 (und damit auch R), denn die Elemente: 
ihres symmetrischen Teils unterscheiden sich von den 
Elementen von Nl,(&) ebenfalls um weniger als d„ (a), 

-so daß auch der symmetrische Teil die Eigenschaft M 
besitzt. und demnach positiv definit ist. 


EEE: 


1 


ER N ee ne 2 = 
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Beide von Parodi angegebenen Verfahren zur Bil- 
dung „positiv definiter'' Matrizen führen somit auf 
Matrizen mit der Eigenschaft 8. Die Aufgabe, solche 
Matrizen herzustellen, wird aber auf die allgemeinste 
und zugleich einfachste Art dadurch gelöst, daß man 


zu einer positiv definiten symmetrischen Matrix, deren 
Herstellung in. bekannter Weise, etwa mittels quadra- 
tischer Formen, leicht geschehen kann, eine beliebige 
schiefsymmetrische Matrix addiert. 


Stuttgart. Erich) Schönhardt. 
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A.0Ostrowski, Vorlesungen über Differen- 
tial- und Integralreehnung. Band II: Differen- 
tialrechnung auf dem Gebiete mehrerer Variablen 
(Lehrbücher und Monographien aus dem Gebiete der 
exaktenWissenschaften, Mathematische Reihe Band 5). 
482 S. mit 55 Abb. Basel 1951. Verlag Birkhäuser. 
Preis brosch. 63,— Schw. Fr., geb 67,— Schw. Fr. 

Das auf drei Bände berechnete Werk ist eine Nie- 
derschrift der Vorlesungen, die der Verfasser seit vie- 
len Jahren regelmäßig an der Universität Basel hält. 
Er willdamit einmal dem vorwiegend an den Anwen- 
dungen interessierten Hörer die Beherrschung des 
mathematischen Kalküls vermitteln, ohne ihn mit 
Suktilitäten zu plagen, andererseits dem Mathema- 
tiker die vermittelten Kenntnisse ausreichend be- 
gründen und eine Darstellung geben, die vor dem wis- 
senschaftlichenGewissen verantwortet werden kann, 

Der erste im Jahre 1945 erschienene Band enthält 
alles das, worauf sich der theoretisch weniger interes- 
sierte Leser beschränken kann. (Aus dem Inhalt: Das 
Wesen der Mathematik, reelle Zahlen, Funktionsbe- 
griff, Grenzwerte, stetige Funktionen, das bestimmte 
Integral, Ableitungen, Technik des Differenzierens 
und Integrierens, Funktionendiskussion, Reihenent- 
wicklungen). Der Integralbegriff wird an die Spitze 
gestellt. - 

Dervorliegendezweite Band behandelt zunächst die 
unendlichen Mengen, Funktionen auf Mengen, unend- 
liche Folgen und Reihen, Ergänzungen zur Differen- 
tialrechnung und Anwendungen auf die Analysis, ins- 
besondere für Funktionen mehrerer Veränderlicher, 
numerische Rechenmethoden. ‘Der Rest des Buches 
ist Anwendungen auf die Differentialgeometrie ge- 
widmet (Anfangsgründe der Kurven- und Flächen- 
theorie). 

Der dritte noch nicht erschienene Band soll Erwei- 
terungen nach der Seite der Integralrechnung hin 
bringen. Differentialgleichungen und Variationsrech- 
nung sind nicht in den Plan des Werkes einbezogen. 

Dem bekannten Verfasser ist es ausgezeichnet ge- 
lungen, Verständlichkeit und Strenge, Originalität 
und fesselnde Darstellung zu vereinen. Durch eine ge- 
wisse Breite ist das Buch auch zum Selbststudium ge- 
eignet; gern läßt der Verfasser den Leser ‚‚einen Blick 
hinter die Kulissen tun‘ und gibt in dem ausführlich 
gehaltenen Text historische Daten sowie Bemerkun- 
gen über die gewählte Terminologie und Bezeichnun- 
gen. Viele Einzelheiten findensich, die sonst selten 
in den Lehrbüchern gebracht werden, so daß gleich- 
zeitig auch ein Nachschlagebuch für den kundigen 
Leser entstanden ist. Besonders muß auf die ungemein 
reichhaltige Aufgabensammlung jedes Bandes hin- 
gewiesen werden; die Aufgaben, die sich nur auf inner- 
mathematische Beispiele beziehen, sind größtenteils 
neu und ausgesprochen interessant. Viele Leser werden 
allerdings dankbar sein, wenn den schwierigeren Auf- 
gaben die Lösungen sowie Hinweise auf dei einzu- 
schlagenden Weg beigegeben würden, selbst wenn da- 
w eine weitere Vermehrung des Umfanges ent- 
8 ©. 


Stuttgart. A Günther Schulz. 
Hermann Hömberger, Einführung in die 
Elektrotechnik. 302 8. mit 277 Abb, Braun- 
Bo nweig/ Berlin Hamburg/ Kiel 1951. Georg Wester- 

„mann Verlag. Preis br. 19,60 DM, geb. 21,80 DM. 
„Dieses Buch ist für die Anfangssemester der Fach- 
richtung Elektrotechnik an den Ingenieurschulen be- 


stimmt. Es will die Grundlagen der Elektrizitätslehre 
vermitteln, die technischen Anwendungen sind weite- 
ren Bänden vorbehalten. Der Stoff gliedert sich nach 
einer Vorbetrachtung über Einheiten, Maßsysteme 
und Gleichungen in die 4 Kapitel: 1. Elektrischer 
Strom (behandelt die elektrischen Grundbegriffe bei 
linienhafter Strömung, dargelegt an festen, flüssigen 
und gasförmigen Körpern), 2. Elektrisches) Feld 
(elektrostaiisches Feld, stationäres Strömungsfeld, 
langsam veränderliches elektrisches Feld), 3. Magneti- 
sches Feld (magnetische Grundbegriffe, stationäres 
magnetisches Feld, langsam veränderliches magneti- 
sches Feld, 4. zusammengesetzter Wechselstromkreis 
(symbolische Rechnung und ihre Anwendung auf 
Stromkreise). 

Abweichend von den üblichen Behandlungs- 
methoden wird mit zeitlich sich ändernden Größen, 
insbesondere mit Sinusgrößen, unmittelbar von An- 
fang an gearbeitet, demzufolge z. B. der Drehstrom 
ei im Anschluß an die Stromverzweigungen ein- 

eführt. 

= Charakteristisch für das Buch sind die Exaktheit 
der Darstellung, die Erläuterung durch zahlreiche 
Bilder und die konsequente Verwendung von Größen- 
gleichungen. Als kleine Anregungen für spätere Auf- 
lagen seien angeführt: Es ist notwendig, die Defi- 
nition des Widerstandes und das Ohmsche Gesetz 
klar auseinander zu halten. Nicht alle Energie- 
umwandlungen sind umkehrbar (S.103). Ein Ein- 
führen der magnetischen Induktion als Flußdichte 
über das Induktionsgesetz würde sicher leichter faß- 
lich sein als über GI]. (280). 

Insgesamt vertritt der Referent die Auffassung, 
daß der gebrachte Stoff, den der Lernende selbst- 
verständlich beherrschen muß, sich leichter einprägen 
würde, wenn das Einheitliche in den 3 Gebieten (Lei- 
ter, Nichtleiter, magnetisches Gebiet) entsprechend 
den neueren Bestrebungen stärker herausgearbeitet 
würde (Strombegriff und Spannungsbegriff im er- 
weiterten Sinn). 


Dresden H. Schönfeld, 

Dr. phil. habil. &. Gebelein und Dr. med. H.-J. Heite, 
Statistische Urteilsbildung. Erläutert 
an Beispielen aus Medizin und Biologie. XVI+ 
192 S. mit 50 Abb. und 20 Beispielen. Berlin-Göttin- 
gen-Heidelberg 1951. Springer-Verlag. Preis kart. 
15,60 DM. 

Die Arbeitsgemeinschaft zwischen Mathematiker 
und Mediziner hat sich für das vorliegende Buch als 
sehr fruchtbringend erwiesen. Esist direkt auf die Be- 
dürfnisse eines Mediziners oder Biologen zugeschnitten, 
der gewillt oder gezwungen ist, sich mit Statistik zu 
beschäftigen, und wahrt dabei die nötige mathemati- 
tische Exaktheit. Die zugehörigen Beweise mußten 
allerdings weggelassen werden oder sind nur ange- 
deutet; Erläuterungen mehr mathematischen Cha- 
rakters sind,im“Kleindruck gegeben. Überhaupt ist 
die Theorie möglichst knapp gehalten, und das Schwer- 
met wird auf die 20 Beispiele gelegt, die z. T. mehr- 

ach herangezogen werden, um eine gewisse Freizügig- 
keit in der Wahl der Methoden zu demonstrieren. 
Allerdings wird der Leser dabei manchmal in Rat- 
losigkeit darüber entlassen, welche Methode nun die 
empfehlenswerteste und rationellste ist. 

„Kurze Inhaltsangabe: Nach einer elementaren Ein- 
führung in die Beärbeitung des statistischen Materials 
(mit Rechenschemata zur 
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schnitt, Streuung, Schiefe und Exzess) werden die 
Gaußsche Verteilung und eine Normalverteilung 
zweiter Art — vom Typ y = exp (— (In @)?) — be- 
handelt. Dabei ist von der Verwendung gleitender 
Durchschnitte weitester Gebrauch gemacht worden, 
z. T. unter Benutzung eigener Ergebnisse von Gebe- 
lein. Auf eine ausführliche und kritische Darstellung 
der Korrelationstheorie folgt die Besprechung der 
Prüffunktionen, wobei die %2- und die i-Verteilung 
durch die oben eingeführte Normalverteilung zweiter 
Art angenähert werden. Die Frage der Zahl der Frei- 
heitsgrade kann nicht unerwähnt bleiben aber auch 
nicht befriedigend aufgehellt werden. Dann werden 
die drei statistischen Schlüsse, wie sie von Gebe- 
leins „Zahl und Wirklichkeit‘ her bekannt sind, 
auseinandergesetzt und auf die Mutungsgrenzen ein- 
gegangen. Die Betrachtung sehr geringer Korre- 
lationen bildet den Schluß. 

Das Buch ist sehr geeignet, das Verständnis für die 
statistische Arbeitsweise in weiten Kreisen zu propa- 
gieren, und esist nur bedauerlich, daß sein stärkstes 
Aktivum, nämlich das medizinisch-biologische Ge- 
wand, der erwünschten Verbreitung gewisse Grenzen 
setzen dürfte. 


Dresden. G. Opitz. 


P. Frauenfelder (Professor am Technikum in Winter- 
thur) und Dr. P. Huber (o. Professor an der Universität 
Basel). EinführungindiePhysik. I. Band, 
492 S..mit 391 Abb. München/Basel 1951. Ernst- 
Reinhard-Verlag. Preis in Leinen 16,50 DM. 

Die Einführung in die Physik, von der der erste 

Band Mechanik, Hydrodynamik und Thermodynamik 
enthält, ist als Ergänzung zur Vorlesung über Ex- 
perimentalphysik gedacht. Allerdings wird ein Haupt- 
teil des Stoffes der Experimentalvorlesung nicht ge- 
bracbt: Versuche werden nicht beschrieben. Das 
Buch soll also keinesfalls etwa die Rolle einer Vor- 
lesungsnachschrift übernehmen. Vielmehr wird in 
systematischem Gang der Lehrstoff — von den Grund- 
prinzipien ausgehend — behandelt. Das Werk kann 
daher als Vorstufe zur theoretischen Physik bezeichnet 
werden. Die Methode ist diejenige der theoretischen 
Physik, jedoch werden nur einfachere Grundprobleme 
behandelt. Die verwendete Mathematik ist bis auf 
einige kleinere Ausnahmen reine Schulmathematik. 
Die Kapitel über die Kreiseleigenschaften oder über 
die Anwendung des ersten und zweiten Hauptsatzes 
könnten sehr gut auch in einem kurzen Lehrbuch der 
theoretischen Physik stehen. Solche Abschnitte, die 
auch kompliziertere mathematische Hilfsmittel ver- 
langen, sind jedoch in Kleindruck wiedergegeben. 
Das Buch dürfte also für jüngere Studenten der Physik 
brauchbar sein. Auch interessierte Studenten der 
Ingenieurwissenschaften oder der Naturwissenschaften 
‚werden es verwenden, wenn sie tiefer in die prinzipielle 
Zusammenhänge physikalischer Fragestellung ein- 
dringen wollen, aber einen Kursus in theoretischer 
Physik nicht absolvieren können. Durch’ die Ziel- 
setzung, die damit gekennzeichnet ist, unterscheidet 
sich das Buch von den übrigen modernen deutschen 
Lehrbüchern und füllt eine Lücke, die nur bei dem 
üblichen System des Physikunterrichts nicht deutlich 
in Erscheinung tritt. 

Sorgfalt und Klarheit der Darstellung verdienen 
Anerkennung. Druck und Ausstattung sind ausge- 
zeichnet. 

Dresden. RBecknagel. 

Prof. E. Lohr, Vektor- und Dyadenrechnung 
für Physiker und Techniker. 2. Aufl. mit 
‘einem Nachtrag. XV + 488 8. m. 38 Abb. Berlin 
1950. Walter de Gruyter & Co. Preis geb. 24,— DM. 
Nachdem in letzter Zeit eine Reihe von Werken 
über Vektorrechnung erschienen sind, ist jetzt auch 

-  erfreulicher Weise das seit längerer Zeit vergriffene 
_ Werk von Lohr wieder zu haben. Die vorliegende 


Pa 


zweite Auflage des vielgebrauchten Buches lag be- 
reits 1944 druckfertig vor. Sie ist ein unveränderter 
Abdruck der ersten, die in Bd. 20 (1940); 8. 182 dieser 
Zeitschrift besprochen wurde. Angehängt sind einige 
in sich abgeschlossene Ergänzungen, die zusammen 
einen Umfang von etwa 75 Druckseiten haben. Sie 
behandeln in oft durchaus eigener-Weise einige Ge- 
biete, deren Hinzufügung nach Ansicht des Verfassers 
wünschenswert war, nämlich: Sphärische Trigono- 
metrie und Vektoren; Derivationen in krummlinigen 
Koordinaten; ebene Vektorrechnung und Funktionen- 
theorie; Verwandte Ergänzungen zur Potentialtheorie 
im Raum; eine Bemerkung über Reihenentwicklung; 
Matrizenmechanik und Integralgleichungen; Grund- 
züge des mathematischen Formalismus der Diracschen 
Theorie des Elektrons. Das zur ersten Auflage ge- 
sagte gilt auch für die zweite. Ihr Erscheinen wird von 
vielen begrüßt werden. 


Dresden. Willers. 


Prof. Dr.-Ing. J.Pirlet, Statik der rahmen- 
artigen Tragwerke. VII + 168S. mit 80 Abb. 
in 5 Taf. Berlin-Göttingen-Heidelberg 1951. Springer- 
Verlag. Preis geb. 24,— DM. 


Der Spannungsnachweis zum durchlaufenden Rah- 
men und durchlaufenden Träger auf elastisch einge- 
spannten Stützen ist eine wichtige Aufgabe des kon- 
struktiven Ingenieurbaues und daher auch in der Li- 
teratur nach den allgemeinen Ansätzen der Elastizi- 
tätstheorie mit zahlreichen Rechenvorschriften wieder-. 
holt bearbeitet worden. Diese unterscheiden sich im 
wesentlichen nur nach den Ansprüchen, die an das Er- 
gebnis gestellt werden und nach dem Zeitaufwand, der 
für die Lösung zur Verfügung gestanden hat. Der 
Verfasser stützt sich bei seiner Arbeit auf den Betrag 
der elastischen Einspannung der Riegel und Pfosten- 
stäbe am Rahmenknoten, wählt also mit dem Ein- 
spannungsgrad denjenigen Begriff, derden Spannungs- 
zustand der einzelene Bauteile in einfacher Weise be- 
stimmt. Dieser Gedankengang ist ausführlich be- 
gründet, für die Aufgaben zu Rechenvorschriften des 
konstruktiven Ingenieurbaues zusammengefaßt, mit 
zahlreichen Tabellen für die Anwendumg vorbereitet 
und durch Rechenbeispiele erläutert worden. Auf 
diese Weise ist eine Arbeit entstanden, die dem Bau- 
ingenieur, vor allem auch dank der hervorragenden 
Ausstattung des Buches, ausgezeichnete Dienste lei- 
sten wird. 


Dresden. K. Beyer. 


H.Reichardt, Der Einfluß der wand- 
nahen Strömung auf den turbulen- 
ten Wärmeübergang. (Mitteilungen aus dem 
Max-Planck-Institut für Strömungsforschung Nr. 3.) 
63 S. mit 13 Abb. u. 9 Tab. Göttingen 1950. Selbst- 
verlag Max-Planck-Institut für Strömungsforschung. 
Preis br. 5,— DM. 


Ineiner früheren Arbeit (diese Zeitschr. Bd.20 (1940), 
S. 297) hatte der Verf. die Geschwindigkeitsverteilung 
einer turbulenten Kanalströmung in großer Wand- 
nähe gemessen und war daraus zu Gesetzen für das 
wandnahe Gebiet gekommen, in dem sowchl turbu- 
lente als auch zähe Reibung vorliegt. Eine Dreiteilung 
der Grenzschicht in einen vollturbulenten Bereich, 
einen laminar-turbulenten und einen rein laminaren 
erwies sich als zweckmäßig für die Berechnung des 
Wärmeübergangs — jedoch insofern unbefriedigend, 
als es keine feste Grenze für die turbulente Bewegung 
und damit auch keine rein laminare Schicht geben 
kann. 

In einer neueren Arbeit (diese Zeitschr. Bd. 31 


- (1951), 8.208—219) gelangt nun der Verf. zu einer 
"allgemeinen Formel für die turbulente Geschwindig- 


keitsverteilung, nach der die turbulente Reibung 
kontinuierlich bis zur Wand abklingt. Damit war die 
Voraussetzung für eine neuerliche Behandlung des 
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Wärmeübergangs ohne jene Schichtgrenzen gegeben, 
die in der vorliegenden Arbeit durchgeführt wurde, 
wobei ein Ansatz mit zwei, aus Strömungsmessungen 
gefolgerten Konstanten (nach der oben genannten 
zweiten Arbeit) und ein früher (obige erste Arbeit) 
bereits bewährtes Näherungsverfahren Verwendung 
fand. 

Die Ergebnisse sind im wesentlichen folgende: 

1. „„Zähe“‘ und ‚‚wärmeleitende Wandschicht‘‘ sind 
verschieden stark. Dabei ist erstere der dimensions- 
lose Wandabstand, in dem der molekulare und der 
turbulente Anteil der Schubspannung gleich ist, die 
letztere entsprechend für die Wärmestromanteile. 

2. Das Temperaturprofil in der Reibungsschicht 
hängt wesentlich von Pr ab. 

3. Eine auf die Maximalwerte der örtlichen Ge- 
schwindigkeit und Temperaturdifferenz in der Rei- 
bungsschicht bezogene ‚‚örtliche Wärmeübergangs- 
zahl‘‘ erweist sich als von der Querschnittsform der 
Reibungsschicht weitgehend unabhängig. 

Die wertvolle Arbeit wird auch dem Experimen- 
tator manche Anregung geben., 
Weilburg. L. Schiller. 

W.Frösel, Experimentelle Unter- 
suchung der kompressiblen Strö- 
mung an und in der Nähe einer ge- 
wölbten Wand. (Mitteilungen aus dem. Max- 


Planck-Institut für Strömungsforschung, Nr. 4). 
105 S. mit 132 Abb. u.4 Tab. Göttingen 1951. Selbst- 
verlag Max-Planck-Institut für Strömungsforschung. 
Preis br. 9,— DM. > 

Die theoretische Behandlung von Strömungen mit 
annähernd Schallgeschwindigkeit stößt im allge- 
meinen auf große Schwierigkeiten. Der Verfasser hat 
es daher unternommen, ein erhebliches experimentelles 
Material zu beschaffen, das, aus rd. 20 000 Einzel- 
messungen des Druckes bzw. der Geschwindigkeit be- 
stehend, in 129 Kurvenbildern niedergelegt und im 
einzelnen ausführlich diskutiert ist. Dazu kommen 
noch die zugehörigen Schlierenaufnahmen (nur bei- 
spielsweise wiedergegeben), als deren Auswertung 
Reibungsschicht und Verdichtungsstoß jeweils in den 
Kurvenbildern eingetragen sind. Die Messungen wur- 
den 1944/45 in dem alten Göttinger Unter-Überschall- 
Kanal durchgeführt und beziehen sich auf die Strö- 
mung an zwei Kreisabschnitten mit R=50 bzw. 
100 mm und f = 5 bzw. 10 mm, wobei die Kanalhöhe 
durch Einbau von Beilagen schrittweise verengt 
werden konnte — bis zur Erzeugung einer ‚‚symmetri- 
schen Düse“. 

Der reichhaltige Inhalt der Schrift dürfte sowohl 
als Anregung für theoretische Behandlung wie auch 
für praktische Auswertung wesentlichen Nutzen 
liefern. 


Weilburg. L. Schiller. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


H.Brandt, Über das quadratische 
Reziprozitätsgesetz. (Berichte über die 
Verhandlungen der Sächs. Akademie der Wissen- 
schaften zu Leipzig. Mathematisch-naturwissenschaft- 
liche Klasse. Bd.99, Heft1.) 178. Berlin 1951. 
Akademie-Verlag. 


Dr. A. U. Huggenberger, Talsperren-Meß- 
technik. Meßverfahren,Instrumente 
und Apparate für die Prüfung der 
Bauwerkein Massenbeton. VII + 1328. 
mit 168 Abb. Berlin-Göttingen-Heidelberg 1951. 
Springer-Verlag. Preis brosch. 22,50 DM. 


Prof. E. Lohr, Vektor- und Dyaden- 
Br für Physiker und Tech- 
niker. 2. Aufl. mit einem Nachtrag. XV + 488 8. 
mit 38 Abb. Berlin 1950. Walter de Gruyter & Co. 
Preis geb. 24,— DM. 


Dr. P. Zühlke (Honorarprof. a. d. Univ. Marburg), 
KonstruktioneninbegrenzterEbene 
(Mathematisch-physikalischeBibliothek Reihel Nr.11.) 
3. Aufl. 42 S. mit 65 Abb. Leipzig 1951. B. G. Teub- 
ner Verlagsgesellschaft. Preis kart. 2,10 DM. 


Pröf. Dr. M. Zacharias (wissenschaftl. Mitarbeiter 
der Deutschen Akademie der Wissenschaften zu 


Berlin, Einführung in die projektive 


Geometrie. (Mathematisch-physikalische Bi- 
bliothek Reihe 1 Nr. 6.) Vierte Aufl. 54 S. mit 21 Abb. 
Leipzig 1951. B. G. Teubner Verlagsgesellschaft. 
Preis kart. 2,— DM. 


Colloque de G&omßtrie differenti- 
elle, tenu & Louvain du 11. au 14. avril 1951. 
(Centre belge de recherches math&matiques.) 235 8. 
Liege et Paris 1951. Georges Thone, Masson & Cie. 
Preis brosch 350 belg. fr., 2450 franz. fr. 


NACHRICHTEN 


Oberrheinische Mathematiker-Zusammenkunft in 
Basel. In der Zeit vom ‘14. bis 16. Dezember 1951 
trafen sich die Mathematiker der Universitäten Basel, 
“Freiburg i. Br. und Straßburg in Basel zu ihrer ersten 
‚„Oberrheinischen Mathematiker-Zusammenkunft‘‘, 
einer Arbeitstagung, die nunmehr jährlich an wech- 
selndem Orte wiederholt werden soll. Die Baseler 
Tagung brachte in zehn Vorträgen der Herren 
Ostrowski und Spieß aus Basel, Chabauty, Deny und 
Ehresmann aus Straßburg, Barner, Bilharz, Görtler 
und Leichtweiß aus Freiburg, sowie von Whitney 
(Harvard Univ., zur Zeit Straßburg) wertvolle neue 
Forschungsergebnisse derreinen und derangewandten 
Mathematik und bot unter der großzügigen Gast- 


freundschaft der Baseler Mathematiker die erstrebte_ 


Gelegenheit zu wissenschaftlicher Aussprache und 
persönlichem Näherkommen, 


Freiburg i. Br. H.Görtler. 


Freiburgi. Br.: Am 11. 2. 52 verstarb im Alter 
von 77 Jahren der emer. a. o. Professor der ange- 
wandten Mathematik an. der Universität Freiburg 
i.Br., Dr. Ernst-August Ansel. 


Zuschrift an den Herausgeber 


Zu W. Haacke: Bemerkungen zur Stabi- 
lisierung eines physikalischen Pendels I. 
Z. angew. Math. Mech. Bd. 31 (1951), S. 161—169 
(zu S. 166 oben). oR | 

Die Lagen ® = Q und ®=arc cos m. Sind nie stabil 


da eg R<mg gefordert werden muß. Sonst istr<0 oder 


® komplex, was beides mechanisch sinnlos ist. Dieses 
Ergebnis wird nicht weiter benutzt. Der Verf. dankt 


Herrn v.d. Vooren, Amsterdam für diesen Hinweis. 
Braunschweig. 
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